融液自由表面波状流动影响因素的研讨.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：053388 融液自由表面波状流动影响因素的研究 林松。敬成君矿蟓田隆夫5 “西南交通大学机械学院，邮编610031．“I==|本东北人学多冗研 (wave 式，通过火规模数值模拟，分析了晶体旋转速度、自由表面热辐射、表面张力搜坩埚壁加热条件对波状 流动的影响。 关键词：Czochralski晶体生长，波状流动，热毛细不稳定性，数值模拟，波数 1．引言 波状流动(wave 备系统中，是一种晕现在坩埚融液自由表面上的耗散结构现象。其发现虽已有20多年， 但由于Cz生长系统中三维时相关热流动现象错综复杂，加之高温环境难以直接实验观 察，故对波状流动的深入研究存在困难，迄今为止主耍停留在冷模型实验和数值模拟。 Jones…曾用水或水与甘油的混合物分别进行过冷模型实验，在室温l=‘研究了波样流动： Leeeta1．”3研究了波状流动与Prandtl数之间的关系；最近，EngeretalI”虽第一 个成功地崩数值方法再现了波样流动，但其物理模型简单，忽略了实际高温情况F自由 表面的热辐射，且将坩埚壁理想化为定温加热条件，其结果仍不能描述实际系统。 本文以实际LiNbOs制备系 统(大尺寸坩埚)为例，通过数 值方法研究了务种因素对波样 流动形成的影响。 2．数值方法 本研究所用物理模型如图1 (a)所示。坩埚静IF，其内径 rc=100mm，高度h=r。。晶体半径 ‘8) (b) 为is,旋转速度为Q；。假定融液 圈1坐标系(a)及分块刚格(b) 1882 假定自由表面的散热以辐射方式散失到温度为T。的环境，坩埚壁的加热条件为定温度 T。或定热流密度q。；在自由表面上，分两种情况进行研究：考虑表面张力和不考虑表面 张力。 ， 描述融液热对流的数学方程如下： f筹+警卜÷嘉f警1_＆ ㈤ I打。馘J如飒L瓯J… 方程中通用变量电无因次参数R旃源项岛，以及坐标变换后的源项s：如F表 胎：” 由于本研究对象尺寸大，所需计算网格多，故采用分块技术和边界适合坐标法”l 进行大规模数值模拟。计算域网格的水平投影如图l(b)所示，总的控制容积数为476，160 个。本研究将方程(1)坐标变换后，利用全守恒控制容积集中格子法离散，空间离散精 速度场；为保证压力和速度正确耦合，，应用了Rhie“的动量插值技术：根据LiNbOs的热 104：所有计算的初始条件均为 物性计算出无因次数为Gr=2．02x107，Pr=13．7，Re／=9．9x Ut=P=T=O。 计算在日本东北大学超级计算机SRS000上完成的。 3．结果 3．1波状流动的数值模拟 算结果显示： 当Q。10rpm，浮力支配流动，流体自壁面上升后在自由表面层向内，然后在中心附 近向下流动完成循环。三维流动复杂但自由表面不出现任何流动图案； 当Q。lOrpm，由旋转引起的向外流动增强，浮力引起的向内流动在自由表面的地方 开始向下，两股流动交汇，形成波状流动。图2所示为以自由表面等温线分布表示的 n。=15rpm时波状流动的演变过程。从均匀的温度场出发，最先形成规则的波数为m=4 的流动图样，后来消失变成不规则m=3的流动图样。随着时间的发展，表面图样在m=3 和m=4间来回变化。经长时间发展，最终表面波数稳定在波数m=3上。 当n。n。。(临界转速)时，旋转引起的向外流支配整个自由表面，波状流动消失。 数值方法结果与实验观察定性一致。 · (a)30s (b)60s (c)120s(d)160s(e)250s(f)380s