采用C语言搭建了碳化硅(SiC)晶体生长炉三维温度场数值模拟平台，基于柱坐标系构建生长炉物理模型，采用有限体积法离散数学模型，利用S2S(Surface to Surface)辐射模型考察了生长室内的辐射换热特性，提出判断辐射面可视性的最短距离法. 模拟了电流强度1250 A、电流频率16 kHz条件下的生长炉温度场，定量揭示了生长室内的辐射换热强度；采用标准偏差法研究了线圈结构对晶体内部温度及温度梯度均匀性的影响. 结果表明，螺旋电磁加热线圈容易导致生长炉内部温度场呈非轴对称特性分布；辐射热流较导热热流大102~103倍，辐射换热促使生长室内温度分布均匀；螺旋线圈的布置方式使晶体截面温度呈非轴对称分布，造成温度梯度的均匀性变差，晶体生长过程中容易产生热应力，影响晶体质量.
Three dimensional numerical simulation platform for SiC crystal growth furnace was established based on C programing language, where physical model of the furnace was built based on cylindrical coordinate, governing equations for electromagnetic and temperature fields were discretized by finite volume method, and radiation characteristic was studied with the help of S2S (Surface to Surface) radiation model. The least distance method was developed, which was used to check whether radiation surfaces were visible with each other or not efficiently. And then the radiation heat transfer in SiC growth chamber and temperature field of SiC growth furnace were studied quantificationally when the current intensity is 1250 A and the current frequency is 16 kHz. The effects of coil structures on crystal temperature field and its gradient distributions were studied by standard deviation method. The results showed that spiral electromagnetic coil generated asymmetrical temperature field easily. The radiation heat flux was 102~103 more than conduction heat flux. The radiation was helpful to increase temperature evenness. The spiral temperature field on the SiC crystal cross-section reduced the poor homogeneity of temperature gradient, which made the crystal to generate large thermal stress.
[1] Shiramomo T, Gao B, Mercier F, et al. Thermodynamical Analysis of Polytype Stability During PVT Growth of SiC Using 2D Nucleation Theory[J]. Journal of Crystal Growth, 2012, 352(1): 177-180.
[2] Hofmann D, Winnacker A, Durst F, et al. On the Sublimation Growth of SiC Bulk Crystals: Development of a Numerical Process Model[J]. Journal of Crystal Growth, 1995, 146(1/4): 214-219.
[3] B?ttcher K, Schulz D. Computational Study on the SiC Sublimation Growth[J]. Journal of Crystal Growth, 2002, 237–239(Part 2): 1196-1201.
[4] Shi Yonggui,Dai Peiyun,Yang Jianfeng, et al. Thermal Stress Analysis of Crystal Morphology in Growth of Silicon Carbide Bulk Crystal by PVT Method[J]. Joural of Xi’An Jiaotong University, 2011, 45(1): 117-121 (in Chinese).
史永贵, 戴培赟, 杨建锋, 等. 物理气相传输(PVT)法生长碳化硅的晶体形态热应力分析[J]. 西安交通大学学报, 2011, 45(1): 117-121.
[5] Meng Daleng. Thermal Field Simulation of the Growth Process of PVT Silicon Carbide[J]. Semiconductor Optoelectronics, , 2014, 35(3): 461-463+501 (in Chinese).
孟大磊. PVT法碳化硅单晶生长系统的热场模拟[J]. 半导体光电, 2014, 35(3): 461-463+501.
[6] Zhang Qunshe, Chen Zhiming, Pu Hongbin, et al. A Hydrodynamics Model of Silicon Carbide Single Crystal Growth by PVT and Finite Element Analysis[J]. Journal of Synthetic Crystals, 2005, 34(5): 69-73 (in Chinese).
张群社, 陈治明, 蒲红斌, 等. SiC晶体PVT生长系统的流体力学模型及其有限元分析[J]. 人工晶体学报, 2005, 34(5): 69-73.
[7] Chen Q S, Zhang H, Prasad V. Heat Transfer and Kinetics of Bulk Growth of Silicon Carbide[J]. Journal of Crystal Growth, 2001, 230(1–2): 239-246.
[8] Klein O, Philip P. Transient Numerical Investigation of Induction Heating During Sublimation Growth of Silicon Carbide Single Crystals[J]. Journal of Crystal Growth, 2003, 247(1–2): 219-235.
[9] Chen Q S, Zhang H, Ma R H, et al. Modeling of Transport Processes and Kinetics of Silicon Carbide Bulk Growth[J]. Journal of Crystal Growth, 2001, 225(2–4): 299-306.
[10] Xue S, Proulx P, Boulos M I. Extended-field Electromagnetic Model for Inductively Coupled Plasma[J]. Journal of Physics D: Applied Physics, 2001, 34(12): 1897-1906.
[11] Zhu Hai-long,Tong Hong-hui,YE Gaoying, et al. Numerical Simulation of Two-dimensional Inductively Coupled Plasma Torch Based on FLUENT[J]. Nuclear Fusion and Plasma Physics, 2012, 32(3): 199-205 (in Chinese).
朱海龙, 童洪辉, 叶高英, 等. 基于FLUENT的二维感应耦合等离子体炬数值模拟[J]. 核聚变与等离子体物理, 2012, 32(3): 199-205.
[12] Tao Wenquan. Numerical Heat Transfer[M], The Second Edition. Xi’An. Xi’An Jiaotong University Press, 2001: 92-99 (in Chinese).
陶文铨. 数值传热学[M], 第2版. 西安: 西安交通大学出版社, 2001: 92-99.
[13] Yang Mo;Wang Yuqing;Fu Yanhong; et al. A Numerical Method for Coupled Conduction and Convection Heat Transfer with Surface Radiation [J]. Joural of Xi’An Jiaotong University, 1992, 26(6): 25-31 (in Chinese).
杨沫, 王育清, 傅燕弘, 等. 具有表面辐射的导热和对流耦合问题的数值计算方法[J]. 西安交通大学学报, 1992, 26(6): 25-31.
[14] B?ttcher K, Cliffe K A. Three-dimensional Thermal Stresses in on-axis Grown SiC Crystals[J]. Journal of Crystal Growth, 2005, 284(3–4): 425-433
