单晶生长炉全局热分析3——自由界面形状的影响.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：073239 单晶生长炉全局热分析木(3) …自由界面形状的影响 贾琛霞，敬成君，高辉 西南交通大学机械工程学院，邮编：61003l 1’el：028 Email：chengjuniing(国126．com 摘要：本文应用已建立的考虑了熔液对流三维性和非定常性的全局热分析模型，通过大量数值模拟， 研究了弯曲自由表面对熔液流动、炉内温度场分布等的影响。研究表明，弯曲自由表面更易引起熔液 流动振荡，降低晶体成长界面发生反转的临界Re数。 关键词：全局热分析，熔体对流，界面反转，弯曲自由表匾 1引言cz熔液法生长单晶体是主要的单晶生长方式。单结晶的品质与结晶成长时炉 内的流动、传热传质等输运现象密切相关。为了生产大尺寸(12英寸以上)7高品质单晶 体，正确把握生长炉内的所有现象是非常必要的。将炉内的所有构成要素全部考虑的流 动、传热传质分析称为全局热分析。然而，由于炉内高温热辐射和熔液非定常三维热对 流等传输现象相互耦合、错综复杂，迄今为止，全局热分析模型一直停留在传统的轴对 称准静态模型的假定上[1]。本文作者发展了传统的模型，在前文报告了考虑熔液流动 三维性和非定常性的全局热分析模型并研究了非定常三维熔液流动的影响。但是，在前 文模型中，熔液自由表面被假定为平坦的表面，这样虽可大大简化自由表面边界条件的 处理，但与实际情况不符，这可能也是计算得到的固一液界面反转临界Re数仍然比实际 值大的原因。根据实际，熔液自由表面与晶体的接触角大约为17。，与坩埚的接触角大 约为90。，因此，溶液自由表面在晶体边缘附近相当陡，在坩埚附近才较为平坦。本文 考虑自由表面弯曲情况，建立更完善的全局热分析模型，并研究弯曲自由表面对炉内温 度场和熔液流动 的影响。 2问题描述 本文研究对 象仍是来自生长 LiNbOs晶体的实 验装置(图1)， 不过，在模型中 考虑熔液自由表 面的弯曲情况。 这时，熔液高度 为50·8mm(平坦 图1 cz生长炉艇况 固2坐标系及格子 +国家自然科学基金资助项目 1288 自由表面假定时为53．3唧)。其 余条件同前文，即重要的参数为 竹=13．6，函=4．67x 105。这里 ·d娟n蘸bn pOInt 主要报告关于自由表面边界条件 的处理技巧和全局热分析的结 果。 _撇t髓e 2．1边界条件的处理 自由 暑u—訇lc霉 表面形状由Young—Laplace方程 芎 求得，局部示意如图3。、根据自 妇I疗eree撵 由表面剪切应力为O，则速度边 界条件为： 0．仃．，2=O，f，．仃．，z=o，矿·n=O， 0 圈3目由表面边界格子 这里，f和撑是自由表面的切向和法向向量，刀是应力两张，y是速度向量。下标孝和 77表示曲线坐标。数值模拟时，一个时间步长内，在假想外部点如(fJZ”约的值给定的情 况下，通过求解控制方程，内部点如(‘工刀舡J)的值被求出，然后根据边界条件更新假想 外部点的值。上述边界条件方