Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Growth of Silicon Carbide Bulk Crystals by Physical Vapor Transport Method and Modeling Efforts in the Process Optimization

Silicon carbide bulk crystals were grown in an induction-heating furnace using the physical vapor transport method. Crystal growth modeling was performed to obtain the required inert gas pressure and temperatures for sufficiently large growth rates. The SiC crystals were expanded by designing a growth chamber having a positive temperature gradient along the growth interface. The obtained 6H-SiC crystals were cut into wafers and characterized by Raman scattering spectroscopy and X-ray diffraction, and the results showed that most parts of the crystals had good crystallographic structures.

Nonlinear Dynamic Responses of the Tensioned Tether under Parametric Excitations

The nonlinear dynamic responses of the tensioned tether subjected to combined surge and heave motions of floating platform are investigated using 2-D nonlinear beam model. It is shown that if the transverse-axial coupling of nonlinear beam model and the combined surge-heave motions of platform are considered, the governing equation is not Mathieu equation any more, it becomes nonlinear Hill equation. The Hill stability chart is obtained by using the Hill's infinite determinant and harmonic balance method. A parameter M, which is the function of tether length, the surge and heave amplitude of platform, is defined. The Hill stability chart is obviously different from Mathieu stability chart which is the specific case as M=0. Some case studies are performed by employing linear and nonlinear beam model respectively. It can be found that the results differences between nonlinear and linear model are apparent.