Three dimensional simulation of melt flow in Czochralski crystal growth with steady magnetic fields

Three-dimensional transient numerical simulations were carried out to investigate the melt convection and temperature fluctuations within an industrial Czochralski crucible. To study the magnetic damping effects on the growth process, a vertical magnetic field and a cusp magnetic field were considered. Due to our special interest in the melt convection, only local simulation was conducted. The melt flow was calculated by large-eddy simulation (LES) and the magnetic forces were implemented in the CFD code by solving a set of user-defined scalar (UDS) functions. In the absence of magnetic fields, the numerical results show that the buoyant plumes rise from the crucible to the free surface and the crystal–melt interface, which indicates that the heat and mass transfer phenomena in Si melt can be characterized by the turbulent flow patterns. In the presence of a vertical magnetic field, the temperature fluctuations in the melt are significantly damped, with the buoyant plumes forming regular cylindrical geometries. The cusp magnetic field could also markedly reduce the temperature fluctuations, but the buoyant plumes would break into smaller vortical structures, which gather around the crystal as well as in the center of the crucible bottom. With the present crucible configurations, it is found that the vertical magnetic field with an intensity of 128 mT can damp the temperature fluctuations more effectively than the 40 mT cusp magnetic field, especially in the region near the growing crystal.