垂直表面自由液体降膜流动表面波的数值模拟研究.pdfVIP

中国工程热物理学会 论文编号:983110 传热传质学学术会议 垂直表面自由液体降膜流动表面波的数值模拟* 叶学民 阎维平 蒋章焰*高正阳 (华 北 电 力 大 学 071003) 中(中科院工程热物理所 髯100080) 摘要 对垂直自由降膜流动表面波建立了数学模型及计算方法，根据实验测定的 在 等温壁面条件下，Re=880时，随时间变化的薄膜厚度，在运动坐标系下，通过采用有 限差分法求解动量方程，连续性方程以及相应的自由界面条件，得到了波动层流液膜 内的速度分布以及压力分布。 关键词:垂直降膜流 速度 压力 表面波 自由界面 前言 液体在重力作用下以薄层形式沿垂直管壁向下流动，这种流动形态称为垂直自由 降膜流动。液体薄膜流具有小流量、小温差、高换热系数、高热流密度、结构简单、 动力消耗小等独特优点，因此，业已作为一项高效传热传质技术在化工、轻工、能 源、航天、石油化工、电子工业等许多工业领域得到广泛应用。譬如，膜蒸发器、膜 吸收器、热管和其他许多两相流和热传递过程中。特别值得注意的是，应用液体薄膜 的独特优点来解决高技术领域中遇到的高热流密度下的热强化问题，越来越引起人们 的关注，这方面的实例有:大规模集成电路的薄膜冷却、第二代核电站安全壳薄膜蒸 发冷却方案以及未来空间动力站应用的热管辐射器、液滴辐射器 (LDR)、新型太阳能 集热器等等。 根据以往对于垂直表面自由降膜的水动力特性的研究，一般认为:Re20，流动为 层流，膜表面是平滑的;20Re1000，流动是具有波动的层流，Re1000}2000,流动 是具有波动的湍流。在工业感兴趣的雷诺数范围内，降膜呈现出非常不规则的波动表 面。图(1)给出了实测的液膜厚度随时间变化的曲线[1]。从图上可以看出，在层流底 层上滑动着大大小小的波团。对于波幅是其层流底层厚度2-5倍的大波，其携带着大 部分流动质量，它对波内、波与壁面、波与外界的传热传质速率，起明显的控制作用 幻【。因此求解膜内的速度分布，压力分布，对于深入分析和描述与垂直降膜有关的各 种传热和传质规律是非常必要的。 由于通过实验测定膜内的速度分布极其困难，因而通过数学模拟获得结果显得极 为重要。已有的大多数分析模型是建立在由Kapitza(1964)发展的速度呈抛物线分布基 *国家自然科学基金资助项 目，谨此致谢. 珊一 157 除了在远低于工业感兴趣的流量下，计算结果并不能准确地 础之上。这些计算模型， Maroneta1(1985)对于孤立波建立分段模型，该模型 代表己测得的薄膜的流动特征。 一-一 一一 、 “. “， 上分 甲行 、山 .尹 占“ 刁 吧比 ，八之 1.1 也不 能解释孤 立炙 菠 幅利 胜仅 幽派 省蟹% 计算结果仅与有限的实验数据吻合，也不能解释A1aaTMTUK1NuliM}a}ZIw 。Kheshigi和Scriven(1987)采用有限元方法及周期性边界条件证实无穷小扰动的发展 但仅限于低雷诺数，也不能模拟充分发展的波形。F.K.Wasdenetal(1989)碘溶 液为介质，应用测得的波形计算了速度场，给出了孤立波及演化中的波的流线图仁2] 。RRUuY丫aann漪g[]ByardD.Wood(1991)采用一种简单的频谱方法求解了周期性波动膜方程 利用现有的线性稳定性分析，其求解结果，在Re50和50Re150的情况下，与实验一 致。但对于较高的Re，该方法并不接近实际情况 5〔]eM.Fang和T.Nosoko(1994)对控 制方程组进行数量级简化后，采用积分法对小雷诺数Re=12和27在动坐标下计算了速 度分布和压力分布[4]0 在己有的相关的研究工作中，分析模型不是限于低雷诺数范围，就是采用了周期 性的边界条件，这与实际情况有一定程度上的差别。而且，工业上应用更