立式降膜闭式冷却塔设计程序研究*

陈二雄，朱冬生，涂爱民，刘世杰，陈杭生

立式降膜闭式冷却塔设计程序研究*

陈二雄1,2,3†，朱冬生1,2,3，涂爱民1,2,3，刘世杰1,2,3，陈杭生1,4

（1. 中国科学院广州能源研究所，广州 510640；2. 中国科学院可再生能源重点实验室，广州 510640；3. 广东省新能源和可再生能源研究开发与应用重点实验室，广州 510640；4. 中国科学院大学，北京 100049）

提出一种立式降膜闭式冷却塔的设计计算方法，按照此方法使用Visual Basic 6.0开发工具设计计算程序，并对程序流程进行详细说明。结果显示，该计算程序适用于立式扭曲管中不同工艺流体（冷却水、喷淋水和空气）进行换热的各种工况。通过对比程序计算和手工计算结果，发现二者吻合程度较好，因此通过计算程序能够提高开发人员的设计效率和准确性，减少迭代计算工作量，方便对不同设计方案结果进行比较。后续可以衍生出针对各种发生相变过程（制冷剂、水蒸气冷凝等）的热力计算。

立式降膜；闭式冷却塔；扭曲管；设计程序

0 引 言

目前在闭式冷却塔设计计算中，还没有简单易用且计算结果准确度高的方法，相关行业设计人员在设计闭式冷却塔的过程中更多的是采用依靠工作经验的方法进行设计，较少一部分设计人员采用的是Excel表格计算方法，此方法需要使用者对表格非常熟悉，并需要多次计算并调整试算初始值，直到假设值和计算结果相吻合方能完成，此过程将消耗大量的时间和人力。此次，笔者考虑尝试编制程序，利用计算机处理试算及迭代的过程，以方便对不同结构参数闭式冷却塔进行设计，提高设计效率，节约人力和降低出错的概率[1]。

利用计算机程序设计空调换热器早有应用，王玉珏等[1]开发了一种管壳式冷凝器计算程序，适用于多种管型和制冷剂，设计算例和手工计算结果表明吻合程度较好。方运惠[2]利用Visual Basic 6.0工具开发套管式、折流杆式和强化换热器设计辅助软件，通过计算程序代替复杂的设计计算，大大减少了迭代计算时间，能够按照工艺条件快速地设计出符合要求的热交换器。

国内外学者对扭曲管应用在换热设备的研究已进行多年，刘世杰等[3]利用FLUENT软件，以水为介质，在雷诺数13 000 << 110 000条件下，研究了机械式蒸汽再压缩（mechanical vapor recompression,MVR）蒸发器内扭曲管的不同结构参数对管内传热与压降性能的影响，拟合出了计算准则关联式。朱冬生等[4]对扭曲管管内湍流换热过程进行了实验和数值研究，采用场协同理论分析了管内强化传热机理，并获得准则关联式。

DZYUBENKO[5-6]对扭曲管管内和管束间传热传质与流阻性能展开了研究，并拟合出关联式。

基于提高闭式冷却塔设计效率，节约人力和降低出错的概率等目的，同时结合国内外学者对扭曲管热力性能计算准则关联式的相关研究，利用Visual Basic语言编制一个用于计算立式降膜闭式冷却塔的计算程序，本文将从计算程序算法的角度对该程序进行分析和介绍，并对工程实例中采用传统方法计算和本程序计算的结果进行对比，分析该程序的可靠性和有效性。

1 程序设计目标

本程序适用于立式降膜闭式冷却塔结构设计，主要针对扭曲管，管内冷却介质为水，管外介质为水和空气，管内冷却水先把热量传递给管外壁液膜，液膜吸热升温，最后液膜把热量传递给空气，空气被风机排到机组外面。

程序在输入必要参数后能够自动进行闭式冷却塔热力性能计算，在满足喷淋水流量、风机风量以及换热器结构等参数后，确定闭式冷却塔整体结构。

2 设计计算思想

（1）目前，闭式冷却塔设计软件不多，主要原因是生产企业都是大多依据经验进行设计，且光滑圆管热力计算公式较为成熟准确。而本程序主要用于扭曲管立式降膜闭式冷却塔设计，通过选用已公开证明有效的换热模型，较为准确地对立式降膜闭式冷却塔进行热力性能计算。

（2）针对不同结构参数的扭曲管，利用相同的换热模式进行热力性能计算，通过合理安排换热器管排数量和单排管数目，得到最佳换热器结构设计。另外，通过计算管外壁换热系数确定喷淋水和空气流量和压降，确定喷淋水泵和风机的参数。

（3）目前商用闭式冷却塔大多采用光滑不锈钢圆管作为换热管，换热系数较扭曲管差，因此采用扭曲管替代光滑圆管有利于减少换热器重量和体积，但是目前的经验对于扭曲管立式降膜闭式冷却塔设计不具有参考价值，需要重新设计一款适用于其结构的软件。

为了计算的简便，本文依据迈克尔焓差理论，假设如下[7]：①换热盘管的管外表面被充分润湿，传热传质界面相同；②认为盘管外喷淋水膜的表面温度和内部温度一致；③所有物性参数为常量，不随空气湿度和水温发生变化；④忽略辐射传热量；⑤出口空气近似认为是饱和空气；⑥忽略水蒸发造成喷淋水质量变化的影响。

3 程序流程

3.1 输入参数

3.1.1 基本参数

输入数据参数包括结构数据、工艺数据和物性数据三大类。结构数据包括扭曲管结构参数、盘管结构参数；工艺参数包括冷却水、喷淋水和空气的温度、流量等参数；物性参数包括冷却水和喷淋水物性参数。

图1为立式降膜闭式冷却塔设计软件的主界面，从图中可以看出，主要包括工艺参数设定、冷却水物性参数设定、喷淋水物性参数设定、扭曲管结构参数设定、盘管结构设定、空气换热量计算结果和计算结果输出几块内容。其中计算结果输出为本设计软件的最后结果。

图1 软件主界面

3.1.2 结构数据参数输入

扭曲管结构参数包括：换热管材质，加工扭曲管所用圆管外径、壁厚，扭曲管管外长轴、管外短轴、扭矩、管长，管材导热系数和管材密度等。

扭曲管是通过对普通圆管的表面进行特殊加工使管内流体扰流程度加强来实现管内强化传热的技术，扭曲管的横截面呈椭圆形状，其长轴和短轴尺寸参数对整个换热器的热力性能有着重要的影响[8]。

盘管结构参数包括换热管片的数量、单片换热管根数和管片间距。换热管片由若干根扭曲管通过弯头连接组成蛇形盘管的结构，管片之间存在一定的间距供喷淋水和空气流通。

3.1.3 工艺数据参数输入

工艺参数包括冷却水、喷淋水和空气的温度、流量等，用户根据实际情况，输入要求值。

其中冷却水和空气的工艺参数所需输入的内容如图2所示，包括冷却水换热量、进口温度和体积流量，空气流量和进风干湿球温度。

图2 参数设置界面

3.2 程序计算流程

主流程见图3。输入必要参数后，程序将首先根据冷却水入口温度计算出口温度及对数平均温差[9]。

图3 计算程序流程图

冷却水出口温度由下式计算：

根据环境干湿球温度和冷却水进口温度，假定一个喷淋水的进口温度（比环境湿球温度高2 ～ 5℃），通过设定喷淋水进出口温差为5℃得到喷淋水出水温度和流量。

在得到上述冷却水和喷淋水进出口温度后计算换热器的对数平均温差：

程序开始前输入已知设计参数，设定闭式冷却塔换热器的结构参数，包括换热管片数量、单片换热管数量和管片间距。先计算管内外传热系数，最后得到总传热系数，根据换热面积和对数平均温差，计算出总换热量，与设计换热量进行比较后，当满足条件（考虑设计余量，范围应在10% ～ 20%）则结束程序并输出计算结果，否则重新调整换热器结构参数并再次运行计算程序（图3）。

3.3 换热系数计算

3.3.1 管内对流换热系数计算

管内冷却水与管壁之间的对流传热系数的计算关系式采用最广泛的Dittus—Boelter关系式[10]：

考虑到扭曲管比普通光滑圆管换热系数高1.5 ～ 2倍，因此上述换热系数计算结果需要乘以合适的系数，系数大小与扭曲管长短轴比有关。

雷诺数计算公式为：

扭曲管管内当量直径计算公式为：

3.3.2 管内压降计算

换热器管内压降主要包括两部分：直管内摩擦压降和各程回弯处的局部阻力降[11]。

（1）扭曲管直管内摩擦阻力引起的压降

（2）局部阻力降

换热器回弯头处局部阻力降通过以下公式计算[10]：

考虑到扭曲管比普通光滑圆管压降大1.2 ～ 1.8倍，因此上述压降计算结果需要乘以合适的系数，系数大小与扭曲管长短轴比有关。

3.3.3 液膜侧界膜换热系数计算

闭式冷却塔喷淋水沿着扭曲管管外传热面流动，形成液膜。沿垂直传热管壁下流的液膜侧界膜导热系数的关联式较多，这里使用的是精度较好的威尔克实验公式[12]。

（1）过渡区的平均界膜导热系数

其中，液膜雷诺数计算公式为：

式中：为管外降膜单位宽度的流量，kg/(m∙h)。

考虑到闭式冷却塔喷淋水在使用时不会全部覆盖到换热管表面，上式在计算时需要乘以一个系数，这里取0.75 ～ 0.95。

（2）流动发达区平均界膜导热系数

当ii<<1 600时：

当1 600 << 3 200时：

当> 3 200时：

当> 1 600时：

当< 1 600时：

另外，ii由下式求得：

4 结果与讨论

为了比较使用Excel表格和VB 6.0程序两种方法计算得到的结果差异，针对其中某个案例分别进行设计，最后进行比较。

本闭式冷却塔设计案例参数如下：换热量150 kW，进水温度45℃，空气流量为30 000 m3/h，进风干、湿球温度分别为32℃、28℃，采用的换热圆管外径16 mm、壁厚1 mm、材质为316L。由以上条件试设计闭式冷却塔换热器。

从表1可以看出，根据VB 6.0程序计算出来的结果与通过Excel表格计算结果对比可知，程序计算总传热系数结果比Excel表格计算结果小2.5%，换热量盈余小11.3%，误差均较小。其余计算结果偏差绝对值均在10%以内。造成偏差的主要原因是迭代计算过程中，手工计算存在输入数值偏差，而计算机程序不存在这个问题。

由此可知，使用Excel表格和VB 6.0程序两种方法计算得到的结果非常接近，基本可以确定能用VB6.0程序替代常规Excel表格方法用于设计闭式冷却塔。

表1 两种方法计算结果比较

5 结束语

采用“编程语言+换热模型”结合的方法开发的闭式冷却塔热力计算程序能够较为准确地计算各种实际工况下的热力性能，能够判断设计方案是否满足实际应用需求，对于大部分使用者，程序易于操作，只要输入设计参数就能马上得到相应的计算结果。另外，可以在编程过程中加入多个不同的换热模型，并将计算结果进行比较分析。

最后，由于目前公开的换热模型均存在计算偏差，与实际情况比较会有一定的误差范围，但是这一误差可以通过大量的工程实验数据对换热模型进行修正，使其计算结果与实际工况相近。

[1] 王玉珏, 黎立新, 季建刚, 等. 管壳式冷凝器设计程序研究[J]. 制冷与空调, 2008, 8(3): 75-78. DOI: 10.3969/ j.issn.1009-8402.2008.03.021.

[2] 方运惠. 基于VB6.0的换热器辅助程序设计[D]. 中国石油大学(华东), 2019.

[3] 刘世杰, 朱冬生, 张洁娜, 等. 三维变形管MVR蒸发器的强化传热数值研究[J]. 高校化学工程学报, 2018, 32(6): 1307-1313. DOI: 10.3969/j.issn.1003-9015.2018. 06.009.

[4] 朱冬生, 石仲璟, 钱泰磊, 等. 扭曲椭圆管换热器的数值模拟及场协同分析[J]. 高校化学工程学报, 2015, 29(1): 64-71. DOI: 10.3969/j.issn.1003-9015.2015.01.009.

[5] DZYUBENKO B V. Influence of flow twisting on convective heat transfer in banks of twisted tubes[J]. Heat transfer research, 2005, 36(6): 449-460. DOI: 10.1615/HeatTransRes.v36.i6.20.

[6] DZYUBENKO B V. Estimation of the thermohydraulic efficiency of heat exchanging apparatuses with twisted tubes[J]. Heat transfer research, 2006, 37(4): 349-363. DOI: 10.1615/HeatTransRes.v37.i4.50.

[7] 张丽霞. 带填料的复合流闭式冷却塔的结构设计与性能分析[D]. 济南: 山东建筑大学, 2018: 11.

[8] 朱冬生, 安冬旭, 李霞, 等. 高效复合强化换热器的管程性能[J]. 化工学报, 2014, 65(2): 453-459. DOI: 10.3969/j.issn.0438-1157.2014.02.012.

[9] 沈维道, 童钧耕. 工程热力学[M]. 4版. 北京: 高等教育出版社, 2007: 393-395.

[10] 杨世铭, 陶文铨. 传热学[M]. 4版. 北京: 高等教育出版社, 2006: 133-134.

[11] 兰州石油机械研究所. 换热器[M]. 2版. 北京: 中国石化出版社, 2013: 33.

[12] 尾花英朗. 交换器设计手册(下册)[M]. 徐中权, 译. 北京: 石油工业出版社, 1984: 379-382.

Design Procedure of Vertical Falling Film Closed Cooling Tower

CHEN Er-xiong1,2,3, ZHU Dong-sheng1,2,3, TU Ai-min1,2,3, LIU Shi-jie1,2,3, CHEN Hang-sheng1,4

(1.Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China;2. CAS Key Laboratory of Renewable Energy, Guangzhou 510640, China;3. Guangdong Provincial Key Laboratory of New and Renewable Energy Research and Development, Guangzhou 510640, China;4.University of Chinese Academy of Sciences, Beijing 100049, China)

A computational design method for the vertical falling film closed cooling tower was proposed, and a corresponding calculation program was designed by Visual Basic 6.0. Results showed that, the calculation program was applicable to various heat transfer conditions of different process fluids (cooling water, spray water and air) in vertical twisted pipes. The program output fits well with the manual computation. Therefore, the calculation program can improve the design efficiency and accuracy of developers, reduce the workload of iterative calculation, and facilitate the comparison of the results of different design schemes. Moreover, it is friendly to beginners. Thermodynamic calculations can then be derived for various phase transition processes (refrigerant, water vapor condensation, etc.).

vertical falling film; closed cooling tower; twisted tube; program design

2095-560X（2021）05-0411-07

TK-9

A

10.3969/j.issn.2095-560X.2021.05.007

陈二雄（1986-），男，硕士，工程师，主要从事强化传热和节能环保方向研究。