板式蒸发式冷凝器两相降膜流动CFD模拟及传热研究.docx

1、板式蒸发式冷凝器两相降膜流动CFD模拟及传热研究华南理工大学学报(自然科学版第36卷第7期Journal of Sou th C hina U n iversity of TechnologyV ol .36N o .72008年7月(N atu ral Science Edition July 2008文章编号:10002565X (20080720006205 收稿日期:20072092103基金项目:广东省自然科学基金重点资助项目(04015950作者简介:朱冬生(19642,男,博士,教授,主要从事热能工程及计算流体力学的研究.E 2mail:cedshzhuscut .edu .cn

2、板式蒸发式冷凝器两相降膜流动CF D 模拟及传热研究3朱冬生张景卫吴治将李元希蒋翔(华南理工大学传热强化与过程节能教育部重点实验室,广东广州510640摘要:建立了气-液两相流二维计算流体力学(CF D 模型,运用Fluent 软件对板式蒸发式冷凝器板束中气-液两相的逆流、并流两种操作进行了模拟,直观地表征了板束中喷淋水流量、风速及风向对水膜流动的影响,并用水蒸气对两种操作进行了传热实验研究.模拟分析得出的结论与实验结果吻合甚好,表明空气与水并流比逆流更有利于利用液体薄膜强化传热的特性.关键词:冷凝器;CF D 模拟;强化传热;降膜;板束中图分类号:T B 61+1;T Q 051.6+1文献

3、标识码:A板式蒸发式冷凝器是利用薄膜传热传质理论和蒸发冷凝传热理论研发出来的一种高效节能节水冷换设备.目前,国内外应用及研究的蒸发式冷凝器的核心部件是换热盘管,其研究主要集中在管型的设计、改进及管表面的处理上123,而关于板式蒸发式冷凝器的研究和报道较少425.板式蒸发式冷凝器中传热传质过程的发生与水膜和空气的流动是同时进行的,它们之间存在着相互影响、相互耦合的关系,对其交换机理的研究是提高板式蒸发式冷凝器热工性能的关键.因此,要研究板式蒸发式冷凝器的传热传质机理,水膜与空气的流动形式及其特点是不可回避的问题,关于下降水膜与空气流动的理论研究已有不少文献进行了报道,采用数值模拟的方法进行研究的

4、报道较少6.文中对板式蒸发式冷凝器的气-液两相降膜流动进行了模拟,分析了喷淋水量、空气流速及风向对水膜流动的影响,并用水蒸气进行了传热实验研究,为板式蒸发式冷凝器的优化设计和传热传质机理的研究提供了重要依据.1模拟及分析文中研究的板式蒸发式冷凝器所用的传热板型为板管式网格板,剖面图与波纹板类似.文中在合理简化实物模型的基础上建立了板式蒸发式冷凝器空气-水两相流流动的二维计算模型,如图1所示.图1计算物理模型及其边界条件Fig .1Physical model and its boundary conditi on模型相对比较简单,因此,采用四边形结构化贴 体网格.在划分网格时,采用垂直于板面方

5、向从气相区向液相区逐层加密的方法.计算采用RNG k 2两方程湍流模型;离散时时间项采用隐格式;对流项采用一阶迎风格式;压力项采用PREST O !算法;压力-速度耦合方程的求解采用P I S O 方法;时间步长在10-6到10-4之间选取.气-液相界面的追踪采用精度较高的Geo 2Reconstruct 界面重构技术7.动量源项由F LUE NT612中编制UDF 接口程序来实现.1.1喷淋水量的影响不同水流量下水膜的流动分布情况如图2所示.从图2可以看出,水膜沿流动方向厚度逐渐变薄,且随着水流量的增大水膜厚度不断增加 .图2不同水流量下水膜的流动分布Fig .2W ater fil m d

6、istributi on at different water fl owrates1.2气相对水膜流动的影响1.2.1气液两相逆流水流量W =0105kg/s 、空气流速u =210m /s 时水膜的流动情况如图3所示 .图3气-液两相逆流时水膜的流动分布Fig .3W ater fil m distributi on in air 2water reverse fl ow从图3可以看出,板面上水膜随板形状形成周期性波动,进风后,由于来流剪切力的作用,水膜受到托滞,1157s 时,水膜明显发生团聚,厚度增大;1174s 时,水膜平均厚度较初始时有明显的增大;10124s 时,水膜自由表面形状

7、与初始时、1142s 时的流动状况相似,水膜流动状态呈周期性变化.水流量W =0105kg/s 、空气流速u =213,215,218,310m /s 时水膜的流动情况如图4所示,气液两相逆流液膜的厚度分布情况如图5所示.结合图3、图4可以看出,随着空气流速的增大,水膜表面的湍动程度逐渐加剧,空气流速从210m /s 增大至215m /s 时,水膜平均厚度比初始状态时均有所增大,且随着空气流速的增大,液膜厚度越来越厚,见图5.当空气流速达到218m /s 时,由于水膜受强大的剪切力作用,水膜开始出现大的液团,液团在下落的过程中被高速吹来的空气吹散,形成液滴, 水图4气-液两相逆流空气流速对水膜

8、流动分布的影响Fig .4Effects of wind vel ocity on distributi on of water fil m in air 2water reverse fl ow7第7期朱冬生等:板式蒸发式冷凝器两相降膜流动CF D 模拟及传热研究 图5气-液两相逆流液膜厚度Fig .5Thickness of water fil m of air 2water reverse fl ow膜由连续相变成分散相;空气流速达到310m /s 时,水膜被空气大量吹脱,板表面几乎完全裸露出来.1.2.2气液两相并流水流量W =0105kg/s,空气流速u =210,213,215,2

9、18m /s 时水膜的流动情况如图6所示.图7为气液两相并流液膜的厚度分布情况.由图6可见,空气与水并流时,由于剪切力作用,使得流速较大的空气带动水膜快速下落,且掠带一部分水膜,水膜的厚度较初始时减薄了,随着风速的不断加大,空气掠带的水团也不断加大,水膜厚度大大减薄 .图6气-液两相并流空气流速对水膜流动分布的影响Fig .6Effects of wind vel ocity on distributi on of water fil m in air 2water parallel fl ow图7气-液两相并流液膜厚度Fig .7Thickness of water fil m of air

10、 2water parallel fl ow2实验研究及分析在空气流速为315m /s 、其它条件相同的情况下,空气与水两相逆流、并流时热流密度q f 随喷淋密度变化的情况如图8所示.由图8可见,喷淋密度较小时,逆流时的热流密度比并流时稍大.因为喷淋密度较小时, 板表面未被图8热流密度随冷却水喷淋密度的变化Fig .8Effects of the s p ray density of the cooling water on theheat strea m density水膜完全覆盖,局部有“干斑”,在板内有热流体传热时,板外水膜受热蒸发,“干斑”面积扩大,在这种情况下,逆流操作使下降水膜受托

11、滞,有利于布膜;而并流操作时,水膜在风的吹送下会变成“溪流”,在蒸发作用下,加剧了“溪流”的生成,更不利于液膜分布.所以,在喷淋密度较小时出现了逆流时的热8华南理工大学学报(自然科学版第36卷流密度比并流时稍大的现象.但在板式蒸发式冷凝器中,板表面要完全被水膜覆盖,才能充分利用蒸发冷却的作用.随着喷淋密度的增大,板表面逐渐被水膜完全覆盖,此时,并流时的热流密度逐渐比逆流时的大.在喷淋密度达到01196kg/(m s 时,并流、逆流的热流密度均随喷淋密度的增大而开始减小,这是因为随着喷淋水量的不断增大,水膜厚度增大,热阻变大,但并流时的热流密度仍比逆流时大.在喷淋密度为01175kg/(m s

12、、其它条件相同的情况下,空气与水两相逆流、并流时热流密度q f 随空气流速u 变化的情况如图9所示 .图9热流密度随空气流速的变化Fig .9Effects of the air vel ocity on the heat strea m density由图9可见,相同的喷淋密度下,并流时的热流密度随空气流速的增大而不断增大,而逆流时的热流密度随空气流速的增大先减小而后又略增大,然后继续减小,并流时的热流密度仍比逆流时大.这是因为喷洒在冷凝板上的水紧贴板壁而下,是典型的连续性降膜,采用并流操作,水膜受剪切力的作用,被空气快速吹下,且掠带一部分水膜,膜厚变小,热阻减小;采用逆流操作,水膜受空气的

13、托滞作用,膜厚变大,热阻增大,空气流速增大时,水膜很容易被“吹脱”,削弱了换热强度.3结论对板式蒸发式冷凝器气液两相降膜流动的CF D 模拟与传热实验表明,模拟分析得出的结论与实验结果吻合甚好.在板式蒸发式冷凝器中,逆流操作在空气流速较小时,水膜厚度增大,热阻变大,空气流速较大时,水膜容易被“吹脱”,使板壁裸露出来,削弱了换热强度;并流操作时水膜厚度随空气流速的增大而减小,热阻也大大减小,强化了换热强度.因此,在板式蒸发式冷凝器中,空气与水并流比逆流更有利于利用液体薄膜强化传热的特性.参考文献:1Kachh waha S S,Dhar P L,Kale S R.Experi m ental s

14、tudiesand nu merical si m ulati on of evaporative cooling of air with a water s p ray J .I nternati onal Journal HeatMass Trans 2fer,1998,41(2:4472474.2J iang Xiang,Zhu Dong 2sheng .Heat transfer enhancementin evaporative condenser C Pr oceedings of the 3rd I n 2ternati onal Sy mposiu m on Heat Tran

15、sfer and Energy Con 2servati on .Guangzhou:South China University of Technol o 2gy Press,2004:102521031.3蒋翔,朱冬生,唐广栋.蒸发式冷凝器管外水膜与空气传热性能与机理的研究J .流体机械,2006,34(8:59262.J iang Xiang,Zhu Dong 2sheng,Tang Guang 2dong .I nvesti 2gati on on heat transfer and resistance perfor mance outside tubes of evaporativ

16、e condenser J .Fluid Machinery,2006,34(8:59262.4周孝清,陈沛霖.板式间接蒸发冷却器的热工计算方法J .华南建设学院西院学报,1998,6(2:9215.Zhou Xiao 2qing,Chen Pei 2lin .The ther mal analysis of the indirect evaporative cooler J .Journal of S outh China Constructi on University:W est Ca mpus,1998,6(2:9215.5张景卫,朱冬生,蒋翔,等.蒸发式冷凝器及其传热分析J .化工机

17、械,2007,34(1:1102114.Zhang J ing 2wei,Zhu Dong 2sheng,J iang Xiang,et al .Evapora 2tive condensers and their heat transfer analysis J .Che m i 2cal Engineering &Machinery,2007,34(1:1102114.6谷芳,刘春江,余黎明,等.气-液两相降膜流动及传质过程的CF D 研究J .高校化学工程学报,2005,19(4:4382444.Gu Fang,L iu Chun 2jiang,Yu L i 2m ing,et al .

18、The CF D si m ulati on of mass 2transfer p r ocess in falling fil m with countercurrent t w o 2phase fl ow J .Journal of Che m ical Engineering of Chinese Universities,2005,19(4:4382444.7H irt C W ,N ichols B D.Volu me of fluid (VOF methodfor the dyna m ics of free boundaries J .Journal of Com 2puta

19、ti onal Physics,1981,39:2012225.9第7期朱冬生等:板式蒸发式冷凝器两相降膜流动CF D 模拟及传热研究01华南理工大学学报(自然科学版第36卷CF D Si m ul ati on and I nvesti gati on i n to Heat Transfer for Falli n g Fil m with Two2Phase Flow i n Pl ate2Type Evaporati ve CondenserZhu D ong2shengZhang J ing2w eiW u Zhi2jiangL i Yuan2x iJ iang X iang(Key

20、 Laborat ory of Enhanced Heat Transfer and Energy Conservati on of the M inistry of Educati on,South China University of Technol ogy,Guangzhou510640,Guangdong,ChinaAbstract:I n this paper,a t w o2di m ensi on computati onal fluid dyna m ics(CF Dmodel of the gas2liquid t w o2phase fl o w was establis

21、hed,and the countercurrent and parallel operati ons of the t w o2phase fl ow in the p late bundles of a p late2type evaporative condenser were si m ulated by means of the Fluent s oft w are,thus directly illustrating the influ2 ences of the water s p ray fl ow,wind vel ocity and wind directi on in t

22、he p late bundles on the fl ow of water fil m.An experi m ental investigati on int o the heat transfer of the t w o operati ons was then perfor med,with water vapor as the heat transfer medium.It is shown that the si m ulated results accord well with the experi m ental ones,meaning that the parallel

23、 fl ow of air and water is superi or t o the countercurrent fl ow in making use of enhanced heat transfer of fluid fil m.Key words:condenser;CF D si m ulati on;enhanced heat transfer;falling fil m;p late bundle“表面热功能结构数字制造及沟槽式微热管关键技术与装备”通过鉴定华南理工大学机械与汽车工程学院汤勇教授主持完成的项目表面热功能结构数字制造及沟槽式微热管关键技术与装备于2008年4月

24、7日通过由广东省科技厅组织、主持的科技成果鉴定。该项目以解决在微电子领域具有普遍意义的高热流密度问题为切入点,以揭示复杂表面热功能结构多尺度和多维数效应的科学规律为突破口,建立复杂表面热功能结构相关热物理、数字化设计及制造理论体系;通过对复杂表面拓扑结构、几何形状及结构参数进行数字化描述,建立了复杂表面热功能结构的定量分析和数字设计的理论,为表面结构拓扑优化设计理论提供了有利支撑;通过对亚结构和微结构典型的犁削、拉削、挤压、旋压等关键技术和加工机理的研究,建立了表面亚结构和伴随生成的微结构滑移和生成过程模型,确定了加工工艺参数对亚结构和微结构的影响规律。课题组在机理研究的基础上,结合表面结构数字化制造技术,提出沟槽式薄壁微热管高速充液旋压成形技术;深入研究分析了微热管的内沟槽充液旋压成形机理,确定了适合于产业化规模生产的工艺参数对沟槽及其亚结构成形的影响规律,成功研制出高性能沟槽式薄壁微热管,为我国微电子散热领域的关键元件制造技术提供了技术理论支撑。于工业生产中利用该项技术研制的微热管具有成本低、重量轻、响应快等优点。鉴定委员会专家经过认真讨论,认为该成果具有自主知识产权,达到国际先进水平,一致同意通过鉴定。