CN109614746B - 一种蒸发式冷凝器结构设计方法 - Google Patents
一种蒸发式冷凝器结构设计方法 Download PDFInfo
- Publication number
- CN109614746B CN109614746B CN201811622756.4A CN201811622756A CN109614746B CN 109614746 B CN109614746 B CN 109614746B CN 201811622756 A CN201811622756 A CN 201811622756A CN 109614746 B CN109614746 B CN 109614746B
- Authority
- CN
- China
- Prior art keywords
- air
- heat exchange
- tube
- heat transfer
- exchange tube
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Air Conditioning Control Device (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
本发明公开了一种蒸发式冷凝器结构设计方法,满足冷水机组用蒸发式冷凝器的结构快速设计,本发明通过自主迭代计算,得到蒸发式冷凝器目标冷量下的需要换热管面积,并可以通过调整输入参数进行结构优化设计。
Description
技术领域
本发明涉及一种蒸发式冷凝器结构设计方法,属于冷水塔设计领域。
背景技术
目前工业领域对蒸发式冷凝器的设计主要靠工程经验或者手动计算,由于计算过程中涉及大量的迭代计算,因而计算效率较低。手动计算容易带来计算错误,精确性得不到保证。另外,手动计算或者工程经验无法考察各结构参数对产品冷量大小的影响,无法进行优化设计。
故,需要一种新的技术方案以解决上述技术问题。
发明内容
为解决现有技术中存在的问题,本发明提供了一种蒸发式冷凝器结构设计方法。
为了解决上述技术问题,本发明提供了一种蒸发式冷凝器结构设计方法,所述蒸发式冷凝器包括出风口、进风口和换热管,包括以下步骤:
S1:给定蒸发式冷凝器需要的冷却负荷、冷凝温度和所处环境的进风口干湿球温度;
S2:初定结构尺寸,包括换热管内外径、换热管管壁厚度、换热管与换热管之间的间距,输入风量和喷水流量;
S3:设定初始水膜温度;
S4:计算空气、水的各点的状态参数,所述状态参数包括空气焓值、空气含湿量和温度;
S5:计算水膜与空气的对流传热系数,初算换热面积;
S6:计算管外与喷水对流传热系数和计算制冷剂在管内对流传热系数,并根据使用情况估算热阻;
S7:计算总传热系数和传热面积;
S8:比较S7中计算得到的传热面积与S5中的初算面积是否相等,若误差小于设定值,则输出需要的管排数,若误差大于等于设定值,则重新假设水膜温度,进行迭代计算;
S9:输出最终的蒸发式冷凝器结果。
进一步的,所述S4中的空气各点的状态参数包括:进风口空气的焓值、进风口空气的含湿量、出风口空气的焓值、出风口空气的含湿量、换热管外空气的平均焓值、换热管外空气的平均含湿量;
所述水各点的状态参数包括:水膜处空气的焓值和水膜处空气的含湿量。
进一步的,所述S5具体包括以下步骤:
S5-1:根据输入风量和初定结构尺寸,计算迎面风速v:
v=G/3600/((L-0.08)·(W-0.08))
式中,G为输入风量,L为长,W为宽;
根据式5计算最窄面处空气流速:
vmax=s/(s-do)·v (5)
式中,s为换热管与换热管之间的间距;d0为换热管外径;v为迎面风速;
S5-2:计算水膜与空气的对流传热系数αwa:
式中;λm为空气平均导热系数,d0为管外径,vmax为最窄面处空气流速,υm为空气平均运动粘度;
S5-3:计算换热管外空气当量对流换热系数αj:
式中,A为水膜温度修正系数,αwa为水膜与空气的对流传热系数,hw为水膜处空气的焓值,hm为换热管外空气的平均焓值,Aw为水膜和空气的接触面积,Ao为换热管外表面积,Cpm为换热管外空气的平均定压比热,tw为水膜温度,tm为换热管外空气的平均温度;
S5-5:初算换热面积A’o:
式中,Qc为机组热负荷。
进一步的,所述S6具体包括以下步骤:
S6-1:计算管外与喷水对流传热系数αw:
式中,tw为水膜温度,Γ为喷淋密度,do为换热管外径;
式中,Mw为蒸发式冷凝器水泵流量;
S6-2:计算制冷剂在管道中的对流换热系数αc·n:
式中,λ为制冷剂导热系数,g为重力加速度,r为制冷剂在该压力下的汽化潜热,μ为制冷剂在该压力下液态时的动力粘度;
S6-3:计算热阻,该热阻包括管壁热阻;
所述管壁热阻Rp根据式15可得:
Rp=δ/λ (15)
式中,δ为换热管管壁厚度,λ为管道导热系数。
进一步的,所述热阻还包括油膜热阻和污垢热阻。
进一步的,在蛇形管中对制冷剂在管道中的对流换热系数αc·n进行修正:
进一步的,所述S7具体包括以下步骤:
S7-1:计算总换热系数K:
式中,Ao为换热管外表面积,Ai为换热管内表面积,A平为换热管内外表面积平均值,Rp为管壁热阻,Roil为油膜热阻,Rfou为污垢热阻,αc·n·s为修正后的管内对流换热系数,αw为换热管外与喷水对流传热系数,αj为换热管外空气当量对流换热系数;
S7-2:计算传热面积A”o:
A"o=Qc/K/(tk-tm) (17)
式中,Qc为机组热负荷,K为总换热系数,tk为制冷机组冷凝温度,tm为换热管外空气的平均温度。
进一步的,所述S8中,管排数n根据式18得到:
式中,一排传热面积=一排中的管数×管长×管道外径周长。
进一步的,所述S8中计算得到的管排数为奇数排,修改管长,使管排数在16~24范围内。
有益效果:本发明与现有技术相比,通过本发明提供的设计方法,能够有效提高设计蒸发式冷凝器结构的效率,节省了人工成本,加快了结构优化的进度。
附图说明
图1为本发明的流程图;
图2为本发明的蒸发式冷凝器结构示意图
其中,冷风机1、收水机2、分水管3、管长4、换热管5;
图3为空气焓湿图;
图4为换热管截面示意图。
具体实施方式
本发明通过迭代计算,得到蒸发式冷凝器目标冷量下的需要换热管面积,并通过调整输入参数进行优化设计。
如图1所示的流程图,本发明具体包括以下步骤:
S1:给定蒸发式冷凝器需要的冷却负荷、冷凝温度和所处环境的进风口干湿球温度;
S2:初定结构尺寸,包括换热管内外径、换热管管壁厚度、换热管与换热管之间的间距,输入风量和喷水流量;
S3:设定初始水膜温度;
S4:计算空气、水的各点的状态参数;该状态参数包括进风口状态参数、水膜处状态参数、出风口状态参数和换热管外空气的平均状态参数;
所述进风口状态参数包括进风口空气的焓值和进风口空气的含湿量,该进风口空气的焓值h1和进风口空气的含湿量d1可由进风口干球温度t1、相对湿度和空气焓湿图可查得;
所述水膜处状态参数包括水膜处空气的焓值和水膜处空气的含湿量,该水膜处空气的焓值hw和水膜处空气的含湿量dw可由假设水膜温度和空气焓湿图查得;
所述出风口状态参数包括出风口空气的焓值h2、出风口空气的含湿量d2和出风口温度t1;
其中,出风口空气的焓值h2由式1可得:
其中,Qc为机组热负荷(机组排热量),ma为空气质量流量,h1为进风口空气的焓值;
由等热湿比线ε=Δh/Δd得:
ε=(h2-h1)/(d2-d1)=(hw-h1)/(dw-d1) (2)
式中,ε为等热湿比;hw为水膜处空气的焓值,d1为进风口空气的含湿量,d2为出风口空气的含湿量,dw为水膜处空气的含湿量;
算出出风口空气的含湿量d2,由空气焓湿图可查出出风口温度t2;
所述换热管外空气的平均状态参数包括换热管外空气的平均焓值、蒸发式冷凝器换热管外空气的平均含湿量和换热管外空气的平均温度tm:
换热管外空气的平均焓值hm由式3可得:
蒸发式冷凝器换热管外空气的平均含湿量dm由式4可得:
ε=(hm-h1)/(dm-d1)=(hw-h1)/(dw-d1) (4)
算出蒸发式冷凝器换热管外空气的平均含湿量dm,由空气焓湿图可查出换热管外空气的平均温度tm;
S5:计算水膜与空气的对流传热系数,初算换热面积;
S5-1:由输入风量G和初定尺寸长L和宽W,计算迎面风速v:
v=G/3600/((L-0.08)·(W-0.08))
计算最窄面处空气流速:
vmax=s/(s-do)·v (5)
式中,s为换热管与换热管之间的间距,d0为换热管外径,v为迎面风速。
S5-2:计算水膜与空气的对流传热系数αwa:
式中,λm为空气平均导热系数,υm为空气平均运动粘度;
S5-3:计算换热管外空气当量对流换热系数αj:
式中,A为水膜温度修正系数,取值范围为:0.94~0.99,hw为水膜处空气的焓值,hm为换热管外空气的平均焓值,Aw为水膜和空气的接触面积,Ao为换热管外表面积,Cpm为换热管外空气的平均定压比热,tw为水膜温度,tm为换热管外空气的平均温度;
式中,αj为管外空气当量对流换热系数,tw为水膜温度,tm为换热管外空气的平均温度;
S5-5:初算换热面积A’o:
S6:计算管外与喷水对流传热系数和计算制冷剂在管内对流传热系数,并根据使用情况估算热阻;
S6-1:管外与喷水对流传热系数αw:
式中,tw为水膜温度,Γ为喷淋密度,do为换热管外径;
式中,Mw为蒸发式冷凝器水泵流量。
S6-2:制冷剂在管内对流换热系数αc·n:
式中,λ为制冷剂导热系数,g为重力加速度,r为制冷剂在该压力下的汽化潜热,μ为制冷剂在该压力下液态时的动力粘度。
在蛇形管中换热修正:
S6-3:热阻
管壁热阻:铜管忽略不计,钢管等需考虑:
Rp=δ/λ (15)
式中,Rp为管壁热阻,δ为换热管壁厚,λ为管道导热系数。
油膜热阻:Roil,对于制冷剂氨取值0.35×10-3~0.6×10-3m2·K/W,对于氟利昂可忽略不计。
污垢热阻:Rfou,空气侧可取0.1×10-3~0.3×10-3m2·K/W。
S7:计算总传热系数和传热面积;
计算总换热系数K
式中,Ao为换热管外表面积,Ai为换热管内表面积,A平为换热管内外表面积平均值,Rp为管壁热阻,Roil为油膜热阻,Rfou为污垢热阻,αc·n·s为修正后的管内对流换热系数,αw为换热管外与喷水对流传热系数,αj为换热管外空气当量对流换热系数;
计算传热面积A”o:
A"o=Qc/K/(tk-tm) (17)
式中,Qc为机组热负荷,K为总换热系数,tk为制冷机组冷凝温度,tm为换热管外空气的平均温度。
S8:比较S7中计算得到的传热面积与S5中得到的初算换热面积是否相等,若误差小于1%,则输出需要的管排数,若误差大于等于1%,则重新假设水膜温度,进行迭代计算;
判断:
若|A"o-A′o|<1%,则进行以下步骤;若不是,则回到S3修改水膜温度,重新计算,直到满足条件。
管排数n
式中,n为管排数;A”o为传热面积;
一排传热面积=一排中的管数×管长×管道外径周长。
S9:假如计算的管排数为奇数排,修改管长使管排数在16~24范围内;
S10:输出最终的蒸发式冷凝器结果。
以下实施例按照上述步骤进行,根据计算结果,可以调整初定结构参数进行优化设计。
实施例1:
本实施例的蒸发式冷凝器要求冷凝负荷30kw,冷凝温度37℃,进口干球温度31℃,进口含湿量14.37g/kg。初定结构尺寸16mm管径铝管,1.5mm壁厚,管间距38mm。冷凝风量为2400m3/h,喷淋水流量1.36kg/s。计算得需要换热面积10.2m2,总传热系数为 624W/m2K,需要管排数28排,计算结构的面积余量为3.1588%。
实施例2:
本实施例蒸发式冷凝器要求冷凝负荷280kw,冷凝温度36℃,进口干球温度31℃,进口含湿量14.37g/kg。初定结构尺寸16mm管径铝管,1.5mm壁厚,管间距32mm。冷凝风量为45000m3/h,喷淋水流量13.89kg/s,管长600mm。计算提示“管排数计算为奇数排,请修改管长”,将管长调整到680mm,计算得需要换热面积54.6m2,总传热系数为1026.7W/m2K,需要管排数20排,计算结构的面积余量为2.935%。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (7)
1.一种蒸发式冷凝器结构设计方法,所述蒸发式冷凝器包括出风口、进风口和换热管,其特征在于:包括以下步骤:
S1:给定蒸发式冷凝器需要的冷却负荷、冷凝温度和所处环境的进风口干湿球温度;
S2:初定结构尺寸,包括换热管内外径、换热管管壁厚度、换热管与换热管之间的间距,输入风量和喷水流量;
S3:设定初始水膜温度;
S4:计算空气、水各点的状态参数,所述状态参数包括空气焓值、空气含湿量和温度;
S5:计算水膜与空气的对流传热系数,初算换热面积;
S6:计算管外与喷水对流传热系数和计算制冷剂在管内对流传热系数,并根据使用情况估算热阻;
S7:计算总传热系数和传热面积;
S8:比较S7中计算得到的传热面积与S5中的初算面积是否相等,若误差小于设定值,则输出需要的管排数,若误差大于等于设定值,则重新假设水膜温度,进行迭代计算;
S9:输出最终的蒸发式冷凝器结果;
所述S4中的空气各点的状态参数包括:进风口空气的焓值、进风口空气的含湿量、出风口空气的焓值、出风口空气的含湿量、换热管外空气的平均焓值、换热管外空气的平均含湿量;
所述水各点的状态参数包括:水膜处空气的焓值和水膜处空气的含湿量;
所述S5具体包括以下步骤:
S5-1:根据输入风量和初定结构尺寸,计算迎面风速v:
v=G/3600/((L-0.08)·(W-0.08))
式中,G为输入风量,L为长,W为宽;
根据式5计算最窄面处空气流速:
vmax=s/(s-do)·v (5)
式中,s为换热管与换热管之间的间距;do为换热管外径;v为迎面风速;
S5-2:计算水膜与空气的对流传热系数αwa:
式中;λm为空气平均导热系数,do为换热管外径,vmax为最窄面处空气流速,υm为空气平均运动粘度;
S5-3:计算换热管外空气当量对流换热系数αj:
式中,A为水膜温度修正系数,αwa为水膜与空气的对流传热系数,hw为水膜处空气的焓值,hm为换热管外空气的平均焓值,Aw为水膜和空气的接触面积,Ao为换热管外表面积,Cpm为换热管外空气的平均定压比热,tw为水膜温度,tm为换热管外空气的平均温度;
S5-5:初算换热面积A’o:
式中,Qc为机组热负荷。
2.根据权利要求1所述的一种蒸发式冷凝器结构设计方法,其特征在于:所述S6具体包括以下步骤:
S6-1:计算管外与喷水对流传热系数αw:
式中,tw为水膜温度,Γ为喷淋密度,do为换热管外径;
式中,Mw为蒸发式冷凝器水泵流量;
S6-2:计算制冷剂在管道中的对流换热系数αc·n:
式中,λ为制冷剂导热系数,g为重力加速度,r为制冷剂在该压力下的汽化潜热,μ为制冷剂在该压力下液态时的动力粘度;
S6-3:计算热阻,该热阻包括管壁热阻;
所述管壁热阻Rp根据式15可得:
Rp=δ/λ’ (15)
式中,δ为换热管管壁厚度,λ’为管道导热系数。
3.根据权利要求2所述的一种蒸发式冷凝器结构设计方法,其特征在于:所述热阻还包括油膜热阻和污垢热阻。
7.根据权利要求6所述的一种蒸发式冷凝器结构设计方法,其特征在于:所述S8中计算得到的管排数为奇数排,修改管长使管排数在16~24范围内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811622756.4A CN109614746B (zh) | 2018-12-28 | 2018-12-28 | 一种蒸发式冷凝器结构设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811622756.4A CN109614746B (zh) | 2018-12-28 | 2018-12-28 | 一种蒸发式冷凝器结构设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109614746A CN109614746A (zh) | 2019-04-12 |
CN109614746B true CN109614746B (zh) | 2021-12-03 |
Family
ID=66010852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811622756.4A Active CN109614746B (zh) | 2018-12-28 | 2018-12-28 | 一种蒸发式冷凝器结构设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109614746B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111278255B (zh) * | 2019-12-31 | 2021-09-21 | 南京理工大学 | 基于凝结传热的相变蓄热装置及其关键参数确定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7065970B2 (en) * | 2003-11-07 | 2006-06-27 | Harpster Joseph W C | Condensers and their monitoring |
CN101414321B (zh) * | 2008-11-14 | 2010-11-10 | 西安交通大学 | 化工用蒸发式冷却/凝器设计方法 |
CN103542621B (zh) * | 2013-09-26 | 2016-01-20 | 西安交通大学 | 一种通用型组合管径空调换热设备流路的设计方法 |
CN107992662B (zh) * | 2017-11-27 | 2021-05-18 | 郑州云海信息技术有限公司 | 重力热管空调系统优化设计及变工况参数的反向计算方法 |
-
2018
- 2018-12-28 CN CN201811622756.4A patent/CN109614746B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109614746A (zh) | 2019-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yi et al. | Visualization study of the influence of non-condensable gas on steam condensation heat transfer | |
He et al. | Comparative study on the performance of natural draft dry, pre-cooled and wet cooling towers | |
Shin et al. | An experimental study of condensation heat transfer inside a mini-channel with a new measurement technique | |
US20090030554A1 (en) | Cooling control device and method | |
Kuang et al. | Numerical simulation of frosting behavior and its effect on a direct-contact ambient air vaporizer | |
Hu et al. | Theoretical investigation of stable dropwise condensation heat transfer on a horizontal tube | |
Wen et al. | Development and experimental study of a novel plate dehumidifier made of anodized aluminum | |
CN109614746B (zh) | 一种蒸发式冷凝器结构设计方法 | |
Nouri-Borujerdi et al. | Experimental study of forced convection heat transfer from a cam shaped tube in cross flows | |
Min et al. | Theoretical analysis of water film evaporation characteristics on an adiabatic solid wall | |
Owen et al. | A numerical investigation of vapor flow in large air-cooled condensers | |
Raiyani et al. | Thermal analysis of hot wall condenser for domestic refrigerator | |
Olcay et al. | Experimental investigation of frost issue on various evaporators having different fin types | |
Tang et al. | Study on the heat and mass transfer performance of a tubular still enhanced by hydrophilic surface modification | |
Liu et al. | Study on surface condensate water removal and heat transfer performance of a minichannel heat exchanger | |
Peyghambarzadeh et al. | Statistical analysis of calcium sulfate scaling under boiling heat transfer | |
CN107729600B (zh) | 蒸发器仿真计算方法 | |
Chang et al. | Modeling and performance simulation of a gas cooler for a CO2 heat pump system | |
Hirbodi et al. | Experimental investigation of natural dehumidification over an annular finned-tube | |
du Plessis et al. | A validated discretized thermal model for application in bare tube evaporative coolers and condensers | |
CN111428184A (zh) | 一种板翅式换热器芯体尺寸计算方法 | |
Liu et al. | Theoretical investigation and experimental verification of a mathematical model for counter-flow spray separation tower | |
Timmermann et al. | Thermal performance of peripheral-finned tube evaporators under frosting | |
CN109959283A (zh) | 高温冷却器热力校核方法及系统 | |
Abraham et al. | Experimental studies on thermal spray-coated horizontal tubes for falling film evaporation in multi-effect desalination system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 211500 No. 8, Longzhong West Road, Longchi street, Liuhe Economic Development Zone, Nanjing, Jiangsu Applicant after: Nanjing Jiuding Environmental Technology Co.,Ltd. Address before: 211500 No. 8, Longzhong West Road, Longchi street, Liuhe Economic Development Zone, Nanjing, Jiangsu Applicant before: NANJING JIUDING REFRIGERATION & AIR-CONDITIONING EQUIPMENT CO.,LTD. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |