Time-dependent simulation of Czochralski silicon crystal growth

We have developed a detailed mathematical model and numerical simulation tools based on the streamline upwind/Petrov-Galerkin (SUPG) finite element formulation for the Czochralski silicon crystal growth. In this paper we consider the mathematical modeling and numerical simulation of the time-dependent melt flow and temperature field in a rotationally symmetric crystal growth environment. Heat inside the Czochralski furnace is transferred by conduction, convection and radiation, Radiating surfaces are assumed to be opaque, diffuse and gray. Hence the radiative heat exchange can be modeled with a non-local boundary condition on the radiating part of the surface. The position of the crystal-melt interface is solved by the enthalpy method. The melt flow is assumed to be laminar and governed by the cylindrically symmetric and incompressible Navier-Stokes equations coupled with the calculation of temperature.