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-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

High-Strength Composite Fibers: Realizing True Potential of Carbon Nanotubes in Polymer Matrix through Continuous Reticulate Architecture and Molecular Level Couplings

Carbon nanotubes have unprecedented mechanical properties as defect-free nanoscale building blocks, but their potential has not been fully realized in composite materials due to weakness at the interfaces. Here we demonstrate that through load-transfer-favored three-dimensional architecture and molecular level couplings with polymer chains, true potential of CNTs can be realized in composites as Initially envisioned. Composite fibers with reticulate nanotube architectures show order of magnitude improvement in strength compared to randomly dispersed short CNT reinforced composites reported before. The molecular level couplings between nanotubes and polymer chains results in drastic differences in the properties of thermoset and thermoplastic composite fibers, which indicate that conventional macroscopic composite theory falls to explain the overall hybrid behavior at nanoscale.

Experimental Study Of Multi-Hole Cooling For Integrally-Woven, Ceramic Matrix Composite Walls For Gas Turbine Applications

In this paper, multi-hole cooling is studied for an oxide/oxide ceramic specimen with normal injection holes and for a SiC/SiC ceramic specimen with oblique injection holes. A special purpose heat transfer tunnel was designed and built, which can provide a wide range of Reynolds numbers (10(5)similar to 10(7)) and a large temperature ratio of the primary flow to the coolant (up to 2.5). Cooling effectiveness determined by the measured surface temperature for the two types of ceramic specimens is investigated. It is found that the multi-hole cooling system for both specimens has a high cooling efficiency and it is higher for the SiC/SiC specimen than for the oxide/oxide specimen. Effects on the cooling effectiveness of parameters including blowing ratio, Reynolds number and temperature ratio, are studied. In addition, profiles of the mean velocity and temperature above the cooling surface are measured to provide further understanding of the cooling process. Duplication of the key parameters for multi-hole cooling, for a representative combustor flow condition (without radiation effects), is achieved with parameter scaling and the results show the high efficiency of multi-hole cooling for the oblique hole, SiC/SiC specimen. (C) 2008 Elsevier Ltd. All rights reserved.

Bubble Behavior and Heat Transfer in Quasi-Steady Pool Boiling in Microgravity

Pool boiling of degassed FC-72 on a plane plate heater has been studied experimentally in microgravity. A quasi-steady heating method is adopted, in which the heating voltage is controlled to increase exponentially with time. Compared with terrestrial experiments, bubble behaviors are very different, and have direct effect on heat transfer. Small, primary bubbles attached on the surface seem to be able to suppress the activation of the cavities in the neighborhoods, resulting in a slow increase of the wall temperature with the heat flux. For the high subcooling, the coalesced bubble has a smooth surface and a small size. It is difficult to cover the whole heater surface, resulting in a special region of gradual transitional boiling in which nucleate boiling and local dry area can co-exist. No turning point corresponding to the transition from nucleate boiling to film boiling can be observed. On the contrary, the surface oscillation of the coalesced bubble at low subcooling may cause more activated nucleate sites, and then the surface temperature may keep constant or even fall down with the increasing heat flux. Furthermore, an abrupt transition to film boiling can also be observed. It is shown that heat transfer coefficient and CHF increase with the subcooling or pressure in microgravity, as observed in normal gravity.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

A load simulation method of piezoelectric actuator in FEM for smart structures

More and more piezoelectric materials and structures have been used for structure control in aviation and aerospace industry. More efficient and convenient computation method for large complex structure with piezoelectric actuation devices is required. A load simulation method of piezoelectric actuation is presented in this paper. By this method, the freedom degree of finite element simulation is significantly reduced, the difficulty in defining in-plane voltage for multi-layers piezoelectric composite is overcome and the transfer computation between material main direction and the element main direction is simplified. The concept of simulation load is comprehensible and suitable for engineers of structure strength in shape and vibration control, thereby is valuable for promoting the application of piezoelectric material and structures in practical aviation and aerospace fields.

Detection of constitutive heterodimerization of the integrin Mac-1 subunits by fluorescence resonance energy transfer in living cells

Macrophage differentiation antigen associated with complement three receptor function (Mac-1) belongs to beta(2) subfamily of integrins that mediate important cell-cell and cell-extracellular matrix interactions. Biochemical studies have indicated that Mac-1 is a constitutive heterodimer in vitro. Here, we detected the heterodimerization of Mac-1 subunits in living cells by means of two fluorescence resonance energy transfer (FRET) techniques (fluorescence microscopy and fluorescence spectroscopy) and our results demonstrated that there is constitutive heterodimerization of the Mac-1 subunits and this constitutive heterodimerization of the Mac-1 subunits is cell-type independent. Through FRET imaging, we found that heterodimers of Mac-1 mainly localized in plasma membrane, perinuclear, and Golgi area in living cells. Furthermore, through analysis of the estimated physical distances between cyan fluorescent protein (CFP) and yellow fluorescent protein (YFP) fused to Mac-1 subunits, we suggested that the conformation of Mac-1 subunits is not affected by the fusion of CFP or YFP and inferred that Mac-1 subunits take different conformation when expressed in Chinese hamster ovary (CHO) and human embryonic kidney (HEK) 293T cells, respectively. (c) 2006 Elsevier Inc. All rights reserved.

Theoretical investigation on the one- and two-photon absorption properties for a series of symmetrical charge transfer fluorene-phenylene and fluorene-thiophene derivatives

One- and two-photon absorption properties of a series of fluorene derivatives with symmetrical charge transfer D-IT-D and A-IT-A structural motifs have been theoretically investigated with ZINDO/S method. The optimized structures and the characterization of frontier molecular orbitals were obtained by using AMI calculations. Two-photon absorption properties of molecules have been studied using three-state model. The calculation results have shown that fluorene-thiophene derivatives exhibit larger two-photon absorption cross-section as compared with other studied molecules. To illustrate the results, the crucial effects of thiophene ring on fluorenethiophene derivatives and the net charge changes on the pi-conjugated bridges are analyzed theoretically. (c) 2006 Elsevier B.V. All rights reserved.

Investigation on nonradiative decay of the (Er3+I13/2)-I-4 -> I-4(15/2) transition in different tellurite glass matrix

Three kinds of Er3+-doped tellurite glasses with different hydroxyl groups are prepared by the conventional melt-quenching method. Infrared spectra are measured to estimate the exact content of OH- groups in samples. The maximum phonon energy in glasses are obtained by measuring the Raman scattering spectra. The strength parameters Omega(t) (t = 2, 4, 6) for all the samples are calculated and compared. The nonradiative decay rate of the Er3+ I-4(13/2) -> I-4(15/2) transition are calculated for the glass samples with different phonon energy and OH- group contents. Finally, the effect of OH- groups on fluorescence decay rate of Er3+ is analysed, the constant KOH-Er Of TWN, TZPL and TZL glasses are calculated to be 9.2 x 10(-19) cm(4)s(-1), 5.9 x 10(-19) cm(4)s(-1), and 3.5 x 10(-19) cm(4)s(-1), respectively.

Structural and optical properties of YAG : Ce3+ phosphors by sol-gel combustion method

High-quality Ce3+-doped Y3Al5O12 (YAG:Ce3+) phosphors were synthesized by a facile sol-gel combustion method. In this sol-gel combustion process, citric acid acts as a fuel for combustion, traps the constituent cations and reduces the diffusion length of the precursors. The XRD and FT-IR results show that YAG phase can form through sintering at 900 degrees C for 2 h. This temperature is much lower than that required to synthesize YAG phase via the solid-state reaction method. There were no intermediate phases such as YAlO3 (YAP) and Y4Al2O9 (YAM) observed in the sintering process. The average grain size of the phosphors sintered at 900-1100 degrees C is about 40 nm. With the increasing of sintering temperature, the emission intensity increases due to the improved crystalline and homogeneous distribution of Ce3+ ions. A blue shift has been observed in the Ce3+ emission spectrum of YAG:Ce3+ phosphors with increasing sintering temperatures from 900 to 1200 degrees C. It can be explained that the decrease of lattice constant affects the crystal field around Ce3+ ions. The emission intensity of 0.06Ce-doped YAG phosphors is much higher than that of the 0.04Ce and 0.02Ce ones. The red-shift at higher Ce3+ concentrations may be Ce-Ce interactions or variations in the unit cell parameters between YAG:Ce3+ and YAG. It can be concluded that the sol-gel combustion synthesis method provides a good distribution of Ce3+ activators at the molecular level in YAG matrix. (c) 2005 Elsevier B.V. All rights reserved.

Cooperative and charge transfer luminescence in Yb3+ -doped yttrium aluminum perovskite (YAlO3)

YAlO3 (YAP) crystals with different Yb3+ concentration have been grown by Czochralski method and cooperative fluorescence of Yb3+ ions in YAP crystal was studied under 940-nm infrared (IR) LD excitation at room temperature. The Yb concentration dependence of absorption intensity of IR and charge transfer bands exhibit different features. The green emission band in the region of 480-520nm was assigned to the cooperative deexcitation of two Yb3+ ions. The remaining upconverted emission bands containing various sharp peaks associated with impurity ions were observed and discussed. Charge transfer luminescence of heavily doped 20at% Yb:YAP is strongly temperature dependent and no concentration quenching of the charge transfer luminescence was found through the investigation of different Yb levels samples. (c) 2006 Elsevier B.V. All rights reserved.

Spectral properties and charge transfer luminescence of Yb3+ : Gd3Ga5O12 (Yb : GGG) crystal

The Yb (10%):GGG and Yb (30%): GGG crystals have been grown by the Czochralski method. The chemical compositions are: Yb1.07Gd1.74Ga5.19O12 and Yb0.33Gd1.47Ga5.2O12. The absorption and emission spectra of Yb:GGG crystal at room temperature have been measured. The spectroscopic parameters of Yb:GGG and Yb:YAG have been compared. Optical absorption spectra of Yb:GGG show 4f-4f transitions related to Gd3+ ion around 300 nm, and also an onset of charge transfer (CT) transitions from oxygen ligands to Gd3+ or Yb3+ cations below 240nm. The CT absorption of Yb3+ is largely overlapped by that of Gd3+ ions. (c) 2005 Elsevier B.V. All rights reserved.