The heat in Czochralski method is transfered by all transport mechanisms such as convection, conduction and radiation and convection is caused by the temperature difference in the molden pool, the rotations of crystal or crucible and the difference of surface tension. This study delvelops the simulation model of Czochralski growth by using the finite difference method with fixed grids combined with new latent heat treatment model. The radiative heat transfer occured in the surfce of the system is treated by calculating the view factors among surface elements. The model shows that the flow is turbulent, therefore, turbulent modeling must be used to simulate the transport phenomena in the real system applied to 8" Si single crystal growth process. The effects of a cusp magnetic field imposed on the Czochralski silicon melt are studied by numerical analysis. The cusp magnetic field reduces the natural and forced convection due to the rotation of crystal and crucible very effectively. It is shown that the oxygen concentration distribution on the melt/crystal interface is sensitively controlled by the change of the magnetic field intensity. This provides an interesting way to tune the desired O concentration in the crystal during the crystal growing.