柴式長晶法生長氧化鋁單晶過程中，為了提升氧化鋁單晶的生長品質，必須了解柴式長晶爐內部的熱傳與熔湯流動情形。由於單晶生長過程無法直接觀察熔湯內的熱流場分佈，所以本研究使用有限體積法之CGSim軟體模擬二維軸對稱氧化鋁單晶生長過程之電磁場、熱場與流場。求解的方式先從RF線圈產生電磁場，感應坩堝計算出熱通量分佈之後，在以此為邊界條件來修正並同時解出整個爐體的熱場與流場，以此溫度場再對晶體、熔湯、坩堝、以及長晶區域上爐室內之氣體來做更進一步的深入計算，最後將晶體的溫度分佈數據匯入至COMSOL做三維非等向熱應力分析。 本研究主要模擬柴式法生長氧化鋁單晶的生長過程，探討生長過程中晶體內部熱輻射對於熱流場、固液界面形狀與熱應力分佈之影響。為了做出更理想的長晶模擬，本研究加入了氬氣的流動、熔湯表面張力、考慮熔湯為紊流模式以及晶體內部熱輻射的影響，並針對表面張力、晶體內部輻射與等效熱傳導係數對於熱流場分佈之影響作分析。最後計算出熱應力的分佈用來分析生長c軸與a軸之氧化鋁單晶的缺陷分佈情形。 The thermal and flow transport in a Czochralski crystal growth furnace plays an important role to effect the single crystal growth quality of sapphire. However, the thermal and flow fields in the melt of the single crystal growth process are difficult to observe in experimental study. This thesis has numerically investigated the thermal and flow transport phenomenons using the finite volume method via CGSim software. The electromagnetic, thermal, and fluid fields during the sapphire single crystal growth process have been investigated. The temperature and flow fields inside the furnace are coupled with the heat generation in the Iridium crucible which is generated by the electromagnetic field using the RF coil. For different stages of Czochralski crystal growth of sapphire, the configuration usually used in a real Czochralski crystal growth process consists of a crucible, induction coil, insulation, crystal, melt and gas. At first the volumetric distribution of heat generation inside the metal crucible is calculated. Using this heat distribution as a source the fluid flow and temperature fields in the whole system are computed. After that, thermal analysis is computed only in the crystallization zone, including the crystal, melt, crucible, and a gas region around the crystal. The computation in the crystallization zone involves the turbulent flow of the sapphire melt, the laminar gas flow, and radiative heat exchange in the semitransparent crystal. Then using the temperature distribution in a single crystal to calculate thermal stress via COMSOL software. The anisotropy of elastic constants and thermal expansion coefficients are considered in the problem. The analyses are performed both for the pulling directions of a-axis and c-axis.
