液膜沿非均匀加热斜面流动传热稳定性.pdfVIP

[收稿日期] 2010- 03- 26 [基金项目] 国家自然基金(资助 [作者简介] 樊小朝( 1979- ) ,男, 华北电力大学硕士研究生,从事热能技术应用研究。 液膜沿非均匀加热斜面流动传热稳定性 樊小朝1,史瑞静2 ( 1. 华北电力大学,河北 保定 071003; 2. 保定科技职业学院,河北 保定 071000) 摘 要 从原始的 N S 方程推导出相应的小扰动方程 ,通过长波摄动法进行了求解,得到了一系列稳定参数在 内的临界雷诺数关系式,然后用 MATLAB 软件数值模拟分析了 Peclet 数、斜面角度和 Magelanni数等对流动稳定 性的影响。 关键词 液膜; 非均匀加热;稳定性; Magelanni数 Stability of the Flow and Heat Transfer of Non Uniform Heated Liquid Film Along the Slope FAN X iao chao 1 , SH I Rui j ing 2 Abstract This paper is derived from the original N S equations and the corresponding equations of small perturbat ions are obtained. T hen t hrough the long wave perturbation method for solving, the cr itical Reynolds number relationship in cluding a ser ies of steady parameters is got. Then the Peclet number , slope angle and Magelanni number effect on the flow stability ar e numerical simulation surveyed by MATLAB softw ar e. Key words liquid film; non uniform heating; stability; Magelanni number !中图分类号? O357 !文献标识码? A !文章编号? 1674- 3229( 2010) 03- 0035- 05 0 引言 液膜流动是工业中较为常见的一种两相流,如 存在于锅炉水冷壁、凝汽器和核电站的冷却系 统 [ 1] 。液膜在自由下降过程中,受小扰动因素的影 响, 将引起流动结构和状态的改变, 导致表面出现 波动。如果引发流动不稳定的小扰动趋于减小,则 流动是稳定的;反之,则处于不稳定状态; 如果扰动 始终不变,则流动处于中性稳定状态 [ 2- 3] 。在工业 应用中经常会遇到衬底加热的情形,平板加热对流 动稳定性产生的影响是我们极想了解的。当一层 水平薄膜底部被加热时, 有两个主要的因素将会影 响流体的稳定性,即浮力和表面张力。当薄膜厚度 很薄时,由表面张力引起的热毛细机制将占主导地 位,经典的贝纳德对流实验事实上就是由热毛细机 制引起的。因此, 我们将前人的研究推广到底部加 热的情况。 在对液膜流动稳定性的研究中, 目前的理论研 究方法分为线性和非线性两种,且流态都局限于小 雷诺数的层流区。并且线性理论得到广泛应用,它 最早开始于 Benjam in [ 4] 和 Yih [ 5, 6] 对等温降膜的线 性稳定性。 定性分析方面, Yih首次提出了临界雷诺数的 概念。对加热液膜流动稳定性的研究上, 印度 Arghya Samanta在 2008 年研究了非惯性液膜沿均 匀加热斜面下滑的稳定性, 他用正则模法分析研究 了下降液膜的线性稳定性 [ 7] , 同年又研究了液膜沿 垂直的非均匀加热斜面流动稳定性 [ 8] 。文中假定壁 面温度是线性变化的,并利用长波摄动法得到非线 性演化方程,然后使用多尺度的方法进行了非线性 稳定分析,从而研究了行波解的演化方程, 得到随着 Peclet数的增加超临界不稳定区域将会增加, 而亚 临界稳定区域将会减少。 1 液膜流动的控制方程及其求解 液膜在重力作用下沿底部加热的倾斜平板下 落,假设周围的环境为无粘滞和静止的气体, 如图 1,平板与水平夹角为 ,液膜垂直平均厚度为 h, 环 境气体温度为 T a , 我们假设自由面温度 T s = T a , 平板上的温度分布为 T w。 #35# 2010 年 6 月 廊坊师范学院学报(自然科学版) Jun. 2010 第 10 卷第 3 期 Journal of Langfang Teachers College( Naturnal Science Edi