融液热对流流动不稳定性分析研究.pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号；061076 融液热对流流动不稳定性分析 刘旭良敬成君 两南交通大学机械工程学院．成都，邮政编码610031 Tcl．028 com．cn E·mail：chengiunjing@yahoo 摘要；奉文班Cz法生氏氧化物晶体系统中融液热对流为对象，在仅考虑浮力和表面张力的情况下， 通过直接数值计算(DNS)求解出不同大小组合的浮力和表面张力作用下的热对流场。通过对DNS结果 的分析。定性得到表面张力和浮力在流场不稳定性中各自起主受作用的参数区域。在此区域内再采用 正交分解法，抽取出流场的最主要基本流场模式，通过对其主要模式的分析比较，得到表面张力和浮 力对流场不稳定性备自起主要作用的情况和前面DNS的分析结果足一致的。结果表明t正交分解法可 用以分析流场的不稳定性。 关键词：正交分解，最基本流场模式，cz晶体生长，数值模拟 1、引言 流动稳定性问题是一个广泛受到关注的问题。比如，在Czochralski单晶氧化物制 备系统中，稳定的融体对流是晶体制备的关键，融体的稳定性直接影响生产工艺控制和 晶体品质。但是由于浮力、表面张力、晶体旋转引起 的离心力等影响因素，使得Cz晶体生跃系统融液内的 热对流现象错综复杂，热对流不稳定的原因难以确定。 因此，要有效抑制热对流，获取高品质晶体，有必要 弄清流动形成时各种不稳定性对热对流的影响。分析 稳定性的方法虽然很多，但它们都不能同答在多种不 稳定性因素存在的情况下，每个不稳定性囡素对形成 对流的贡献是多少。本文的目的在于探索回答上述问 题的新方法一正交分解法一的可行性。 2、物理数学模型 图1简化模型 根据研究目的，本文以一个简单的物理模型展开研究。图1所示轴对称cz单品体 生长系统中融体对流的驱动力只包含浮力和表面张力。假设融液为不可压缩．满足 Boussinessiq假设的牛顿流体；流动为层流；粘性耗敞忽略。丁是，描写融体热对流的 控制方程可表达如卜．： 生+生+丝。0 f1) R引i8Z 鲁+一等+旷等。一篆+c一号+去等+豢+·等，㈠ 毕套项H：尉家自然科学缺会资助项日(、。年¨教育部留学川罔堆金资助 42 c。， 鲁+一案+屹等■善+唼等+豢+参，+ar(T—D 瓦aT+‘器+圪詈一百1(_．1丽aT+紊+争aT ’aR 。az Pr、RaRaR2az。 ct， 由于简化模型中朱考虑坩埚和晶体的旋转，敌周方向速度方程不出现。上述控制方程是 以，c、v／r,、pv 2／c、f／’，和乙为特征长度、速度，压力、时间和温度无量纲 化后得到的，其中P为融液密度，L为晶体融点，v为融液运动粘性系数。控制方程 中出现的无量纲特征数的定义为格拉晓夫数GrI蝉，普兰特数丹-兰，其中4为 V‘ 口 热扩散系数，g为重力加速度，口为热膨胀系数。 本研究边界条件如下：底部和中心轴上绝热；周边边界定温疋；结晶区定温L。 自由表面向空间辐射换热，满足一罟．‰口·一F)，其中R。。!号乏。A为融 ^ a厶 液导热系数，s为融液辐射发射系数，0为斯蒂藩～玻耳兹曼常数，L为坩埚外