Comparison of thermal stress computations in Czochralski and Kyropoulos growth of sapphire crystals

• Published in: Journal of crystal growth 2018-10, Vol.499, p.77-84
• Main Authors: Stelian, Carmen, Sen, Gourav, Duffar, Thierry
• Format: Article
• Language: English
• Subjects: A1. Computer simulation; A1. Stresses; A2. Kyropoulos method; Asymmetry; B1. Sapphire; Comparative analysis; Computer simulation; Convection furnaces; Crystal growth; Crystals; Czochralski method; Dislocation density; Elastic anisotropy; Elastic properties; Engineering Sciences; Ingots; Materials; Mathematical models; Sapphire; Simulation; Stress analysis; Stress concentration; Stress distribution; Stress state; Thermal stress; Three dimensional analysis; Two dimensional models
• ISSN: 0022-0248; 1873-5002
• DOI: 10.1016/j.jcrysgro.2018.08.002

•3D numerical computations of stress in sapphire crystals grown by Cz and Ky methods are performed.•Kyropoulos grown ingots exhibit high thermal stress only in a thin region at the crystal periphery.•Simulations of Czochralski technique show high stress in almost whole volume of the crystal.•Numerical results are in agreement with experimental observations about dislocation densities.•Computations with temperature dependent elastic constants provide better results. Thermal stress computations during sapphire growth are compared between a resistive Czochralski furnace and a Kyropoulos inductive furnace. 2D – axisymmetric global simulations are performed to compute the thermal field in the furnace and the convection in the melt. Temperatures carried out from global modeling at a given stage of the growth process, are used for thermal stress computations in the crystal. Three-dimensional stress analysis, which takes into account the anisotropic elastic constants of sapphire, shows nearly axisymmetric von Mises stress distribution in the crystal. It is shown that applying 2D – axisymmetric modeling of thermal stress may result in significant errors. Computations performed for a crystal of 10 cm in diameter grown in a Kyropoulos furnace show only a thin region of 2–3 mm with high thermal stress located at the crystal periphery. The model predicts very low thermal stresses in the central part of the crystal. Simulations of an ingot grown by Czochralski technique, show higher thermal stresses in almost the whole volume of the crystal. Numerical computations are in agreement with our previous measurements of dislocation density in sapphire crystals grown by the Kyropoulos method. The present numerical results can explain the experimental observations showing that sapphire crystals grown by Czochralski technique have a much higher dislocation density than Kyropoulos grown ingots.