Effect of thermocapillary convection in an industrial Czochralski crucible: numerical simulation

Abstract The influence of thermocapillary or Marangoni convection on the growth of silicon crystals is investigated in an industrial Czochralski crucible using a quasi direct numerical simulations approach. An optimized parallel-vector block-structured finite volume Navier–Stokes equations solver is extended to include the effects of thermocapillary convection. Owing to the presence of surface tension gradients, the magnitude of the radial velocity towards the crystal becomes high at the free surface. Consequently, the temperature along the free surface is increased. The thermocapillary convection reduces the temperature fluctuations below the crystal and the free surface of the melt. However, below the crystal, the turbulent kinetic energy, k , is higher at different melt depths when the thermocapillary effect is taken into account. It is necessary to study the above-mentioned effects of surface tension driven convection since they influence the quality of the crystal.

[1]  Biswajit Basu,et al.  Three-dimensional simulation of flow and thermal field in a Czochralski melt using a block-structured finite-volume method , 2000 .

[2]  K. Yi,et al.  Spoke patterns on molten silicon in Czochralski system , 1994 .

[3]  K. Kakimoto,et al.  Natural and forced convection of molten silicon during Czochralski single crystal growth , 1989 .

[4]  D. Schwabe,et al.  Marangoni convection in open boat and crucible , 1981 .

[5]  K. Kakimoto,et al.  Oxygen distribution at a solid–liquid interface of silicon under transverse magnetic fields , 2000 .

[6]  W. Rhim,et al.  Thermophysical properties measurement of molten silicon by high-temperature electrostatic levitator: density, volume expansion, specific heat capacity, emissivity, surface tension and viscosity , 2000 .

[7]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[8]  D. Schwabe,et al.  Experiments on buoyant, thermocapillary, and forced convection in Czochralski configuration , 1983 .

[9]  D. Schwabe,et al.  Convection in a Czochralski crucible – Part 1: non-rotating crystal , 2001 .

[10]  B. Basu,et al.  Numerical Simulation of Fluid Flow and Heat Transfer in an Industrial Czochralski Melt Using a Parallel-Vector Supercomputer , 2000 .

[11]  Vishwanath Prasad,et al.  Turbulent transport of oxygen in the Czochralski growth of large silicon crystals , 1999 .

[12]  Adrian D. Jones.  Scaling analysis of the flow of a low Prandtl number Czochralski melt , 1988 .

[13]  W. Wilcox,et al.  Inhomogeneities due to thermocapillary flow in floating zone melting , 1975 .

[14]  Michael Schäfer,et al.  Efficient Simulation of Incompressible Viscous Flows on Parallel Computers , 1996 .

[15]  T. Fukuda,et al.  Observation of crystal-melt interface shape in simulated Czochralski method with model fluid , 1993 .

[16]  N. Imaishi,et al.  Three-dimensional numerical simulation of oxide melt flow in Czochralski configuration , 2000 .

[17]  M. Eguchi,et al.  The role of surface-tension-driven flow in the formation of a surface pattern on a Czochralski silicon melt , 2001 .

[18]  K. Yi,et al.  Structure of temperature and velocity fields in the Si melt of a Czochralksi crystal growth system , 1995 .

[19]  Brian Cabral,et al.  Imaging vector fields using line integral convolution , 1993, SIGGRAPH.

[20]  N. Imaishi,et al.  Three-dimensional numerical simulation of spoke pattern in oxide melt , 1999 .

[21]  Nobuyuki Kobayashi,et al.  Steady state flow in a Czochralski crucible , 1995 .

[22]  Robert A. Brown,et al.  Comparison of three turbulence models for simulation of melt convection in Czochralski crystal growth of silicon , 1999 .

[23]  J. Lumley,et al.  Instabilities, transition and turbulence in the Czochralski crystal melt , 1992 .

[24]  Georg Müller,et al.  Effects of various magnetic field configurations on temperature distributions in Czochralski silicon melts , 2001 .

[25]  T. Tsukada,et al.  Global simulation of a silicon Czochralski furnace , 2002 .

[26]  H. Wilke,et al.  Influence of Marangoni convection on the flow pattern in the melt during growth of Y3Al5O12 single crystals by the Czochralski method , 2000 .

[27]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[28]  S. Ostrach,et al.  Convective effects in crystals grown from melt , 1981 .

[29]  D. Schwabe,et al.  Coupling and separation of buoyant and thermocapillary convection , 1989 .

[30]  H. L. Stone.  ITERATIVE SOLUTION OF IMPLICIT APPROXIMATIONS OF MULTIDIMENSIONAL PARTIAL DIFFERENTIAL EQUATIONS , 1968 .