环形浅液池内热毛细对流热力学特性.pdfVIP

中国一程热物理学会 工程热力学与能源利用 2005、年学术会议论文 编号：051110 环形浅液池内热毛细对流的热力学特性’ 李友荣唐经文黄旭方 曾丹苓 重庆太学动力工程学院i熏庆400044 1 TEL：023-65 12284，E—mail：liyourong@yahoo．corn 摘要：为丁了解微重力下水平温度梯度作用时环形浅液池内硅熔体热毛剀对流的热力学特性，利用有 限差分法进行7非稳态三维数值模拟。结果表Ⅲj，当温度梯度较小时，流动为轴对憨稳定流动，系统 总熵产较小；随着温度梯度的增加≯流动将失去儿稳定性，首先转化为径向脉动波，系统总熵产呈周 期性变化；再增加温度梯度时。流动转化为热流体波，系统总熵产较大，但不再随时间变化。 关键词：环形浅液池；热毛绷对流：数值模拟；熵产 1．引言 当自由表面存在温度不均匀时，会产生表面张力的不均匀，从而驱动流体运动，这 种由丁二表面张力驱动的流动称为热毛细对流。热毛细对流广泛存在于薄膜喷涂…、液滴 蒸发“1、晶体生妖”1等工程技术领域，并对实际生产过程有着重要的影响，网此，数十 年来已进行了较多研究。Smith和Davis”’完成了水平温度梯度作用下县有自由界面的无 限大水平薄液层内热毛细对流的线性稳定性分析，结果发现了两种不稳定性，～种是稳 了热流体波的存在“…。另一方面，在}扭壁加热的圆柱形容器里，Kamotani等…在微重力 下热毛细对流试验中观察到了2．3片叶状图案；Schwabe等”“’在地面和人造卫星上进行 的系列环形浅液池内热毛细对流试验中，观察到了螺旋状的短波和长波两种热流体波模 式，并发现重力对热毛细对流有稳定作J；：|j；Gamier等““1观察到了螺旋臂形热流体波纹 和靠近冷却内壁的脉动靶形波纹；同时，Sim等…I”和LI等…“2。通过数值模拟也证实了环 形液池内热流体波和I驻波的存在。然而，到目前为止，对液池内热毛细对流的热力学特 性的研究非常少，本文拟通过三维数值模拟探讨流烈变化与系统熵产的关系，为揭示热 毛细对流的机理提供一条新的途径。 2．物理数学模型 5 物理模型如图1所示，内径为一=1mm，外径为r。=50mm，深度为d=3mm，上部 底部平lj项部均绝热。在自由表面考虑热毛细力的作刚，在其它液一吲界面满足无滑移条 性都为常数；2)流体速度较小，流动为层流；3)自由表面为平面且不变形。 在上述假定条件F，连续性方群、动量方稃和传热方程分别为： ’ 目家自然科学捧每资助项日(N 123 Vtv=0 (1) ’ 宴+，．珊：一上十胛zv， (2) 就p ar+v．Vr：dV2T (3) Ot 图1物理模型 其中，v为流体速度矢量，f为时间，础硅熔体密度，P为压力，劝动量扩散率，日为 热扩散率。热毛细对流过程熵产率为； ‘ 驴砉I蚓2+手伊 (4) 式中嘞粘性耗散，第一项为不可逆传热引起的熵声， 第二项为粘性流动引起的熵产。 系统的总熵产为 sz=』jgdV． 。 (5) 边界条件为： s 自由表面0司rirro，0OZn)， ∥誓=埽iOT，乒虿0vo=YT焉∥吉。埽i一乒虿2丽 偿 底部臼=o，rj，‰