数值模拟环形液层浮力对流与热毛细对流的相互作用.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093299 数值模拟环形液层浮力对流与热毛细对流 的相互作用1 石万元，李友荣，彭岚 (重庆大学动力工程学院，重庆400044) Tel ：023； Email ： shiwy@cqu.edu.cn 摘要：为了理解浮力对热毛细对流的影响，二维数值模拟了 10mm 深环形硅熔体(Pr=0.011)和硅油 (Pr=6.7)液层内浮力对流、热毛细对流及浮力－热毛细对流，结果表明：对硅熔体，浮力对流较强， 热毛细对流较弱，浮力强化了热毛细对流，使热毛细对流更容易失稳；对硅油，热毛细对流较强，浮 力对流较弱，浮力削弱了热毛细对流，使热毛细对流趋于稳定。数值模拟结论澄清了浮力对流与热毛 细对流的作用机理，揭示了浮力对不同普朗流体热毛细对流及其稳定性影响的根源。 关键词: 热毛细对流, 浮力, 硅熔体, 硅油 0. 引言 在地面常重力环境下，直拉法单晶体材料生长过程中浮力和热毛细力都是导致熔体对 流的重要驱动力，它们对所生长的单晶体材料质量具有不可忽视的影响[1] ，数十年来， 国内外不少学者对浮力对流或热毛细对流以及它们对单晶体材料制备过程的影响进行 了广泛深入的研究，有的讨论浮力对流，有的讨论热毛细对流，近来，一些学者讨论 了浮力－热毛细对流。例如，Schwabe 等[2-4]在微重力和常重力环境下，实验研究了环 形液池内硅油 (0.65 cSt) 的热毛细对流和浮力－热毛细对流，发现浮力对热毛细对流具 [5] [6] 有稳定作用。Peng 等 通过数值模拟验证了Schwabe 等人的结论。我们 对 1mm 深的 浅液池内热毛细对流数值模拟也发现浮力使热毛细对流趋于稳定。 然而，对于其中的机理，目前还缺乏清楚的解释。譬如，浮力本身是竖直向上 的，高温壁附近的热毛细对流也是竖直向上的，二者方向相同，浮力似乎应该使流动 趋于更强烈才对，然而，为何实验和数值模拟均发现浮力对热毛细对流具有稳定作用 呢？为了回答这个问题，本文数值模拟了具有不同普朗特数的硅熔体(Pr=0.011)和硅油 (Pr=6.7)在环形液层内的浮力对流、热毛细对流、以及浮力－热毛细对流，分析了浮力 对热毛细对流及其稳定性影响的根源。 1. 物理数学模型及计算方法 [2-4] 我们的物理模型与 Schwabe 等 的实验模型一样：一内径 Ri 为 20mm 、外径 Ro 为 40mm 、深度为 10mm 的环形池外壁维持高温 T ，内壁低温 T ，上表面为自由面，底 h c 面为固壁，假设底面和自由面绝热。流体为硅熔体和硅油，它们的物性参数见表 1。假 设流体为不可压缩流体，满足 Boussinesq 假设，对热毛细对流和浮力－热毛细对流， 考虑自由表面的热毛细力。由于在较小的内外壁温差下热毛细对流是轴对称的，而且 我们的主要研究目的侧重于浮力对热毛细对流的影响，稳定的轴对称流动更便于分 析，因此，本文的数值模拟以二维轴对称对流为研究对象，几何模型如图 1 所示，它 表示环形池沿径向的剖面图。 数学模型的无量纲形式表示为： 1 基金项目：重庆市自然科学基金项目(CSTC2008BB4301), 中国博士后基金项目(20080440694), 重庆大 学青年骨干教师创新能力培育基金项目(CDCX017) 图1 几何模型 表1 硅熔体和硅油的物性参数 硅熔体 硅油