[发明专利]化工用蒸发式冷却/凝器设计方法

摘要

一种化工用蒸发式冷却/凝器的设计方法，根据化工流程中被冷却/凝混合气体计算沿程放热量，即热力计算部分，以及根据上述结果计算相应的冷却/凝器设计参数，即结构设计部分。针对化工领域常见的25种原料气成分建立可扩充物性数据库，针对换热管材料建立相应的可扩充数据库；将冷却/凝过程按进出口温差沿程进行若干等份，对每一段内分别进行热力计算以获得沿程数据分布，并以此结果进行结构设计；将得到的数据输入到输入模块、热力计算模块、结构设计模块由输出模块输出设计结果。本发明采用分段式处理和迭代计算思想，可有效提高计算精度；集成了蒸发式冷却/凝器的热力计算与结构设计过程，可有效降低企业开发周期与成本，提高设计效率。

主权项

1、一种化工用蒸发式冷却/凝器的设计方法，其特征在于：包括如下步骤：1)根据化工用蒸发式冷却/凝器的原料气成分建立可扩充物性数据库，可扩充物性数据库涵盖化工原料气成分的物性参数；2)根据换热管材料建立可扩充热物性数据库，该可扩充热物性数据库涵盖工程换热管材的传热参数；3)根据热力工况输入参数以及结构设计输入参数分别建立工况参数输入与结构参数输入的图形化模块；根据热力计算输出参数与结构设计输出参数分别建立热力计算输出模块与结构设计输出模块；4)根据热力工况输入参数以及混合气体冷却/凝过程，计算相应的总换热量以及沿程换热量分布；热力计算过程：所需计算参数如下：工质进口温度：                        T₁，(K)工质出口温度：                        T₂，(K)进口压力：                            p₁，(MPa)出口压力：                            p₂，(MPa)体积流量：                            Q_v，(Nm³/h)第i成分的体积百分数：(1)混合工质进出口相态判定具体步骤如下：(a)根据进口压力p₁及第i成分的体积百分数计算第i成分的分压力(b)根据进出口温度T₁、T₂分别采用线性插值得到该成分对应的饱和压力(c)对第i成分气体进行判定：若ps,1i>ps,2i>p1i,]]>则该成分发生冷却，进出口工质均为气体；若ps,1i>p1i>ps,2i,]]>则该成分发生部分冷凝，进口工质为气体、出口工质为过冷液体与饱和气体；若p1i>ps,1i>ps,2i,]]>则降温过程中该成分在冷凝器进口处发生部分冷凝，部分气体凝结为液体，剩余气体变为饱和气体，在冷凝器内部继续发生相变，出口工质为过冷液体与饱和气体；(2)第i成分气体的换热量由下式确定：Qi=Qm,1iM(h1i-h2i)/3600---(1)]]>其中，Qⁱ为该成分对应的换热量，单位为kW；为该成分的进口气体摩尔流量，单位为kmol/h；M为该成分的分子量，单位为kg/kmol；与分别为该成分的进出口气体比焓，单位为kJ/kg，总的换热量为各成分产生的换热量累加：Q=Σ1mQi,]]>m为混合工质物质种类总数      (2)(3)为获得管热管沿程热量分布并以此进行结构设计，按温度分段热力计算，即按照进出口温度求得温差ΔT：ΔT＝T₁-T₂                            (3)然后将温差进行n等分，则每一个温度段内的进出口温度为：T1j=T1-(j-1)(ΔT/n),j=1,2,3...,n-1---(4)]]>T2j=T1j-ΔT/n,j=1,2,3...,n-1---(5)]]>针对每一个温度段，返回热力计算步骤(1)～(2)，此处每一温度段内的出口参数，包括压力、流量，也即下一温度段的入口参数，将各个温度段内计算得到的分段换热量Q^j进行累加得到全段总热量Q，即：Q=Σ1nQj,]]>n为按温度分段的总段数(6)5)热力计算模块针对每一种混合气体成分提供了自动与手动两种相变判定方式；6)将设计所需参数，包括环境参数、结构参数、其它参数，结合由热力计算模块得到的沿程换热量分布一同传递至结构设计输入模块；7)结构设计模块根据输入参数采用迭代计算，内迭代求解换热系数，外迭代在满足换热要求前提下保证压降满足压力损失要求，对混合气体发生部分相变条件下的换热系数求解，即对应内迭代；根据换热量即热力计算提供及换热系数，求解换热面积及相应结构，若解满足换热要求不满足压力要求，调整参数，迭代计算同时满足换热与压力要求的解，即外迭代，热力计算与结构设计模块的计算结果即混合气体以及各组分气体的总换热量及沿程换热量分布，换热管排数与每排根数，配风量、配水量以及换热器外形尺寸，即长、宽、高传送至输出模块；结构设计模块的计算过程如下：所需计算参数如下：工质进口温度：                      T₁，(K)工质出口温度：                      T₂，(K)操作压力：                          p，(Pa)体积流量：                          Q_V，(Nm³/h)第i成分的体积百分数：空气进口比焓：                      i₁，(kJ/kg)空气出口比焓：                      i₂，(kJ/kg)横断面管长：                         l°，(m)管子内径：                           d_i，(m)管子外径：                               d_o，(m)迎面风速：                               w_f，(m/s)空气进口温度：空气出口温度：(饱和水温度)(K)喷淋水温度：配风系数：                               ε^G，(m³/MJ)配水系数：计算中使用到的其他参数定义见下：(1)假设喷淋水温度TH2O=303K,]]>并给定初始配风系数、配水系数ε^G＝15m³/MJ，ϵH2O=0.012m3/MJ;]]>(2)采用如下公式计算总换热系数：K=1(1/hi+Ri)dodi+Ro+Rw(dodm)+1ho---(7)]]>式(7)中，h_i为冷却盘管内混合工质换热系数，其计算方法如下：对于冷凝管内混合工质中冷凝成分，将其定压降温过程分为过热气冷却至饱和气(g-g)、饱和气冷凝为饱和液(g-1)、饱和液冷却至过冷液(l-1)三个阶段分别进行换热系数计算；(a)过热气体冷却对流换热系数h_i，a采用格林尼斯基公式：对于气体有：ct=(TfTw)0.45(TfTw=0.5~1.5)]]>以及公式：hi=λfdiNuf---(9)]]>计算得到，其中管内流体定性温度tf,i=t1+t22]]>(t₁、t₂分别为冷却器进出口处温度，℃)；管内湍流流动的达尔西阻力系数f＝(1.821g Re-1.64)^-2；(b)相变换热系数h_i，b采用分段热力计算结果，在一定温度间隔内h_i，2采用下面公式计算：hi,bjhvj=QjQj-Qlhj---(10)]]>其中，表示不凝性气体的对流换热系数，其计算方法与上述(a)部分完全相同；为第j温度段内的潜热换热量；(c)过冷液对流换热系数h_i，c：其计算方法与上述(a)类似，区别仅在对于液体有：ct=(PrfPrw)0.11(PrfPrw=0.05~20)]]>(d)根据换热面积关系可以确定管内混合工质平均对流换热系数h_i：QhiΔtm=Qahi,aΔtm,a+Qbhi,bΔtm,b+Qchi,cΔtm,c---(11)]]>其中，Q表示总换热量，Q_a表示(a)过程换热量，Q_b表示(b)过程换热量，Q_c表示(c)过程换热量；Δtm=(T1+T2)/2-TH2O,]]>Δt_m，a表示(a)过程中工质与管外水膜的平均温差，Δt_m，b表示(b)过程中工质与管外水膜的平均温差，Δt_m，c表示(c)过程中工质与管外水膜的平均温差；式(7)中，R_w为冷却盘管导热热阻，采用下式计算得到：Rw=δλ---(12)]]>式(7)中，h_o为管外冷却水的对流换热系数，采用下式计算得到：ho=118(mdo)13---(13)]]>其中，式(7)中，R_i与R_o分别为管内外壁污垢热阻，取R_i＝0.002m²·K/W；(3)冷却水与空气间换热采用传质系数K_oga进行计算，即Koga=3.62×104(ReG)0.9(Ref,o)0.15(do)-2.6(s1do)-1---(14)]]>上式中各系数确定如下：Re_f，o＝4m/u_f、Re_G＝d_oG_m/u_G、并可以计算出以下各参数：s=QmGρGwf,]]>n=Bs1;]]>(4)计算主要结构参数管排数N，则需要对下述两式进行求解，得到Z₁与Z₂：Z1=1Ψ2·ln[b2(i1-i1*)+(Ψ1+b1)(TH2O-T2)b2(i2-i1*)+(Ψ1+b1)(TH2O-T2)]---(15)]]>Z2=1Ψ1·ln[b2(i1-i1*)+(Ψ2+b1)(t1w-T2)b2(i2-i1*)+(Ψ2+b1)(t1w-T2)]---(16)]]>其中，Ψ₁和Ψ₂是下述方程的两个根，Ψ²+(b₁+b₄)Ψ+(b₁·b₄-b₂·b₃)＝0其中，参数由式(7)与式(14)的计算结果以及如下各式确定：b₁＝(a₁+a₃)，a1=K·a′·sQmi·cpi;]]>b₂＝-a₂，a2=Koga·sQmH2OcpH2O;]]>b₃＝αa₃，a3=K·a′·sQmi·cpi;]]>b₄＝αa₂-a₄，a4=koga·sQmG,]]>α=di*dt;]]>(5)判断Z₁与Z₂之差是否小于工程许可值，计算中取|Z₁-Z₂|≤0.02|Z₁+Z₂|作为收敛条件，如果满足收敛条件则用公式N=2Z/3s1]]>得到管排数N，其中Z＝(Z₁+Z₂)/2，进而计算得到外形几何参数；如果不满足收敛条件，重新假定结构设计部分步骤1中的按步骤(2)中(a)～(d)的顺序迭代求解各个变量，直到满足内迭代收敛条件为止；(6)计算冷却/凝盘管进出口压损，利用勃拉休斯公式计算沿程摩擦系数f：f=0.3164Re0.25---(17)]]>利用下面公式计算冷却/凝盘管进出口压损：Δp=(f·l0×Ndi+(N-1)·ξ)·12ρV2---(18)]]>其中，ξ取值为0.27，如果压损Δp大于工质进出口允许压损，则增大开始输入的配风系数，对步骤(1)～(6)重复计算，直到满足外迭代的压损条件为止；8)输出模块将热力计算与结构设计模块计算结果以Excel电子表格形式输出。