分离结晶过程中熔体热毛细对流的数值模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：073403 分离结晶过程中熔体热毛细对流的数值模拟木 彭岚范菊艳李友荣 (晕庆大学动力T程学院，重庆400044) 1 Tel：023—6511283，E—mail：pen91an@cqu．edu．cn 摘要：为j，J，解新型分离结晶生K中熔体热毛细对流的蕞本特征，利用有限差分法进行’r数值模拟，熔体深径比 某一临界值后，流动将转化为1r稳态流动。 关键词：热毛细对流，分离结晶，数值模拟，临界条件 I引言 耸得到极火提高，但关丁．此方法的研究目前只是在微重力条件‘卜-取得了实验上的成功【l叫。在地 面条件下，Sz0纳m及其合作者【5，6j川安瓿进行直径1厘米左右的晶体生长实验时，就没能完全 实现坩埚不浸润(cruc．bJede．w酣in2)现象，实验过样中晶体重新粘附在安瓿壁上的现象时有 发生，其根本原冈就在丁-缺乏系统的理论研究，没有真正掌握分离结品时熔体内的流动和传输 规律。本文报道了一组分离结晶Bndgman生长过程中熔体热毛细对流的数值模拟结果，重点考 察了流型的演变过程及流动的失稳机制。 2物理数学模型 物理模型如图l所示，一晶体?卜径为，f，坩埚半径为％，熔体高度为厶，顶部为州蹙，结晶界 面和坩埚外接分别维挣噎定温度％和死(乃，功。由于沿自由表面存在温度梯度，冈此，在表 面张力梯度的驱动F，在熔体内会诱发热毛细对流。为简化起见，假定：1)熔体为不可压缩的 牛顿流体；2)流速较低，流动为二维层流：3)在自由界面考虑热毛细力的作刚，．在其它闻譬 满足无滑移条件；4)表面张力是温度的线性函数。 分别将％7¨％，吮用忉吮，2(沩流体动量扩散率，础密度)作为控制方程的无量纲参 考时间、长度、速度和压力，则无鼙纲化后的控制方程为 V．y=0， (1) 奠■，r !二+y。Vy：一V尸+V2y， (2) af 里竺+y．vp：上V2p． (3) ar 厅 ． 式中，y为无冈次速度矢鼙，沩无冈次时间，尸为无冈次压力，6目zZ，∥力瓦一乃≯为无冈次温度， 办=妇为熔体普朗特数，口为热扩散率，。 +恭金项日：周家自然科学耩会资助项目(No．50676】12) 1424 边界条件如下： 坩埚壁，胆l： 驴卿，伊1． (4∽) 中心处，R卸： (5a-c) 忙0，at／强=0，ap／棘=0． 顶部固壁处，蒯： ‘瞄陉O，ap／彪=0． (6a．c) 、结晶界面：z-旬，0S飚R： 睁垮=o，p=0． (7a．c) 图1物理模型 自由表面，Z_0，鄹1： 痞o，等：兰票，筹=0． 、(8a-c)。 ‘；8zB．PraR；aZ 其中，R和Z分别为无因次坐标，萨“∞忆为无因次自由表面宽度，^勿=乃(％乙№吨岫回 为Ma瑚goni(呦数，为；也怕r为表面张力温度系数，乒沩粘性系数，n沩热扩散率。 定义无因次流函数油 y：一』旦．丝矿：土旦竺 r