浮力对环形浅液池内热毛细对流影响.pdfVIP

中固T程热物理学会 传热传质学 2005年学术会泌论文 编蛩：05316 浮力对环形浅液池内热毛细对流的影响+ 石万元a李友荣“··曾丹苓4今石宣之6 ‘重庆大学动力．T-程学院．重庆400044．6日本九州大学先导物质化学珂f究所 电话：023E-maih liyourong@yahoo．COm 摘要：通过对外壁受热、内壁冷却，厚度为I、2mm的环形浅液池内硅油热毛细对流的_三维数值模拟． 从表面热流体波和内外壁面的对流换热两方面分析，浮力的影响。一般地，祚浮力作用下，随液池厚 度的增加．浮力的影响增加，外壁处波动振幅高于内壁，热流体波的波数有所增加。但总体上热流体 波的振幅下降，内外壁面局部Nu数的涨落降低，表明重力对热流体激具有稳定作用。 关键词：数值模拟．浮力，热毛纠对流．硅油，热流体波 1．引言 随着对半导体单品生长质量要求越来越高，以及微米、纳米尺度流动传热现象等研 究的深入，热毛细对流越来越引起学者们的注意。尤其是当界面的温差超过临界值时， 单晶生长过程中，热毛细对流对晶体质量有不利的影响，多年来相关研究报道较多【14】。 出环形液池内沿周向传播的热流体波和驻波口】。Li等模拟了硅熔体和硅油的热毛细对 流，发现在较浅的液池里热流体波早弯曲的轮辐状，而在较厚的熔体里为真线轮辐状 ”“…。在低于2mm厚的液池里，重力对热毛细对流的影响的数值模拟鲜有报道，为此， 这里通过数值模拟讨论l、2mm厚硅油环形池内重力作用下的热毛细对流。 2．模型及数值方法 mm、深度为d=-I、 物理模型如图1所示，环形液池|_i：|径为Yi=20mm、外径为ro=40 绝热，内、外壁分别维持恒温咒雨I 顶部自由表面考虑热毛细力的作用，而在所有同鼙满足无滑移条件。 餐纲参考时间、|丈度、速度雨I压力，则无量纲化后的控制方程为： +同家白然科学举会资助项M(No ¨联系人 VV=0， (1) 娑+Vw：研+v：y+G舰：，(2) dr 一 望+y．v。：土甲：。 、。(3) ar Pr 边界条件如F： 自由界面(Z=d／ro，O．5(R1，0≤出2州 图1物理模型 里：Re塑，篓：Re旦，矿：o，盟：o． H az aRlaz R00’一 aZ 底部(z=oJ0．5尺1，0≤02叫 ¨：Vo：眨：o．鲁：o． 6 内壁(R=O5，0≤Z≤d仫0≤02叫 ■=go=K=0．0=0，=0 嫡 外壁(尺=1，0≤z≤d／ro,0s02叫 一=％=圪=0，0=0。=I． 0 其中，V为无冈次速度矢量．P为无冈次压力，日=rF疋』吖功一TO为无困次温度，沩无 表面张力温度系数，口为热扩散率，助体积膨胀系数，g为重力水平，勘为重力加速度。 200C时硅油的物性参数如表1所示。 袤1 200C时硅油的物性参数 计算区域被划分为82‘X(16．21)2X1838个非均匀(周向除外)的交错网格。控制方程 X10。秒至2X