水平温差下环形浅液池内热毛细对流的渐近解 刘佳.docVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：143538 水平温差下环形浅液池内 热毛细对流的渐近解 刘佳 彭岚 张全壮 （低品位能源利用技术及系统教育部重点实验室，重庆大学动力工程学院，重庆 400044 Tel: 023 E-mail: penglan@cqu.edu.cn 摘要：为了研究水平温差下环形浅液池内热毛细对流的特征，采用渐近线方法求得了环形浅液池内热毛细对流的近似解析解，得到主流区域的速度场、温度场以及自由表面形状。环形池保持内壁低温、外壁高温、自由表面与空气对流换热。 关键词：热毛细对流； 0 前言 热毛细对流是当存在界面时，由不均匀的界面张力拖曳流体引起的。由于热毛细对流在工程应用中广泛存在，对它的研究一直方兴未已。数值模拟虽然已成为研究热毛细对流的主要方法之一，但由于分析解便于探讨物理参数对流动的影响，在此方面也引起了广泛的研究[1-4]。 Cormack，Leal和Imberger[5]用渐近法分析了边壁具有不同温度的矩形浅液池内流体的自然对流。Sen和Davis[6]对矩形池内热毛细对流展开研究，分别得到了不同界面条件下中心区域和边壁区域的近似解。Prakash和Koster[7]采用渐近法、数值模拟及实验方法，对双层流体内热毛细对流进行了分析，结果表明中心区域的渐近解与数值模拟结果具有高度的一致性。 目前在考虑表面变形的情况下，针对热毛细对流的渐近解的研究多数集中于矩形腔[8-10]，本文拟将其拓展至环形池的研究，以期获得环形浅液池内热毛细对流的速度场、温度场以及界面形状。 1物理数学模型 物理模型如图1所示：环形液池内径为ri，外径为ro=ri+l，l为液池宽度，深度为d，内、外壁分别维持恒温Tc和Th（ThTc）=h(r)表示，且维持不变，底部为绝热固壁。假设流体为不可压缩牛顿流体，除表面张力外，满足Boussinessq近似，表面张力是温度的线性函数。定义纵横比为：；定义为：；定义为。 假定水平温差较小，流动为轴对称二维流动，分别将d，和(为流体动力黏度，为表面张力温度系数，为温差)作为控制方程组的无量纲参考长度、速度和压力，并引入无因次流函数和涡量： ，， 图1 物理模型 (1.1) (1.2) (1.3) 边界条件如下： ： (1.4a-c) ：； (1.5a-c) ： (1.6a-c) ： (1.7a) (1.7b) (1.7c) (1.7d) 此外，由于进出环形池的流体净质量流量为零，流体应满足方程 (1.8) 另，流体为不可压缩流体，体积恒定，故满足条件： (1.9) 式中，p为无因次压力，为无因次温度。毛细雷诺数，Marangoni数，毛细数，表征自由表面变形的程度，取。Biot数，μ为动力黏度，为表面张力温度系数，α为热扩散率，ν为运动黏度，λa、λl 在浅液池内，，。为简单起见，将r用x替代，即, 。 2 中心区域的解 在中心区域，作如下变形：，，，，，，各物理量按下式展开：，并代入到变形后的控制方程和边界条件中，得到各级方程组和边界条件。 零级近似控制方程组及边界条件如下： ，， (2.1a-c) ： (2.2a-c) ：； (2.3a,b) ；