Numerical investigation of heat transport and fluid flow during the seeding process of oxide Czochralski crystal growth Part 1: non‐rotating seed

For the seeding process of oxide Czochralski crystal growth, the flow and temperature field of the system as well as the seed‐melt interface shape have been studied numerically using the finite element method. The configuration usually used initially in a real Czochralski crystal growth process consists of a crucible, active afterheater, induction coil with two parts, insulation, melt, gas and non‐rotating seed crystal. At first the volumetric distribution of heat inside the metal crucible and afterheater inducted by the RF coil was calculated. Using this heat source the fluid flow and temperature field were determined in the whole system. We have considered two cases with respect to the seed position: (1) before and (2) after seed touch with the melt. It was observed that in the case of no seed rotation (ω_seed = 0), the flow pattern in the bulk melt consists of a single circulation of a slow moving fluid. In the gas domain, there are different types of flow motion related to different positions of the seed crystal. In the case of touched seed, the seed‐melt interface has a deep conic shape towards the melt. It was shown that an active afterheater and its location with respect to the crucible, influences markedly the temperature and flow field of the gas phase in the system and partly in the melt. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of active afterheater on the induction heating process in oxide Czochralski systems

Different shapes and orientations of an active afterheater for oxide Czochralski crystal growth systems are considered and corresponding results of electromagnetic field and volumetric heat generation have been computed using a finite element method (Flex‐PDE package). For the calculations, the eddy current in the induction coil (i.e. the self‐inductance effect) has been taken into account. The calculation results show the importance of an active afterheater, its shape as well as its geometry and position with respect to the crucible on the heat generation distribution in a CZ growth system. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of thermal and mass transport during Czochralski crystal growth of sapphire

The thermal and flow transport in an inductively heated Czochralski crystal growth furnace during a crystal growth process is investigated numerically. The temperature and flow fields inside the furnace, coupled with the heat generation in the iridium crucible induced by the electromagnetic field generated by the RF coil, are computed. The results indicate that for an RF coil fixed in position during the growth process, although the maximum value of the magnetic, temperature and velocity fields decrease, the convexity of the crystal‐melt interface increases for longer crystal growth lengths. The convexity of the crystal‐melt interface and the power consumption can be reduced by adjusting the relative position between the crucible and the induction coil during growth. (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Controlling the axial temperature gradient in inductively heated Czochralski systems

A heating system is described that allows to increase the temperature gradients in rf‐heated Czochralski setups during the growth. This system applies an active afterheater, operated by a separate induction coil parallel to that heating the crucible. By changing the inductivity of an additional coil located outside the growth chamber, the ratio of the rf currents flowing through the crucible and afterheater and, consequently, the geometry of the heat input to the setup can be altered. The efficiency of the heating systems has been demonstrated in various experiments.     

Simulation of temperature and flow fields in an inductively heated melt growth system

The goal of the research presented here is to apply a global analysis of an inductively heated Czochralski furnace for a real sapphire crystal growth system and predict the characteristics of the temperature and flow fields in the system. To do it, for the beginning stage of a sapphire growth process, influence of melt and gas convection combined with radiative heat transfer on the temperature field of the system and the crystal‐melt interface have been studied numerically using the steady state two‐dimensional finite element method. For radiative heat transfer, internal radiation through the grown crystal and surface to surface radiation for the exposed surfaces have been taken into account. The numerical results demonstrate that there are a powerful vortex which arises from the natural convection in the melt and a strong and large vortex that flows upwards along the afterheater side wall and downwards along the seed and crystal sides in the gas part. In addition, a wavy shape has been observed for the crystal‐melt interface with a deflection towards the melt. (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three‐dimensional numerical investigation of the effects of an open window on the convective heat transfer of an oxide Czochralski system

A three‐dimensional numerical analysis was carried out for a real Czochralski crystal growth furnace containing only gas and without any melt and crystal in order to investigate the effects of a small observation window on the temperature and flow field of the system. For this approach, the induction heating equations, the Navier‐Stokes equation with Boussinesq approximation, the continuity and energy equations have been solved in cylindrical coordinates using the finite element method. It has been found that the flow and thermal fields in the system are obviously three‐dimensional and non‐axisymmetric. The gas enters the system through the window is directed towards the opposite side wall where it is divided into two parts of vertical direction as well as expands in horizontal direction. Consequently, there is a spiral gas flow in the crucible and afterheater which rotates upwards in azimuthal direction along the walls. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A modified heater system for RF-Czochralski equipments

The temperature gradients in a Czochralski-type crystal growth equipment with induction heating can be advantageously influenced by additionally heating the growing crystal and the crucible bottom during the pulling procedure. An experimental solution is presented for both afterheater and bottom heater coupled to the same induction coil as the crucible.     

Numerical investigation of heat transport and fluid flow during the seeding process of oxide Czochralski crystal growth Part 2: rotating seed

In this paper, the role of seed rotation on the characteristics of the two‐dimensional temperature and flow field in the oxide Czochralski crystal growth system has been studied numerically for the seeding process. Based on the finite element method, a set of two‐dimensional quasi‐steady state numerical simulations were carried out to analyze the seed‐melt interface shape and heat transfer mechanism in a Czochralski furnace with different seed rotation rates: ω_seed = 5‐30 rpm. The results presented here demonstrate the important role played by the seed rotation for influencing the shape of the seed‐melt interface during the seeding process. The seed‐melt interface shape is quite sensitive to the convective heat transfer in the melt and gaseous domain. When the local flow close to the seed‐melt interface is formed mainly due to the natural convection and the Marangoni effect, the interface becomes convex towards the melt. When the local flow under the seed‐melt interface is of forced convection flow type (seed rotation), the interface becomes more concave towards the melt as the seed rotation rate (ω_seed) is increased. A linear variation of the interface deflection with respect to the seed rotation rate has been found, too. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two and three‐dimensional numerical modeling of induction heating in oxide Czochralski systems

In this article we report on a set of two‐dimensional and three‐dimensional numerical calculations for three different oxide Czochralski configurations in order to compare the results of the electromagnetic fields and the heat generation distribution. Two configurations without and with a gap between the crucible and active afterheater have axisymmetric conditions while the configuration with an open observation window in the afterheater is characterized by a non‐axisymmetric geometry. It has been found that in the non‐axisymmetric configuration under the influence of the observation window is located at the crucible side wall after a short distance behind the window cut out. Besides this influence the volumetric distribution of heat generation in the system is about symmetric. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The role of inner and internal radiation on the melt growth of sapphire crystal

In this paper, for an inductively heated Czochralski furnace used to grow sapphire single crystal, influence of the inner (wall‐to‐wall) and crystal internal (bulk) radiation on the characteristics of the growth process such as temperature and flow fields, structure of heat transfer and crystal‐melt interface has been studied numerically using the 2D quasi‐steady state finite element method. The obtained results of global analysis demonstrate a strong dependence of thermal field, heat transport structure and crystal‐melt interface on both types of radiative heat transfer within the growth furnace.     

Numerical optimization of czochralski sapphire single crystal growth using orthogonal design method

Crystal quality during Czochralski (Cz) growth is influenced significantly by the convexity of solid/liquid (S/L) interface, which is related to operating conditions, such as Radio‐Frequency (RF) coil position, crystal rotation and crucible rotation. Generally, a flat interface shape is preferred for high‐quality crystal growth. It is difficult to achieve the optimized conditions even from numerical modeling due to the large computational load from examining all of the affecting factors. Orthogonal design/test method, fortunately, provides an efficient way to organize the study of multiple factors with the minimization of computational load. In the paper, this method is adopted to investigate the affecting factors of Cz‐sapphire single crystal growth based on the coupled calculation of thermal field and melt flows. The orthogonal analysis clearly reveals the optimized growth conditions to achieve a relative flat S/L interface under the current ranges of the parameters.     

Global simulation of a Czochralski furnace for silicon crystal growth against the assumed thermophysical properties

In order to understand the effects of the thermophysical properties of the melt on the transport phenomena in the Czochralski (Cz) furnace for the single crystal growth of silicon, a set of global analyses of momentum, heat and mass transfer in small Cz furnace (crucible diameter: 7.2 cm, crystal diameter: 3.5 cm, operated in a 10 Torr argon flow environment) was carried out using the finite‐element method. The global analysis assumed a pseudosteady axisymmetric state with laminar flow. The results show that different thermophysical properties will bring different variations of the heater power, the deflection of the melt/crystal interface, the axial temperature gradient in the crystal on the center of the melt/crystal interface and the average oxygen concentration along the melt/crystal interface. The application of the axial magnetic field is insensitive to this effect. This analysis reveals the importance of the determination of the thermophysical property in numerical simulation. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Combined effects of crucible geometry and Marangoni convection on silicon Czochralski crystal growth

In order to understand the influence of crucible geometry combined with natural convection and Marangoni convection on melt flow pattern, temperature and pressure fields in silicon Czochralski crystal growth process, a set of numerical simulations was conducted. We carry out calculation enable us to determine temperature, pressure and velocity fields in function of Grashof and Marangoni numbers. The essential results show that the hemispherical geometry of crucible seems to be adapted for the growth of a good quality crystal and the pressure field is strongly affected by natural and Marangoni convection and it is more sensitive than temperature. (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Simulation of Temperature and Flow Field in the Melt for the Vapour‐pressure‐ controlled Czochralski Growth of GaAs

The influence of the melt flow on the temperature field and interface during the vapour‐pressure‐controlled growth of GaAs was studied numerically with the commercial general‐purpose program FIDAP^TM. The thermal boundary conditions for the domain of seed, crystal, boron oxide and crucible were taken from a global calculation for an equipment used at the IKZ to grow 6″ crystals. Due to the large melt volume the buoyancy forces become rather strong and have to be counteracted by reasonable rotation rates. Preliminary results have been obtained for iso‐ and counter‐rotation showing that the flow field exhibits structures on small scales. High rotation rates are needed to counteract the buoyancy flow efficiently and to achieve a smooth flat interface. Even if the the flow structure is not resolved in detail, the interface shape can be deduced form the calculations.     

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.     

Numerical investigation of the effect of heat shield shape on the oxygen impurity distribution at the crystal–melt interface during the process of Czochralski silicon crystal growth

In this study, a numerical simulation is performed to investigate the effect of the shape of the heat shield on the oxygen concentration in the melt. The results show that the oxygen concentration in the melt can be significantly decreased by increasing the speed of the argon gas near the crucible wall. This can be achieved by enlarging the horizontal length of the heat shield. The oxygen concentration at the melt–crystal interface varies with the length of the crystal growth. In the initial stage, there is a significant decrease in the oxygen concentration as the growth length increases. There is also a significant reduction in the emission of oxygen from the crucible wall due to the lower melt depth and crucible temperature. The transportation of oxygen impurity towards the melt–crystal interface is suppressed by the vortex motion in the melt. When the crystal exceeds a certain length, the oxygen concentration in the melt–crystal interface starts to increase with increasing crystal length, due to the drop in vortex motion in the melt.     

Influence of crucible geometry and position on the induction heating process in crystal growth systems

A series of 2D finite element numerical simulations of induction heating process for an oxide Czochralski crystal growth system has been done for different shapes and locations of a metal crucible. Comparison between the computational results shows the importance of crucible shape, geometry and its position with respect to the RF-coil on the electromagnetic field and heat generation distribution in the growth setup.     

3D simulation of the effects of growth parameters on the growth of sapphire crystals using heat exchanger method

3D simulations using the commercial CFDRC and FIDAP code, which are based on finite element techniques, were performed to investigate the effects of anisotropic conductivity on the convexity of the melt–crystal interface and the hot spots of sapphire crystal in a heat‐exchanger‐method crystal growth system. The convection boundary conditions of both the energy input to the crucible by the radiation as well as convection inside the furnace and the energy output through the heat exchanger are modeled. The cross‐sectional flow pattern and the shape of the melt–crystal interface are confirmed by comparing the 3‐D modeling results with previous 2D simulation results. In the 3D model, the “hot spots” in the corners of the crucible are donut shaped, and the shape changes with the value of the conductivity of anisotropic crystal. The outline of the crystal becomes more convex as the conductivity in the z direction (k_sz) increases. The outline of melt–crystal interface is elliptical when the anisotropic conductivity is moving in the radial direction (k_sx and k_sy). The portion at the outline touching the bottom of the crucible is smaller than the maximum outline of the crystal, meaning that the shape at the “hot spot”, changes with the value of the conductivities of anisotropic crystal. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three dimensional numerical study of the effect of large Grashof number on HEM crystal growth

Heat transfer and fluid flow in HEM crystal growth of silicon in cylindrical cavity is studied numerically. The walls of the crucible are heated to a fixed temperature. The exchanger that causes and induces natural convection is seated at the middle‐bottom of the crucible. The finite‐volume method is employed to solve the governing equations with proper boundary conditions. The effects of transport mechanism on the temperature distribution, melt flow, pressure and stream function are presented. We focus our work on the pressure field which has not yet been studied in HEM crucible. Also, we extend our work on a wide range Grashof number and for large numbers until 10¹² not yet studied in HEM furnace. It is found that the onset of flow fluctuations appears at Gr = 10¹⁰. Uniform temperature is observed in the entire melt at high Grashof number with development of a thermal boundary layer close to the exchanger. The thermal boundary layer thickness is calculated for strong buoyancy regime. Besides, for very high Gr number, buoyancy has less effect on temperature and then on melt‐crystal interface shape. During enlarging Gr, pressure evolution is related to temperature variation more than flow pattern.     

数值模拟自然对流对直拉单晶硅的影响

在直拉单晶硅生长的过程中,自然对流对晶体界面的形状、温度场及应力分布影响很大。本文采用二维模型对熔体内自然对流对单晶硅的影响作了数值模拟,在低雷诺数时采用层流模型,高雷诺数时采用紊流模型,Gr的变化范围从3×106到3×1010,这样涵盖了从小尺寸到大尺寸的直拉单晶硅生长系统。数值结果表明熔体的流动状态不仅与熔体的Gr有关,还与熔体高度和坩埚半径的比值密切相关。当Gr>108时,熔体内确实存在紊流现象,层流模型不再适合,随着Gr的增大,紊流现象加剧,轴心处的等温线变得更为陡峭,不利于晶体生长。