专利名称：：化工用蒸发式冷却/凝器设计方法
技术领域：
：本发明属于制冷空调
技术领域：
，具体涉及一种采用沿程分段处理的化工用蒸发式冷却/凝器设计方法。
背景技术：
：传统的换热器设计历史悠久，设计步骤趋于成熟，相关技术人员己经积累了丰富的设计经验。以管壳式换热器为例，其设计过程己经得至鹏范化，设计步骤发展非常完善，国内已经出现了针对该种换热器的计算机设计软件并已申请了发明专利(200410066599.5)。蒸发式冷却/凝器作为一种新型的换热设备，具备节水、节能、结构简单紧凑、安装维修方便等优点，与传统设备相比，可使换热效率进一步提高10%以上。我国在能源利用方面存在严重的利用率低，消耗量过大的问题，尤其体现在占总能耗七成的工业生产部门中。化学工业是我国重要的原材料工业和支柱产业，与钢铁、建材形成高能耗工业前三位，其能耗占全国工业总會巨耗的15%以上。可以说，要提高我国能源利用率，在很大程度上取决于化工生产中节能技术的推广和使用。蒸发式冷却/凝技术是一种利用汽化潜热的高效换热手段，在化工生产领域的应用可以有效降低成本、提高产品竞争力。然而，蒸发式冷却/凝器的设计方法目前发展还并不成熟，国内厂家大多依靠经验或参照传统的设计方法进行设计，存在设计效率低下，计算精度不足等缺点。而将其引用于化工领域又会遇到如下困难被冷却/凝对象为混合气体，各组分随温度降低各自发生凝结，换热系数难以确定，进而使得换热面积不易确定，导致计算得到的结构设计参数可靠度降低。
发明内容本发明的目的在于提供一种采用分段式处理和迭代计算方式化工用蒸发式冷却/凝器设计方法。为达到上述目的，本发明采用的技术方案是1)根据化工用蒸发式冷却/凝器的原料气成分建立可扩充物性数据库，可扩充物性数据库涵盖化工原料气成分的物性参数；2)根据换热管材料建立可扩充热物性数据库，该可扩充热物性数据库涵盖工程换热管材的传热参数；3)根据热力工况输入参数以及结构设计输入参数分别建立工况参数输入与结构参数输入的图形化模块；根据热力计算输出参数与结构设计输出参数分别建立热力计算输出模块与结构设计输出模块；4)根据热力工况输入参数以及混合气体冷却/凝过程，计算相应的总换热量以及沿程换热量分布；热力计算过程所需计算参数如下工质进口鹏7i，(K)工质出口驢『2，(K)进口压力A，(MPa)出口压力jv(MPa)体积流量a，(Nm3/h)第/成分的体积百分数伊:，(％)(1)混合工质进出口相态判定具体步骤如下(a)根据进口压力a及第/成分的体积百分数^计算第/成分的分压力(b)根据进出口温度石、r2分别采用线性插值得到该成分对应的饱和压力&、A;(c)对第/成分气体进行判定若乂一A〉A、则该成分发生冷却，进出口工质均为气体；若式一PX2，则该成分发生部分冷凝，进口工质为气体、出口工质为过冷液体与饱和气体；若A、A〉乂2，则降温过程中该成分在冷凝器进口处发生部分冷凝，部^体凝结为液体，剩余气体变为饱和气体，在冷凝器内部继续发生相变，出口工质为过冷液体与饱和气体；(2)第/成分气体的换热量由下式确定<formula>formulaseeoriginaldocumentpage13</formula>(l)其中，2'为该成分对应的换热量，单位为kW;2二为该成分的进口气体摩尔流量，单位为kmol/h;M为该成分的分子量，单位为kg/kmd;"与《分别为该成分的进出口气体比焓，单位为kJ/kg，总的换热量为各成分产生的换热量累加-<formula>formulaseeoriginaldocumentpage13</formula>为混合工质物质种类总数(2)(3)为获得管热管沿程热量分布并以itkiS行结构设计，按温度分段热力计算，即按照进出口温度求得温差Ar:<formula>formulaseeoriginaldocumentpage13</formula>(3)然后将温差进行"等分，贝晦一个驢段内的进出口鹏为<formula>formulaseeoriginaldocumentpage13</formula>(4)<formula>formulaseeoriginaldocumentpage13</formula>(5)针对每一个温度段，返回热力计算步骤(1)~(2)，此处每一温度段内的出口参数，包括压力、流量，也即下一温度段的入口参数，将各个温度段内计算得到的分段换热量^进行累加得到全段总热量G，艮P..g-^G、"为按温度分段的总段数(6)5)热力计算模±央针对每一种混合气体成分提供了自动与手动两种相变判定方式；6)将设计所需参数，包括环境参数、结构参数、其它参数，结合由热力计算模±央得到的沿程换热量分布一同传递至结构设计输入模块；7)结构设计模块根据输入参数采用迭代计算，内迭代求解换热系数，外迭代在满足换热要求前提下保证压降满足压力损失要求，对混合气体发生部分相变条件下的换热系数求解，即对应内迭代；根据换热量即热力计算提供及换热系数，求解换热面积及相应结构，若解满足换热要求不满足压力要求，调整参数，迭代计算同时满足换热与压力要求的解，即外迭代，热力计算与结构设计模块的计算结果即混合气体以及各组分气体的总换热量及沿程换热量分布，换热管排数与每排根数，配风量、配水量以及换热器外形尺寸，即长、宽、高传送至输出模块；结构设计模块的计算过程如下所需计算参数如下工质进口温度《，(K)工质出口温度r2，(K)操作压力P，(Pa)体积流量-a，第/成分的体积百分数W，(%)空气进口比烚/，(kJ/kg)空气出口比焓'"(kJ/kg)横断面管长/。，(m)管子内径《，(m)管子外径《，(m)迎面风速wf，(m/s)空气进口温度7;°，(K)空气出口温度7f，(饱和水、温度)(K)喷淋水温度"2。，(K)配风系数sc，(m3/MJ)配水系数fH2°，(m3/MJ)计算中使用到的其他参数定义见下:<table>tableseeoriginaldocumentpage15</column></row><table><table>tableseeoriginaldocumentpage16</column></row><table>(1)假设喷淋水温度"2°=303&，并给定初始配风系数、配水系数sG=15m3/MJ，sH2O=0.012m3/MJ;(2)采用如下公式计算总换热系数<formula>formulaseeoriginaldocumentpage16</formula>(7)式(7)中，A为冷却盘管内混合工质换热系数，其计算方法如下对于冷凝管内混合工质中冷凝成分，将其定压降温过程分为过热气冷却至饱和气(g-g)、饱和气冷凝为饱和液(g-l)、饱和液冷却至过冷液(1-1)三个阶段分别进行换热系数计算；(a)过热气体冷却对流换热系数、采用格林尼斯基公式<formula>formulaseeoriginaldocumentpage16</formula>(8)对于气体有:<formula>formulaseeoriginaldocumentpage17</formula>以及公式-<formula>formulaseeoriginaldocumentpage17</formula>(9)计算得到，其中管内流体定性温度ff，i-，G,、^分别为冷却器进出口处温度，°C);管内湍流流动的达尔西阻力系数/=(1.8218勘-1.64)—2;(b)相变换热系数&采用分段热力计算结果，在一定温度间隔内&采用下面公式计算:<formula>formulaseeoriginaldocumentpage17</formula>(10)其中，W表示不凝性气体的对流换热系数，其计算方法与上述(a)部分完全相同；Q(为第/温度段内的潜热换热量；(c)过冷液对流换热系数、其计算方法与上述(a)类似，区别仅在对于液体有:<formula>formulaseeoriginaldocumentpage17</formula>(d)根据换热面积关系可以确定管内混合工质平均对流换热系数^:<formula>formulaseeoriginaldocumentpage17</formula>(11)其中，S表示总换热量，0>表示(a)过程换热量，込表示(b)过程换热量，a表示(C)过程换热量；<=(7;+7^/2-"2°，A^表示(a)过程中工质与管外水膜的平均温差，A"表示(b)过程中工质与管外水膜的平均温差，，。表示(C)过程中工质与管外水膜的平均温差；式(7)中，iw为冷却盘管导热热阻，采用下式计算得到及w-^"2)式(7)中，A。为管外冷却水的对流换热系数，采用下式计算得到:其中，H,O4"/°式(7)中，《与i。分别为管内外壁污垢热阻，取i^0.002m2.K/W;(3)冷却水与空气间换热采用传质系数尺。^进行计算，艮卩Y1=3.62x104(i^广(&。广5(《广I，(14)o乂上式中各系数确定如下ief。=4m/Wf、i^=《Gm/M(J、Gffl=/n气、'。；并可以计算出以下各参数^"^-，5=^，"=S;(4)计算主要结构参数管排数AT，则需要对下述两式进行求解，得到Z,与4:屮,In1(15)(16)其中，^和4^是下述方程的两个根，f+a^+bJT+Ovbrb2.b3)=0其中，参数由式(7)与式(14)的计算结果以及如下各式确定b,=(ai+a3)，a,=~:——~;b7=-a，，a7="un;b,=aa,,a,=~^—~;b4=aa，_a4，a4=6「1g,,，2220H2OH20,333.C'，4244gGa=(5)判断Zi与Z,之差是否小于工程许可值，计算中取IZ，-Z^0.02lZ^^I作(13)/z00011I*II为收敛条件，如果满足收敛条件则用公式iV-2Z/^^得到管排数iV，其中Z"Z1+Z2)/2，进而计算得到外形几何参数；如果不满足收敛条件，重新假定结构设计部分步骤l中的7^g，按步骤(2)中(a)~(d)的顺序迭代求解各个变量，直到满足内迭代收敛条件为止；(6)计算冷却/凝盘管进出口压损，利用勃拉休斯公式计算沿程摩擦系数/-<formula>formulaseeoriginaldocumentpage19</formula>利用下面公式计算冷却/凝盘管进出口压损<formula>formulaseeoriginaldocumentpage19</formula>其中，《取值为0.27，如果压损Ap大于工质进出口允许压损，则增大开始输入的配风系数，对步骤(1)~(6)重复计算，直到满足外迭代的压损条件为止；8)输出模块将热力计算与结构设计模±央计算结果以Excel电子表格形式输出。本发明的可扩充物性数据库涵盖氩气、一氧化碳、二氧化碳、氢气、水蒸汽、氮气、氨气、氧气、甲醇、乙醇、二甲醚、空气、甲烷、乙烷、丙烷、丁烷、异丁烷、戊烷、己烷、庚烷、辛烷、壬烷、葵烷、乙烯、丙烯，数据内容包含气体名称、化学式、分子量、临界M、临界压力、三相点温度、饱和鹏、饱和压力、比烚、比定压热容、导热系数、粘度、表面张力，数据范围为温度-6.85°0226.85°(:、压力0MPa100MPa，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充；可扩充热物性数据库涵盖铝、不锈钢、黄铜、青铜、10#碳钢、20#碳钢的传热系数，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充；沿程换热量分布包括潜热换热量、显热换热量与总换热量，是将气体的冷却凝过程按进出口温差进行若干等份，对每一段内分别进行热力计算得到的；相变判定方式，在饱和点附近的工况通过手工判定得到的计算结果与自动判定计算结果对照的方式保证结果的可靠性；结构设计模块采用内外双层迭代计算方式求解，计算结果包括，热力参数压力损失、换热温差、换热管面积、总换热系数;结构参数换热管排数、每排换热管数；其它参数配风量、酉己水量。本发明使用方便、集成化程度高、计算精度高、可扩充性强，使有关设计人员迅速高效地完成蒸发式冷却/凝器的设计工作，从而大大降低企业的开发成本与周期、提高设计效率。具体实施例方式1)根据化工用蒸发式冷却/凝器的原料气成分建立可扩充物性数据库，可扩充物性数据库涵盖氩气、一氧化碳、二氧化碳、氢气、水蒸汽、氮气、氨气、氧气、甲醇、乙醇、二甲醚、空气、甲烷、乙烷、丙烷、丁烷、异丁烷、戊烷、己烷、庚烷、辛烷、壬烷、葵烷、乙烯、丙烯，数据内容包含气体名称、化学式、分子量、临界温度、临界压力、三相点温度、饱和温度、饱和压力、比焓、比定压热容、导热系数、粘度、表面张力，数据范围为温度-6.85i:226.85t:、压力0MPa100MPa，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充；2)根据换热管材料建立可扩充热物性数据库，该可扩充热物性数据库涵盖铝、不锈钢、黄铜、青铜、画碳钢、20#碳钢的传热系数，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充；3)根据热力工况输入参数以及结构设计输入参数分别建立工况参数输入与结构参数输入的图形化模块；根据热力计算输出参数与结构设计输出参数分别建立热力计算输出模块与结构设计输出模块；4)根据热力工况输入参数以及混合气体冷却/凝过程，计算相应的总换热量以及沿程换热量分布，该沿短换热量分布包括潜热换热量、显热换热量与总换热量，是将气体的冷却/凝过程按进出口温差进行若干等份，对每一段内分别进行热力计算得到的；热力计算过程所需计算参数如下工质进口温度7;，(K)工质出口温度r2，(K)进口压力A，(MPa)出口压力/72，(MPa)体积流量gv，(Nm3/h)第/成分的体积百分数^，(％)(1)混合工质进出口相态判定具体步骤如下-(a)根据进口压力A及第/成分的体积百分数^计算第z'成分的分压力(b)根据进出口温度;、r2分别采用线性插值得到该成分对应的饱和压力&、(c)对第Z成分气体进行判定若A一/4〉A、则该成分发生冷却，进出口工质均为气体；若K夕/^X，则该成分发生部分冷凝，进口工质为气体、出口工质为过冷液体与饱和气体；若PX"^则降温过程中该成分在冷凝器进口处发生部分冷凝，部分气体凝结为液体，剩余气体变为饱和气体，在冷凝器内部继续发生相变，出口工质为过冷液体与饱和气体；(2)第Z'成分气体的换热量由下式确定-g《,M("-g)/3600(1)其中，"为该成分对应的换热量，单位为kW;2"为该成分的进口气体摩尔流量，单位为kmol/h;M为该成分的分子量，单位为kg/kmol;"与g分别为该成分的进出口气体比烚，单位为kJ/kg。总的换热量为各成分产生的换热量累加e=|;《，附为混合工质物质种类总数(2)1(3)为获得管热管沿程热量分布并以此进行结构设计，按温度分段热力计算，即按照进出口温度求得温差Ar:△r=7;-r2(3)然后将温差进行"等分，贝梅一个温度段内的进出口温度为.-r/=7;—(j.—l)(A77w)，J=l，2，3.."7_l(4)r/=r/-A77w,风2，3…，w-1(5)针对每一个温度段，返回热力计算步骤(1)~(2)，此处每一温度段内的出口参数，包括压力、流量，也即下一温度段的入口参数。将各个温度段内计算得到的分段换热量^进行累加得到全段总热量e，艮P-0=t^，"为按温度分段的总段数(6)15)热力计算模i央针对每一种混合气体成分提供了自动与手动两种相变判定方式；该相变判定方式在饱和点附近的工况通过手工判定得到的计算结果与自动判定计算结果对照的方式保证结果的可靠性；6)将设计所需参数，包括环境参数、结构参数、其它参数，结合由热力计算模块得到的沿程换热量分布一同传递至结构设计输入模i央；7)结构设计模块根据输入参i^用迭代计算，该结构设计模±央采用内外双层迭代计算方式求解，计算结果包括，热力参数压力损失、换热温差、换热管面积、总换热系数；结构参数换热管排数、每排换热管数；其它参数配风量、配水量；内迭代求解换热系数，外迭代在满足换热要求前提下保证压降满足压力损失要求，对混合气体发生部分相变条件下的换热系数求解，即对应内迭代；根据换热量即热力计算提供及换热系数，求解换热面积及相应结构，若解满足换热要求不满足压力要求，调整参数，迭代计算同时满足换热与压力要求的解，即外迭代，热力计算与结构设计模块的计算结果即混合气体以及各组分气体的总换热量及沿程换热量分布，换热管排数与每排根数，配风量、配水量以及换热器外形尺寸，即长、宽、高传送至输出模块；结构设计计算过程所需计算参数如下工质进口温度工质出口温度操作压力体积流量第/成分的体积百分数:空气进口比烚-空气出口比烚横断面管长管子内径管子外径迎面风速空气进PI鹏空气出口鹏喷淋水温度配风系数配水系数:t;，(k)r2，P，(Pa)Gr，(Nm3/h)W，(o/o)/，(kJ/kg)z'2，(kJ/kg)/。，(m)t/j，(m)《，(m)Wf，(m/s)t;g，(k)r2G，(饱称Jc鹏)(K)"20，(K)eG，(m3/MJ)C(m3/MJ)计算中使用到的其他参数定义见下:名称单位<table>tableseeoriginaldocumentpage24</column></row><table>P、空气密度kg/m3(1)假设喷淋水温度(如严。=30310,并给定初始配风系数、配水系数(如sG=15m3/MJ，sH2O=0.012m3/MJ)。(2)采用如下公式计算总换热系数《=-1+《)&+及0+及4Q、(7)义+-式(7)中，A为冷却盘管内混合工质换热系数，其计算方法如下对于冷凝管内混合工质中冷凝成分，将其定压降温过程分为过热气冷却至饱和气(g-g)、饱和气冷凝为饱和液(g-l)、饱和液冷却至过冷液(1-1)三个阶段分别进行换热系数计算；(a)过热气体冷却对流换热系数A,,采用格林尼斯基公式(//8)(ie-1000)/Vf1+12.7V77^(/Vf2/3-l)1+《、2/3对于气体有:/、0.45)、r、乙以及公式:^如(9)计算得到，其中管内流体定性温度^=^(f,、^分别为冷却器进出口处温度，°C);管内湍流流动的达尔西阻力系数/=(1.82^^-1.64)一2。(b)相变换热系数&采用分段热力计算结果，在一定温度间隔内、采用下面公式计算:e乂(10)其中，W表示不凝性气体的对流换热系数，其计算方法与上述(a)部分完全相同；^为第J'温度段内的潜热换热量。(c)过冷液对流换热系数\。其计算方法与上述(a)类似，区别仅在对于液体有、ou、0,05~20(d)根据换热面积关系可以确定管内混合工质平均对流换热系数A,23*m,S，a^""m，a"i,b^""m，b"i，c盯(11)其中，^表示总换热量，辽表示(a)过程换热量，0>表示(b)过程换热量，a表示(c)过程换热量；=(r1+r2)/2-rH2。，a^表示(a)过程中工质与管外水膜的平均温差，aq表示(b)过程中工质与管外水膜的平均温差，a"表示(c)过程中工质与管外水膜的平均温差。式(7)中，&为冷却盘管导热热阻，采用下式计算得到及w=:(12)式(7)中，h。为管外冷却水的对流换热系数，采用下式计算得到:^=118/、i附、《乂(13)其中，'H204"Z式(7)中，《与及。分别为管内外壁污垢热阻，取《-0.002m2.K/W(3)冷却水与空气间换热采用传质系数iC。^进行计算，即^g"-3.62x104(仏g)。9(年。)。15(《)_261A义(14)上式中各系数确定如下^f。=4w/"f、ieG=《Gm/W。、Gw=/ng;、,。；并可以计算出以下各参数&，，S=^，<formula>formulaseeoriginaldocumentpage27</formula>(4)计算主要结构参数管排数iV，则需要对下述两式进行求解，得到Z,与Z,1<formula>formulaseeoriginaldocumentpage27</formula>(15)<formula>formulaseeoriginaldocumentpage27</formula>(16)其中，^和^是下述方程的两个根，甲2+(b,+b4)甲+(Vb4—b2.b3)=0其中，参数由式(7)与式(14)的计算结果以及如下各式确定b广(a,+^)，a,=~^~t";b，=-a，，a，=_n;b,=aa,,a,=~^;b4=aa2—a4，a4=.s—a=-《.kr加cp0'c'(5)判断Z,与^之差是否小于工程许可值，计算中取IZ广Z^0.02lZ,+Z2l作为收敛条件。如果满足收敛条件则用公式iV-2Z/^^得到管排数iV，其中Z=(Z1+Z2)/2，进而计算得到外形几何参数；如果不满足收敛条件，重新假定结构设计部分步骤l中的rw，按步骤(2)中(a)~(d)的顺序迭代求解各个变量，直到满足收敛条件为止(内迭代)。(6)计算冷却凝盘管进出口压损，利用勃拉休斯公式计算沿程摩擦系数/:<formula>formulaseeoriginaldocumentpage27</formula>(17)禾U用下面公式计算冷却/凝盘管进出口压损:<formula>formulaseeoriginaldocumentpage27</formula>(18)其中，《取值为0.27。如果压损Ap大于工质进出口允许压损，则增大开始输入的配风系数，对步骤(1)~(6)重复计算，直到满足压损条件为止(外迭代)。278)输出模±央将热力计算与结构设计模±央计算结果以Excel电子表格形式输出。下面结合具有代表性的甲醇合成段冷凝器设计工况，对本发明的化工用蒸发式冷却/凝器计算机设计方、法作进一步说明。热力计算工况参数为气体成分及体积含量氢气73.89%、氮气9.34%、甲烷2.6%、一氧化碳5.4%、二氧化碳3.64%、水蒸汽1.13%、甲醇3.99%;进口参数温度100。C、压力5.65MPa;出口参数温度40。C、压力5.60MPa;允许压降0.05MPa;体积流量28670Nm3/h。结构设计要求参数为，环境参数海拔高度50m、干球，3(TC、湿球，25°C;结构参数(差排布置)纵向间距75mm、横向间距75mm、管长3m、管外径30mm、管内径26mm、管壁热阻0.0002m2K/W、换热管材料为不锈钢；其它参数配风系数(初值)15m3/MJ、配水系数(初值)0.012m3/MJ、迎风风速3m/s。具体设计步骤如下(1)根据25种化工领域常见的原料气成分建立可扩充物性数据库，内容包括气体名称、化学式、分子量、临界温度、临界压力、三相点温度、饱和温度、饱和压力、比烚、比定压热容、导热系数、粘度、表面张力等。温度与压力的覆盖范围分别为-6.85。C226.85T:与0MPa100MPa。(2)根据换热管材料建立可扩充热物性数据库，内容主要为换热管材料名称与导热系数等。(3)根据工况输入参数以及结构输入参数在输入界面进行输入操作，程序进行相应的参数传递。(4)根据热力计算所需参数，在数据库中查找对应比烚值计算得到总换热量为1241,348kW，其中总显热换热量为667.408kW，总潜热换热量为573.940kW。沿程换热量计算采用将气体的冷却/凝过程按进出口温差进行20等份(可在1~20范围内输入任意整数值)，对每一段内分别进行热力计算得到。由于数据量过大，此处仅列出混合气体计算结果，如下表所示:<table>tableseeoriginaldocumentpage29</column></row><table>(5)热力计算模±央针对每一种混合气体成分提供了自动与手动两种方式的相变判定方式，可有效避免状态点与饱和点过于接近而产生的误差。此处采用自动判断方式进行，结果显示水蒸汽与甲醇两种成分在降温过程中发生冷凝，冷凝开始的M段为88-C85X:。(6)将结构设计要求参数，包括环境参数、结构参数以及其它参数，结合由热力计算模i央得到的沿程换热量分布一同传递至结构设计输入模块中。(7)结构设计模块根据输入参数对问皿tm代计算求解，以保证换热能力满足要求的同时压力损失也低于允许压降。计算过程中得到的从管内混合气体至管外水膜的沿程换热系数分布如下鹏段('c)从管内混^H体至管外水膜的沿程换热系数(W/m2K)<table>tableseeoriginaldocumentpage30</column></row><table>同时考虑管外侧水膜与空气间的传热传质关系并进行迭代计算，求得最终的总换热系数96.807W/m2K。其他主要计算结果包括，热力参数:压力损失0.009848MPa、换热温差34.000057-C、换热管面积379.918m2;结构参数换热管排数40、每排换热管数34;其它参数配风量73806.547m3/h、配水量58.7831^/h等。由上述参数设计得到的蒸发式冷凝器外观尺寸约为长乂宽><高=3.00111><2.60111><2.57111。(8)热力计算与结构设计模块均制定有完整详细的输入输出参数列表，并采用Excel电子表格形式呈现，方便对数据的进一步分析处理。权利要求1、一种化工用蒸发式冷却/凝器的设计方法，其特征在于包括如下步骤1)根据化工用蒸发式冷却/凝器的原料气成分建立可扩充物性数据库，可扩充物性数据库涵盖化工原料气成分的物性参数；2)根据换热管材料建立可扩充热物性数据库，该可扩充热物性数据库涵盖工程换热管材的传热参数；3)根据热力工况输入参数以及结构设计输入参数分别建立工况参数输入与结构参数输入的图形化模块；根据热力计算输出参数与结构设计输出参数分别建立热力计算输出模块与结构设计输出模块；4)根据热力工况输入参数以及混合气体冷却/凝过程，计算相应的总换热量以及沿程换热量分布；热力计算过程所需计算参数如下工质进口温度T1，(K)工质出口温度T2，(K)进口压力p1，(MPa)出口压力p2，(MPa)体积流量Qv，(Nm3/h)第i成分的体积百分数(1)混合工质进出口相态判定具体步骤如下(a)根据进口压力p1及第i成分的体积百分数计算第i成分的分压力(b)根据进出口温度T1、T2分别采用线性插值得到该成分对应的饱和压力(c)对第i成分气体进行判定若则该成分发生冷却，进出口工质均为气体；若则该成分发生部分冷凝，进口工质为气体、出口工质为过冷液体与饱和气体；若则降温过程中该成分在冷凝器进口处发生部分冷凝，部分气体凝结为液体，剩余气体变为饱和气体，在冷凝器内部继续发生相变，出口工质为过冷液体与饱和气体；(2)第i成分气体的换热量由下式确定其中，Qi为该成分对应的换热量，单位为kW；为该成分的进口气体摩尔流量，单位为kmol/h；M为该成分的分子量，单位为kg/kmol；与分别为该成分的进出口气体比焓，单位为kJ/kg，总的换热量为各成分产生的换热量累加m为混合工质物质种类总数(2)(3)为获得管热管沿程热量分布并以此进行结构设计，按温度分段热力计算，即按照进出口温度求得温差ΔTΔT＝T1-T2(3)然后将温差进行n等分，则每一个温度段内的进出口温度为针对每一个温度段，返回热力计算步骤(1)～(2)，此处每一温度段内的出口参数，包括压力、流量，也即下一温度段的入口参数，将各个温度段内计算得到的分段换热量Qj进行累加得到全段总热量Q，即n为按温度分段的总段数(6)5)热力计算模块针对每一种混合气体成分提供了自动与手动两种相变判定方式；6)将设计所需参数，包括环境参数、结构参数、其它参数，结合由热力计算模块得到的沿程换热量分布一同传递至结构设计输入模块；7)结构设计模块根据输入参数采用迭代计算，内迭代求解换热系数，外迭代在满足换热要求前提下保证压降满足压力损失要求，对混合气体发生部分相变条件下的换热系数求解，即对应内迭代；根据换热量即热力计算提供及换热系数，求解换热面积及相应结构，若解满足换热要求不满足压力要求，调整参数，迭代计算同时满足换热与压力要求的解，即外迭代，热力计算与结构设计模块的计算结果即混合气体以及各组分气体的总换热量及沿程换热量分布，换热管排数与每排根数，配风量、配水量以及换热器外形尺寸，即长、宽、高传送至输出模块；结构设计模块的计算过程如下所需计算参数如下工质进口温度T1，(K)工质出口温度T2，(K)操作压力p，(Pa)体积流量QV，(Nm3/h)第i成分的体积百分数空气进口比焓i1，(kJ/kg)空气出口比焓i2，(kJ/kg)横断面管长l°，(m)管子内径di，(m)管子外径do，(m)迎面风速wf，(m/s)空气进口温度空气出口温度(饱和水温度)(K)喷淋水温度配风系数εG，(m3/MJ)配水系数计算中使用到的其他参数定义见下(1)假设喷淋水温度并给定初始配风系数、配水系数εG＝15m3/MJ，(2)采用如下公式计算总换热系数式(7)中，hi为冷却盘管内混合工质换热系数，其计算方法如下对于冷凝管内混合工质中冷凝成分，将其定压降温过程分为过热气冷却至饱和气(g-g)、饱和气冷凝为饱和液(g-1)、饱和液冷却至过冷液(l-1)三个阶段分别进行换热系数计算；(a)过热气体冷却对流换热系数hi，a采用格林尼斯基公式对于气体有以及公式计算得到，其中管内流体定性温度(t1、t2分别为冷却器进出口处温度，℃)；管内湍流流动的达尔西阻力系数f＝(1.821gRe-1.64)-2；(b)相变换热系数hi，b采用分段热力计算结果，在一定温度间隔内hi，2采用下面公式计算其中，表示不凝性气体的对流换热系数，其计算方法与上述(a)部分完全相同；为第j温度段内的潜热换热量；(c)过冷液对流换热系数hi，c其计算方法与上述(a)类似，区别仅在对于液体有(d)根据换热面积关系可以确定管内混合工质平均对流换热系数hi其中，Q表示总换热量，Qa表示(a)过程换热量，Qb表示(b)过程换热量，Qc表示(c)过程换热量；Δtm，a表示(a)过程中工质与管外水膜的平均温差，Δtm，b表示(b)过程中工质与管外水膜的平均温差，Δtm，c表示(c)过程中工质与管外水膜的平均温差；式(7)中，Rw为冷却盘管导热热阻，采用下式计算得到式(7)中，ho为管外冷却水的对流换热系数，采用下式计算得到其中，式(7)中，Ri与Ro分别为管内外壁污垢热阻，取Ri＝0.002m2·K/W；(3)冷却水与空气间换热采用传质系数Koga进行计算，即上式中各系数确定如下Ref，o＝4m/uf、ReG＝doGm/uG、并可以计算出以下各参数(4)计算主要结构参数管排数N，则需要对下述两式进行求解，得到Z1与Z2其中，Ψ1和Ψ2是下述方程的两个根，Ψ2+(b1+b4)Ψ+(b1·b4-b2·b3)＝0其中，参数由式(7)与式(14)的计算结果以及如下各式确定b1＝(a1+a3)，b2＝-a2，b3＝αa3，b4＝αa2-a4，(5)判断Z1与Z2之差是否小于工程许可值，计算中取|Z1-Z2|≤0.02|Z1+Z2|作为收敛条件，如果满足收敛条件则用公式得到管排数N，其中Z＝(Z1+Z2)/2，进而计算得到外形几何参数；如果不满足收敛条件，重新假定结构设计部分步骤1中的按步骤(2)中(a)～(d)的顺序迭代求解各个变量，直到满足内迭代收敛条件为止；(6)计算冷却/凝盘管进出口压损，利用勃拉休斯公式计算沿程摩擦系数f利用下面公式计算冷却/凝盘管进出口压损其中，ξ取值为0.27，如果压损Δp大于工质进出口允许压损，则增大开始输入的配风系数，对步骤(1)～(6)重复计算，直到满足外迭代的压损条件为止；8)输出模块将热力计算与结构设计模块计算结果以Excel电子表格形式输出。2、根据权利要求1所述的蒸发式冷却/凝器的计算机设计方法，其特征在于所说的可扩充物性数据库涵盖氩气、一氧化碳、二氧化碳、氢气、水蒸汽、氮气、氨气、氧气、甲醇、乙醇、二甲醚、空气、甲垸、乙烷、丙垸、丁垸、异丁烷、戊烷、己烷、庚烷、辛烷、壬烷、葵烷、乙烯、丙烯，数据内容包含气体名称、化学式、分子量、临界温度、临界压力、三相点温度、饱和温度、饱和压力、比烚、比定压热容、导热系数、粘度、表面张力，数据范围为温度-6.85r226.85t:、压力0MPa100MPa，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充。3、根据权利要求1所述的蒸发式冷却/凝器的计算机设计方法，其特征在于所说的可扩充热物性数据库涵盖铝、不锈钢、黄铜、青铜、10#碳钢、20#碳钢的传热系数，且所有数据通过外部数据文件导入，通过外部数据文件进行数据扩充。4、根据权利要求1所述的蒸发式冷却/凝器的计算机设计方法，其特征在于所说的沿程换热量分布包括潜热换热量、显热换热量与总换热量，是将气体的冷却/凝过程按进出口温差进行若干等份，对每一段内分别进行热力计算得到的。5、根据权利要求1所述的蒸发式冷却/凝器的计算机设计方法，其特征在于所说的相变判定方式，在饱和点附近的工况通过手工判定得到的计算结果与自动判定计算结果对照的方式保证结果的可靠性。6、根据权利要求1所述的蒸发式冷却/凝器的计算机设计方法，其特征在于所说的结构设计模块采用内外双层迭代计算方式求解，计算结果包括，热力参数压力损失、换热温差、换热管面积、总换热系数；结构参数换热管排数、每排换热管数；其它参数配风量、配水量。全文摘要一种化工用蒸发式冷却/凝器的设计方法，根据化工流程中被冷却/凝混合气体计算沿程放热量，即热力计算部分，以及根据上述结果计算相应的冷却/凝器设计参数，即结构设计部分。针对化工领域常见的25种原料气成分建立可扩充物性数据库，针对换热管材料建立相应的可扩充数据库；将冷却/凝过程按进出口温差沿程进行若干等份，对每一段内分别进行热力计算以获得沿程数据分布，并以此结果进行结构设计；将得到的数据输入到输入模块、热力计算模块、结构设计模块由输出模块输出设计结果。本发明采用分段式处理和迭代计算思想，可有效提高计算精度；集成了蒸发式冷却/凝器的热力计算与结构设计过程，可有效降低企业开发周期与成本，提高设计效率。文档编号G06F17/50GK101414321SQ200810232269公开日2009年4月22日申请日期2008年11月14日优先权日2008年11月14日发明者何雅玲,冯秀荣,刘怀亮,杰孙,董晓强申请人:西安交通大学