多效蒸发工艺设计计算

目录第一章前言§1·1概述`第二章蒸发工艺设计计算§2·1蒸浓液浓度计算§2·2溶液沸点和有效温度差的拟定§2·2·1各效由于溶液的蒸汽压下降所引起的温度差损失/§2·2·2各效由于溶液静压强所因引起的温度差损失§2·2·3由经验不计流体阻力产生压降所引起的温度差损失§2·3加热蒸汽消耗量和各效蒸发水量的计算§2·4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数K的拟定§2·5温差的重新分派与试差计算§2·5·1重新分派各效的有效温度差，§2·5·2反复上述计算环节§2·6计算结果列表第三章NaOH溶液的多效蒸发优化程序部分§3·1具体的拉格朗日乘子法求解过程§3·2程序内部变量说明§3·3程序内容：§3·4程序优化计算结果§3·5优化前后费用比较第四章蒸发器工艺尺寸计算§4·1加热管的选择和管数的初步估计§4·1·1加热管的选择和管数的初步估计§4·1·2循环管的选择§4·1·3加热室直径及加热管数目的拟定§4·1·4分离室直径与高度的拟定§4·2接管尺寸的拟定§4·2·1溶液进出§4·2·2加热蒸气进口与二次蒸汽出口§4·2·3冷凝水出口第五章、蒸发装置的辅助设备§5·1气液分离器§5·2蒸汽冷凝器§5·2·1冷却水量§5·2·2计算冷凝器的直径§5·2·3淋水板的设计§5·3泵选型计算§5·4预热器的选型第六章重要设备强度计算及校核§6·1蒸发分离室厚度设计§6·2加热室厚度校核第七章小结与参考文献：

符号说明希腊字母：c——比热容，KJ/(Kg.h)α――对流传热系数，Ｗ/m2.℃d——管径，mΔ――温度差损失，℃D——直径，mη――误差，D——加热蒸汽消耗量，Kg/hη――热损失系数，f——校正系数，η――阻力系数，F——进料量，Kg/hλ――导热系数，Ｗ/m2.℃g——重力加速度，9.81m/s2μ――粘度，Pa.sh——高度，mρ――密度，Kg/m3H——高度，mk——杜林线斜率K——总传热系数，W/m2.℃∑――加和L——液面高度，mφ――系数L——加热管长度，mL——淋水板间距，m下标：n——效数1,2,3――效数的序号n——第n效0――进料的p——压强，Pai――内侧q——热通量，W/m2m――平均Q——传热速率，Wo――外侧r——汽化潜热，KJ/Kgp――压强R——热阻，m2.℃/Ws――污垢的S——传热面积，m2w――水的t——管心距，mw――壁面的T——蒸汽温度，℃u——流速，m/sU——蒸发强度，Kg/m2.h上标：V——体积流量，m3/h′：二次蒸汽的W——蒸发量，Kg/h′：因溶液蒸汽压而引起的W——质量流量，Kg/h〞：因液柱静压强而引起的x——溶剂的百分质量,％：因流体阻力损失而引起的

第一章前言§1·1概述1蒸发及蒸发流程蒸发是采用加热的方法，使具有不挥发性杂质（如盐类）的溶液沸腾，除去其中被汽化单位部分杂质，使溶液得以浓缩的单元操作过程。蒸发操作广泛用于浓缩各种不挥发性物质的水溶液，是化工、医药、食品等工业中较为常见的单元操作。化工生产中蒸发重要用于以下几种目的：1获得浓缩的溶液产品；2、将溶液蒸发增浓后，冷却结晶，用以获得固体产品，如烧碱、抗生素、糖等产品；3、脱除杂质，获得纯净的溶剂或半成品，如海水淡化。进行蒸发操作的设备叫做蒸发器。蒸发器内要有足够的加热面积，使溶液受热沸腾。溶液在蒸发器内因各处密度的差异而形成某种循环流动，被浓缩到规定浓度后排出蒸发器外。蒸发器内备有足够的分离空间，以除去汽化的蒸汽夹带的雾沫和液滴，或装有适当形式的除沫器以除去液沫，排出的蒸汽如不再运用，应将其在冷凝器中加以冷凝。蒸发过程中经常采用饱和蒸汽间壁加热的方法，通常把作热源用的蒸汽称做一次蒸汽，从溶液蒸发出来的蒸汽叫做而次蒸汽。2蒸发操作的分类按操作的方式可以分为间歇式和连续式，工业上大多数蒸发过程为连续稳定操作的过程。按二次蒸汽的运用情况可以分为单效蒸发和多效蒸发，若产生的二次蒸汽不加运用，直接经冷凝器冷凝后排出，这种操作称为单效蒸发。若把二次蒸汽引至另一操作压力较低的蒸发器作为加热蒸气，并把若干个蒸发器串联组合使用，这种操作称为多效蒸发。多效蒸发中，二次蒸汽的潜热得到了较为充足的运用，提高了加热蒸汽的运用率。按操作压力可以分为常压、加压或减压蒸发。真空蒸发有许多优点：（1）、在低压下操作，溶液沸点较低，有助于提高蒸发的传热温度差，减小蒸发器的传热面积；（2）、可以运用低压蒸气作为加热剂；（3）、有助于对热敏性物料的蒸发；（4）、操作温度低，热损失较小。在加压蒸发中，所得到的二次蒸气温度较高，可作为下一效的加热蒸气加以运用。因此，单效蒸发多为真空蒸发；多效蒸发的前效为加压或常压操作，而后效则在真空下操作。3蒸发操作的特点从上述对蒸发过程的简朴介绍可知，常见的蒸发时间壁两侧分别为蒸气冷凝和液体沸腾的传热过程，蒸发器也就是一种换热器。但和一般的传热过程相比，蒸发操作又有如下特点:(1)

沸点升高蒸发的溶液中具有不挥发性的溶质，在港台压力下溶液的蒸气压较同温度下纯溶剂的蒸气压低，使溶液的沸点高于纯溶液的沸点，这种现象称为溶液沸点的升高。在加热蒸气温度一定的情况下，蒸发溶液时的传热温差必然小于加热唇溶剂的纯热温差，并且溶液的浓度越高，这种影响也越显著。(2)

物料的工艺特性蒸发的溶液自身具有某些特性，例如有些物料在浓缩时也许析出晶体，或易于结垢；有些则具有较大的黏度或较强的腐蚀性等。如何根据物料的特性和工艺规定，选择适宜的蒸发流程和设备是蒸发操作彼此必须要考虑的问题。(3)

节约能源蒸发时汽化的溶剂量较大，需要消耗较大的加热蒸气。如何充足运用热量，提高加热蒸气的运用率是蒸发操作要考虑的另一个问题。4蒸发设备

蒸发设备的作用是使进入蒸发器的原料液被加热，部分气化，得到浓缩的完毕液，同时需要排出二次蒸气，并使之与所夹带的液滴和雾沫相分离。

蒸发的主体设备是蒸发器，它重要由加热室和蒸发室组成。蒸发的辅助设备涉及：使液沫进一步分离的除沫器，和使二次蒸气所有冷凝的冷凝器。减压操作时还需真空装置。兹分述如下：

由于生产规定的不同，蒸发设备有多种不同的结构型式。对常用的间壁传热式蒸发器，按溶液在蒸发器中的运动情况，大体可分为以下两大类：

（1）循环型蒸发器

特点：溶液在蒸发器中做循环流动，蒸发器内溶液浓度基本相同，\o"接近于完毕液浓度：完毕液浓度的大小取决于蒸发器的热流密度和溶液在蒸发器中的停留时间。"接近于完毕液的浓度。操作稳定。此类蒸发器重要有

a.中央循环管式蒸发器，

b.悬筐式蒸发器

c.外热式蒸发器，

d.列文式蒸发器

e.强制循环蒸发器。

其中，\o"溶液的循环流动由密度差引起，由于结构不同，溶液的循环速度不同，传热效果和适应场合也不相同"前四种为自然循环蒸发器。

（2）单程型蒸发器

特点：溶液以液膜的形式一次通过加热室，不进行循环。

优点：溶液停留时间短，故特别合用于热敏性物料的

蒸发；\o"循环型蒸发器加热管中溶液浓度接近于完毕液浓度，单程型蒸发器中的溶液平均浓度接近预料液浓度和完毕液浓度的平均值。由于溶液浓度较低，因此由于沸点升高所引起的温度差损失也较小。"温度差损失较小，表面传热系数较大。

缺陷：设计或操作不妥时不易成膜，热流量将明显下降；不合用于易结晶、结垢物料的蒸发。

此类蒸发器重要有

a.升膜式蒸发器，

b.降膜式蒸发器，

c.刮板式蒸发器

本次设计采用的是中央循环管式蒸发器：结构和原理：其下部的加热室由垂直管束组成，中间由一根直径较大的中央循环管。当管内液体被加热沸腾时，中央循环管内气液混合物的平均密度较大；而其余加热管内气液混合物的平均密度较小。在密度差的作用下，溶液由中央循环管下降，而由加热管上升，做自然循环流动。溶液的循环流动提高了沸腾表面传热系数，强化了蒸发过程。这种蒸发器结构紧凑，制造方便，传热较好，操作可靠等优点，应用十分广泛，有"标准蒸发器"之称。为使溶液有良好的循环，中央循环管的截面积，一般为其余加热管总截面积的40%～100%；加热管的高度一般为1～2m；加热管径多为25～75mm之间。但事实上，由于结构上的限制，其循环速度一般在0.4～0.5m/s以下；蒸发器内溶液浓度始终接近完毕液浓度；清洗和维修也不够方便。

第二章蒸发工艺设计计算§2·1蒸浓液浓度计算多效蒸发的工艺计算的重要依据是物料衡算和、热量衡算及传热速率方程。计算的重要项目有：加热蒸气（生蒸气）的消耗量、各效溶剂蒸发量以及各效的传热面积。计算的已知参数有：料液的流量、温度和浓度，最终完毕液的浓度，加热蒸气的压强和冷凝器中的压强等。蒸发器的设计计算环节多效蒸发的计算一般采用试算法。（1）

根据工艺规定及溶液的性质，拟定蒸发的操作条件（如加热蒸气压强及冷凝器的压强），蒸发器的形式、流程和效数。（2）

根据生产经验数据，初步估计各效蒸发量和各效完毕液的浓度。（3）

根据经验假设蒸气通过各效的压强降相等，估算个效溶液沸点和有效总温差。（4）

根据蒸发器的焓衡算，求各效的蒸发量和传热量。（5）

根据传热速率方程计算各效的传热面积。若求得的各效传热面积不相等，则应按下面介绍的方法重新分派有效温度差，反复环节（3)至（5），直到所求得各效传热面积相等（或满足预先给出的精度规定）为止。F==4557.3kg/h总蒸发量:W=F(1-)=4557.3(1-)=3038.2kg/h并流加料蒸发中无额外蒸汽引出,可设:W2:W3=1:1.1:1.2而W=W1+W2+W3=3038.2kg/h由以上三式可得:W1=920.7kg/h;W2=1012.7kg/h;W3=1104.8kg/h;X1==0.125;X2==0.174;X3=0.3§2·2溶液沸点和有效温度差的拟定设各效间的压强降相等，则总压强差为：=P1-PK/=501.3-30.4=470.9KPaP=式中P---各效加热蒸汽压强与二次蒸气压强之差KPa，----第一次加热蒸气的压强KPa-----末效冷凝器中的二次蒸气的压强KPa各效间的压强差可求得各效蒸发室的压强即P1/=P1-Pi=501.3-470.9/3=344.3KPaP2/=P1-2Pi=501.3-2470.9/3=187.4KPaP3/=Pk/=30.4KPa

由各效的二次蒸汽压强从手册中查得相应的二次蒸汽温度和汽化潜热列与下表中:

第一效第二效第三效二次蒸汽压强Pi/(KPa)

344.3

187.4

30.4二次蒸汽温度Ti/(℃)(即下一效加热蒸汽温度)

138.2

117.9

66.8二次蒸汽的汽化潜热(即下一效加热蒸汽的ri/)

2154.1

2197.5

2332.8多效蒸发中的有效传热总温度差可用下式计算：

有效总温度差式中-----有效总温度差，为各效有效温度差之和，℃。-----第一效加热蒸气的温度，℃。-----冷凝器操作压强下二次蒸气的饱和温度，℃。-------总的温度差损失，为各效温度差损失之和，℃，=/+//+///式中/---由于溶液的蒸汽压下降而引起的温度差损失，℃，//---由于蒸发器红溶液的静压强而引起的温度差损失，℃，///----由于管道流体阻力产生压强降而引起的温度差损失，℃，

§2·2·1各效由于溶液的蒸汽压下降所引起的温度差损失/杜林规则（dnhring’srule）:某种溶液的沸点和相同压强下标准液体（一般为水）的沸点呈线性关系。在以水的沸点为横坐标，该溶液的沸点为纵坐标并以溶液的浓度为参数的直角坐标图上，可得一组直线，称为杜林直线。附录五位NaOH水溶液的杜林线图。运用杜林线图，可根据溶液的浓度及实际压强下水的沸点查出相同压强下溶液的沸点，从而可以得出△/值。根据杜林规则也可以计算液体在各种压强下沸点的近似值。此法的义举世：某液体在两种不同压强下两沸点之差（），与水灾同样压强爱慕；两沸点之差（），其比值为一常数，即（（）/（）=，求得k值，其他任意压强下的沸点可一由下式求得，即=-()；所以，不要杜林线图也可以计算出溶液的△/值。我们的设计是根据Ti/(即相间压强下水的沸点)和各效完毕液浓度xi由附录NaOH水溶液的杜林线图查得各效溶液的沸点tAi分别为:tA1=143.8℃;tA2=126.2℃;tA3=84.8℃;则=143.8-138.2=5.6℃=126.2-117.9=8.3℃=84.8-66.8=18.0℃所以=5.6+8.3+18.0=31.9℃§2·2·2各效由于溶液静压强所因引起的温度差损失由于蒸发器中溶液静压强引起的温度差损失某些蒸发器在操作时，器内溶液需维持一定的液位，因而蒸发器中溶液内部的压强大于液面的压强，致使溶液内部的沸点较液面处的为高，两者之差即为因溶液静压强引起的温度差损失，为简便起见，日夜内部的沸点可按液面和底层的平均压强来查取，平均压强近似按静力学方程估算：pm=p/+式中pm—蒸发器中液面和底层的平均压强，pap/—二次蒸气的压强，即液面处的压强，，pa—溶液的平均密度，-液层高度g-重力加速度，根据pm=p/+取液位高度为1米由NaOH水溶液比重图可得下列数据:NaOH水溶液密度(Kg/m3)Pm1=344.3+=349.6KPaPm2=187.4+=193.0KPaPm3=30.4+=36.7KPa根据各效溶液平均压强查得相应的饱和溶液温度为:T=138.8℃;T=118.9℃;T=72.2℃根据=式中--根据平均压强求取的水的沸点℃，--根据二次蒸气压强求得水的沸点℃所以1=T-T=138.8-138.2=0.6℃2=T-T=118.9-117.9=1.0℃3=T-T=72.2-66.8=5.4℃=0.6+1.0+5.4=7.0℃§2·2·3由经验不计流体阻力产生压降所引起的温度差损失由于管道流体阻力产生的压强降所引起的温度差损失在多效蒸发中末效以前各效的二次蒸汽流到次一效的加热室的过程中由于管道阻力使其压强减少蒸汽的饱和温度也相应减少由此引起的温度差损失即为，根据经验其值可以省略。=0℃根据以估算的各效二次蒸汽压强及温度差损失△，即可由下式估算溶液各效溶液的沸点t所以总的温度差损失为=++=31.9+7.0+0=38.9℃溶液的沸点ti=Ti/+5.6+0.6+0=6.2℃℃℃所以各效溶液沸点：t1=138.2+6.2=144.4℃，t2=117.9+9.3=127.2℃，t3=66.8+23.4=90.2℃§2·3加热蒸汽消耗量和各效蒸发水量的计算第i效的焓衡算式为：有上式可求得第i效的蒸发量.若在焓衡算式计入溶液的能缩热及蒸发器的热损失时，尚需考虑热运用系数一般溶液的蒸发，可取得0.94-0.7△x（式中△x为溶液的浓度变化，以质量分率表达）。第i效的蒸发量的计算式为式中------第i效的加热蒸汽量，当无额外蒸汽抽出时=------第i效加热蒸气的汽化潜热------第i效二次蒸气的汽化潜热-----------原料液的比热---------水的比热，--------分别为第i效及第i-1效溶液的沸点-----------第i效的热运用系数无因次，对于加热蒸气消耗量，可列出各效焓衡算式并与式（3-2）联解而求得。第一效的焓衡量式为:W1=由相关手册查得cp0=3.8644KJ/(Kg.℃)，cpw=4.2035KJ/(Kg.℃)代入W1==0.9225D1(a)同理第二效的热衡算式为:W2=]0.94-0.7(0.174-0.125)=0.9057所以W2=第三效的热衡算式为:W3=W3==0.7456W2-0.0568W1+237.93(c)又W1+W2+W3=3038.2(d)联立(a),(b),(c),(d)式,解得:W1=1057.9Kg/hW2=1032.6Kg/hW3=947.7Kg/hD1=1169.0Kg/h§2·4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数K的拟定任意一效的传热速率方程为Si=式中---第i效的传热速率，W。----第i效的传热系数，W/（m2,℃）.---第i效的传热温度差，℃Si-------第i效的传热面积，m2在三效蒸发中，为了便于制造和安装，通常采用各效传热面积相等的蒸发器，即S1=S3=S3=S若求得的传热面积不相等，应根据各效面积相的原则重新分派各校的有效温度差。方法如下：总的传热系数K；称为公式（A）式中，---管外蒸汽冷凝传热系数与管内液体沸腾传热系数，----管外和管内的污垢热阻，-----管壁的导热系数，----管壁的厚度，，，-----加热管的外径、平均直径与内径，计算K0值重要在于求取关内溶液的沸腾传热系统。该值受溶液的性质、蒸发器的类型、沸腾传热的形式伊基蒸发操作的条件等许多因素的影响。因此，一般沸腾传热膜系数关联式的准确度较差。Q1=D1r1=984.32113=1.94Kg/h=T1-t1=151.8-145.8=6.0℃查得RS0=0.86×10-4m2.K/W,RSi=3.44×10-4m2K/W,取规格为的加热管则d0=38mm,di=38-2.5=33mm,dm=,b=2.5mm第一效：蒸汽在垂直管外的冷凝对流传热系数假设壁温tw1=150℃，特性温度t=（151.8+150）/2=150.9℃，=151.8-150=1.8℃，由t值查得冷凝液的导热系数=0.6839，密度=916.14Kg/m3，粘度=0.1848×10-3，又=2112.9Kj/Kg1==1.13)0.25=10806.6；验算:Re==398.5<2023，所以1=10806.6由经验公式:(*)其中,：---沸腾液体与其蒸汽间的表面压力,N/m;----液体的导热系数,W/m.k;P-----二次蒸汽压力,Pa;cL------液体的比热,J/Kg.K;r-----液体的汽化潜热J/Kg;----液体的密度Kg/m3;-----蒸汽的密度Kg/m3-----液体的动力粘度Pas;-----水的动力粘度Pa.s;-----管壁与沸腾液体间的传热温差K;C-----与表观液面高度比例h有关的系数;h:-----加热管内含纯液柱高度占总管长的比例;对一效查得：=0.07943N/m，P=344.3KPa，C=0.33，cL=3813.1J/Kg.K;=0.7036W/m.k;r=2154.1KJ/Kg;=1090.5Kg/m3;=1.87807Kg/m3;=0.1944×10-3Pa.s;=0.5345×10-3Pas;=150-144.4=5.6K;带入上式解得：1=5903.9带入公式（A）得K1=1082.3检查，带入数据解得tw=149.9℃，与假设基本相同，故不必重算。所以K1=1082.3第二效：假设tw=135℃，特性温度t=（135+138.2）/2=136.6℃，=138.2-135=3.2℃，由t值查得冷凝液的导热系数=0.6853，密度=929.06Kg/m3，粘度=0.2068×10-3Pas，又=2145.1Kj/Kg由公式：2=得2=1.13)0.25=9221.7；验证：Re==529.9<2023，所以2=9221.7由经验公式:由于t=（135+127.2）/2=131.1℃对二效查得：=0.08359N/m，P=187.4KPa，C=0.33，cL=3738.0J/Kg.K;=0.7090W/m.k;r=2197.5KJ/Kg;=1142.4Kg/m3;=1.06027Kg/m3;=0.2233×10-3Pa.s;=0.6925×10-3Pas;=135-127.2=7.8K;带入上式解得：2=5443.8带入公式（A）得K2=1045.7检查，带入数据解得tw=135.3℃，重新假设tw=135.3℃，所以t=(135.3+138.2)/2=136.25℃△

t=138.2-135.3=2.9℃由t=136.25℃查得=0.6853，密度=928.93Kg/m3，粘度=0.2065×10-3Pas，又=2154.1Kj/Kg根据式2=得：2=9454.3；验证：验证：Re==493.1<2023，所以2=9454.1由于t=（135.3+127.2）/2=131.25℃对二效查得：=0.08359N/m，P=187.4KPa，C=0.33，cL=3738.0J/Kg.K;=0.7090W/m.k;r=2197.5KJ/Kg;=1142.4Kg/m3;=1.06027Kg/m3;=0.2233×10-3Pa.s;=0.6925×10-3Pas;=135.3-127.2=8.1K;根据公式：得2=5653.1由公式（A）计算得到K2=1057.3检查带入数据解得tw=135.2℃，与假设基本相同，故不必重算。所以K2=1057.3第三效根据式3=，Re<2023，验证：Re=假设tw3=110℃，特性温度t=（117.9+110）/2=113.95℃，=112.9-110=7.9℃，由t值查得冷凝液的导热系数=0.6846，密度=947.88Kg/m3，粘度=0.2503×10-3Pas，又=2197.5Kj/Kg由公式：3=得3=1.13)0.25=7114.6；验证：Re==817.5<2023，所以2=7114.6由t=（110+90.2）/2=100.1℃对三效查得：=0.0855N/m，P=30.4KPa，C=0.33，cL=3638.30J/Kg.K;=0.7090W/m.k;r=2332.8KJ/Kg;=1275Kg/m3;=0.19326Kg/m3;=0.2832×10-3Pa.s;=1.555×10-3Pas;=110-90.2=19.8K;经验公式:带入上式解得：3=4383.8带入公式（A）得K3=961.6检查带入数据解得tw=110.0℃，与假设基本相同，故不必重算。所以K3=961.6即K1=1082.3,K2=1057.3,K3=961.6。由于参数太多，值查的不准，计算的K值差别不大，优化效果不明显，故用经验值，K1=1860,K2=1280,K3=700分别带入公式：Si=得到S1==m2S2==S3==m2误差为：1--误差较大，故应调整各效的有效温度差，反复上述环节。§2·5温差的重新分派与试差计算§2·5·1重新分派各效的有效温度差，平均传热面积:S==38.26m2重新分派有效温度差:==℃==℃==℃§2·5·2反复上述计算环节：§2·5·2·1由所求各效蒸汽量求各效溶液浓度X1=;X2=，X3=0.3，§2·5·2·2计算各效溶液沸点末效溶液沸点和二次蒸汽压强保持不变,各种温度差损失可视为衡值，故末效溶液的沸点=90.2.而=23.53℃,则第三效加热蒸汽温度（即第二效二次蒸汽温度）为T3=T2/=90.2+23.53=113.75℃,由于第二效二次蒸汽温度为T2/=113.75℃,再由X2=0.185查杜林曲线得:tA2=121.74℃,所以t2=tA2+=121.74+1.0=122.74℃同理：由t2=122.74℃，=12.93℃则T2=T1/=t2+=122.74+12.93=135.67℃再由T1/=135.67℃，X1=0.130，查杜林曲线得:tA1=141.62℃t1=141.62+0.6=142.22℃；也可以由t1=T1=151.8-9.64=142.16℃说明溶液的温度差损失变化不打，不必重新计算，故有效温度差为46.1℃以上计算结果总结如下:效数第一效第二效第三效溶液沸点ti℃142.22122.7490.2加热蒸汽的温度(℃)151.8134.5107.93有效温度差(℃)9.4612.9323.53§2·5·2·3各效焓衡算:二次蒸汽的汽化潜热和二次蒸汽温度Ti/(℃)(即下一效加热蒸汽)温度如下表：效数第一效第二效第三效二次蒸汽温度Ti/(℃)(即下一效加热蒸汽)温度135.67113.7366.8二次蒸汽的汽化潜热(即下一效加热蒸汽的ri/)2161.32222.32332.8第一效：=0.94－0.7×（0.130－0.10）=0.919W1=0.919Kg/h第二效：W2=第三效：=0.94－0.7×(0.30－0.185)=0.8595W3=0.8595[]Kg/hW3=211.14－0.0482W1+07705W2由于：W1+W2+W3=3038.2Kg/h得：W1=1047.1Kg/h，W2=1023.4Kg/hW3=967.7Kg/h，D=1165.4Kg/h与第一次热量恒算所得结果：W1=1057.9Kg/h，W2=1032.6Kg/h，W3=947.7Kg/h比较并计算误差得=0.0103，=0.009，=0.0021相对误差均在误差允许范围之内股计算得各效蒸发面积合理。其各效溶液浓度无明显变化不必再算。§2·5·2·4计算各效传热面积S1==m2S2==m2S3==m2误差计算得：1-=1－36.27/38.15=0.049<0.05；所以误差允许取平均面积的S=37.18m2§2·6计算结果列表效数123冷凝器加热蒸汽温度(℃)151.8135.67113.7366.8操作压强Pi/(KPa)319.6162.630.430.4溶液沸点ti℃142.22122.7490.2

完毕液浓度(%)13.018.530.0

蒸发水量WiKg/h1047.11023.4967.7

生蒸汽量DKg/h1165.4

传热面积Sim237.1837.1837.18

表中Pi/按T1/=T2=135.67℃查得,P2/按T2/=T3=113.73℃查得,Pk/≈P3/=30.4KPa第三章NaOH溶液的多效蒸发优化程序部分§3·1具体的拉格朗日乘子法求解过程具体的拉格朗日乘子法求解过程如下：（1）先输入笔算得相关数据（与优化计算的过程相比可视为常量的或者不变的量）a0=0.018b0=0.000146,B=40000，c0=4180,c1=0.5,c2=1,c=0，e=0.00001,f=4557.3,fc=0.15，k1=1860,k2=1280,k3=700，o=7680,p=1000,R0=2112.9,r3=2332.8,y=0.6,x0=0.1,x3=0.3,t0=20,t3=66.8,T1=151.8,Tk=66.8,Tpm[3]={138.8,118.9,72.2},Ta[3]={143.8,126.2,84.8},T[3]={138.2,117.9,66.8};（2）

根据公式w=f*(1-x0/x3);计算总的蒸发量（3）

根据Tpm[3]，Ta[3]，T[3]德具体数值计算相关量b视为不便的量。（4）

在输入笔算的第一效和第二效蒸气的汽化潜热为程序内部循环作准备（5）

根据公式w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2;D1=w3*r3/R0;计算分别算出各效蒸发量并求出生蒸汽的量（6）

根据公式x1=f*x0/(f-w1);x2=f*x0/(f-w1-w2);进一步求出第一效和第二效的浓度（7）

根据根据拉格朗日乘子法，求出各效温度差损失dt1、dt2、dt3。(8)运用经验公式R1=(2466904.9-1584.3*T1-4.9*T1*T1)/1000;R2=(2466904.9-1584.3*t1-4.9*t1*t1)/1000;求变化后的第一效和第二效蒸气的汽化潜热。当满足条件时退出循环并开始计算面积，同时求出费用，假如不满足条件就自动传递变量第一效和第二效蒸气的汽化潜热，重新进入循环直至满足条件（9）分别计算年生蒸汽费用，蒸发器年折旧维修费用和整孔系统费用。并考虑节省的费用。§3·2程序内部变量说明程序内部变量说明：a0------计算年蒸汽费用时的常量0.018,a--------双精度变量,A[3]----面积的一维变量,b------计算年蒸汽费用时的常量00.000146,b--------双精度变量,B--------蒸发器用不锈钢市场钢材价格,40000,c0--------水的比热容,c1---------电费0.5,c2-------动力费附加系数,c----------双精度变量,,d---------双精度变量,D1-------生蒸汽-双精度变量dt1,dt2,dt3,------分别为温度差变量e------循环控制常量0.00001,E---------双精度变量,f----------溶液流量,fc-------蒸发器念折旧及维修率0.15,h---------双精度变量,H[3]--------增资系数的一维变量,k1,k2,k3-------各效总的传热系数L---------双精度变量,M---------双精度变量,N---------双精度变量,m---------双精度变量,n---------双精度变量,o---------蒸发器年工作时间7680,p---------水的通常密度,R0---------生蒸汽的汽化潜热,r1,r2,r3--------各效汽化潜热的双精度变量,w-------总的蒸发量,w1,w2,w3-------各效蒸发量,vw--------双精度变量,y----------泵的效率,j,j1,j2,j3-------------费用的双精度变量，t0---------------冷凝水的进口温度20,t1,t2,t3-------------各效的溶液的沸点,T1------------生蒸汽的温度,Tk---------------冷凝器的蒸汽温度,R1,R2-------第一效和第二效蒸气的汽化潜热双精度变量x0,x1,x2,x3---------原料液浓度和各效内部的浓度,Tpm[3]---------溶液静平均压强相应的饱和温度的一维变量,Ta[3]-----------由杜林曲线查出的各小溶液沸点的一维变量,T[3]-------------二次蒸汽温度的一维变量;§3·3程序内容：#include<stdio.h>#include<math.h>main(){inti;doublea0=0.018,a,A[3], b0=0.000146,b,B=40000, c0=4180,c1=0.5,c2=1,c=0, d=0,D1,dt1,dt2,dt3, e=0.00001,E,f=4557.3,fc=0.15,h,H[3], k1=1860,k2=1280,k3=700, L,M=0,N,m,n,o=7680,p=1000, R0=2112.9,R1,R2,r1,r2,r3=2332.8, w,w1,w2,w3,vw,y=0.6, j,j1,j2,j3, x0=0.1,x1,x2,x3=0.3, t0=20,t1,t2,t3=66.8, T1=151.8,Tk=66.8,Tpm[3]={138.8,118.9,72.2}, Ta[3]={143.8,126.2,84.8},T[3]={138.2,117.9,66.8};for(i=0;i<3;i++){c+=Ta[i]-T[i];d+=Tpm[i]-T[i];}b=T1-Tk-c-d;scanf("%f,%f",&r1,&r2);w=f*(1-x0/x3);for(;;){w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2;D1=w3*r3/R0;x1=f*x0/(f-w1);x2=f*x0/(f-w1-w2);m=D1*R0/k1;n=w1*r1/k2;h=w2*r2/k3;a=sqrt(m)+sqrt(n)+sqrt(h);dt1=b*sqrt(m)/a;dt2=b*sqrt(n)/a;dt3=b*sqrt(h)/a;R1=(2466904.9-1584.3*T1-4.9*T1*T1)/1000;t1=t3+dt3+dt2;R2=(2466904.9-1584.3*t1-4.9*t1*t1)/1000;if(fabs((R1-r1)/r1)<=e&&fabs((R2-r2)/r2)<=e)break;else{r1=R1;r2=R2;}}t1=T1-dt1;t2=t3+dt3;printf("xiangguanshuju:\n");printf("D1=%f,\nw=%f,w1=%f,w2=%f,\nR1=%f,R2=%f,\ndt1=%f,dt2=%f,dt3=%f,\nT1=%f,t1=%f,t2=%f,t3=%f,\nx0=%f,x0=%f,x1=%f,x2=%f,x3=%f\n", D1,w,w1,w2,w3,R1,R2,dt1,dt2,dt3,T1,t1,t2,t3,x0,x1,x2,x3);printf("xiamianjisuanmianji:\n");A[0]=m/dt1/3.6;A[1]=n/dt2/3.6;A[2]=h/dt3/3.6;printf("mianjishi:\n");printf("A[0]=%f,\nA[1]=%f,\nA[2]=%f",A[0],A[1],A[2]);printf("\njisuanfeiyong:\n");j1=o*D1*(a0+b0*T1)/10000;for(i=0;i<3;i++){if(A[i]<=100)H[i]=1;elseif(A[i]<=200&&A[i]>100) H[i]=1.2; else H[i]=1.5;}for(i=0;i<3;i++)printf("%f,%f\n",A[i],H[i]);for(i=0;i<3;i++)M+=(4400+B-620)*1.2*(0.667+0.0287*A[i])*H[i];j2=fc*M/10000;vw=w3*(r3*1000+c0*5)/p/c0/(Tk-5-t0)/3600;j3=c1*c2*o*(21+1450*vw*vw)*vw*p/102/y/10000;j=j1+j2+j3;printf("j=%f,\nj1=%f,\nj2=%f,\nj3=%f\n",j,j1,j2,j3);L=j1/j;E=j2/j;N=j3/j;printf("L=%f,E=%f,N=%f\n",L,E,N);return;}§3·4程序优化计算结果程序运营结果：D1=1069.216649,w=3038.202300,w1=1068.914989,w2=1000.496284,w3=968.427580,A[0]=28.070734,A[1]=33.838044,A[2]=45.757387,Dt1=12.019197.Dt2=14.488617,Dt3=19.592186,R1=2113.496284；R2=2257.212455；优化后的：总的费用：J=37.469867万元,各种费用所占的比例：年生蒸汽费用：J1=32.980020万元,L=0.880174,蒸发器年折旧维修费用：J2=4.011927万元,E=0.107071真空系统费用：J3=0.477921万元,N=0.012755§3·5优化前后费用比较等面积的各项费用：总的费用：J’=40.6329万元各种费用所占的比例：年生蒸汽费用：J1’=35.9468万元,L’=0.8847蒸发器年折旧维修费用：J2’=4.0988万元E’=0.1009真空系统费用：J3’=0.5873万元N=0.01445两者比较：总的节省费用：J0=J’-J=40.6329-37.469867=3.1631万元，节省费用比例:7.78%年生蒸汽费用：J10=J1’-J1=35.9486-32.980020=2.9668万元,L0=0.9379,年折旧维修费用:J20=J2’-J2=4.0988-4.0119=0.0869万元,E0=0.02747,真空系统费用：J30=J3’-J3=0.5873-0.4779=0.1094万元,N0=0.03459,第四章蒸发器工艺尺寸计算蒸发器的重要结构尺寸（以下均以第一效为计算对象）我们选取的中央循环管式蒸发器的计算方法如下。§4·1加热管的选择和管数的初步估计§4·1·1加热管的选择和管数的初步估计蒸发器的加热管通常选用38*2.5mm无缝钢管。加热管的长度一般为0.6—2m，但也有选用2m以上的管子。管子长度的选择应根据溶液结垢后的难以限度、溶液的起泡性和厂房的高度等因素来考虑，易结垢和易起泡沫溶液的蒸发易选用短管。根据我们的设计任务和溶液性质，我们选用以下的管子。可根据经验我们选取：L=2M，38*2.5mm可以根据加热管的规格与长度初步估计所需的管子数n’，=124（根）式中S=----蒸发器的传热面积，m2，由前面的工艺计算决定（优化后的面积）；d0----加热管外径，m；L---加热管长度，m；因加热管固定在管板上，考虑管板厚度所占据的传热面积，则计算n’时的管长应用（L—0.1）m.§4·1·2循环管的选择循环管的截面积是根据使循环阻力尽量减小的原则考虑的。我们选用的中央循环管式蒸发器的循环管截面积可取加热管总截面积的40%--100%。加热管的总截面积可按n’计算。循环管内径以D1表达，则所以mm对于加热面积较小的蒸发器，应去较大的百分数。选取管子的直径为：循环管管长与加热管管长相同为2m。按上式计算出的D1后应从管规格表中选取的管径相近的标准管，只要n和n’相差不大。循环管的规格一次拟定。循环管的管长与加热管相等，循环管的表面积不计入传热面积中。§4·1·3加热室直径及加热管数目的拟定加热室的内径取决于加热管和循环管的规格、数目及在管板撒谎可以的排列方式。加热管在管板上的排列方式有三角形排列、正方形排列、同心圆排列。根据我们的数据表加以比较我们选用三角形排列式。管心距t为相邻两管中心线之间的距离，t一般为加热管外径的1.25—1.5倍，目前在换热器设计中，管心距的数据已经标准化，只要拟定管子规格，相应的管心距则是定值。我们选用的设计管心距是：拟定加热室内径和加热管数的具体做法是：先计算管束中心线上管数nc，管子安正三角形排列时，nc=1.1*；其中n为总加热管数。初步估计加热室Di=t(nc-1)+2b’,式中b’=(1—1.5)d0.然后由容器公称直径系列，试选一个内径作为加热室内径并以该内径和循环管外景作同心圆，在同心圆的环隙中，按加热管的排列方式和管心距作图。所画的管数n必须大于初值n’，若不满足，应另选一设备内径，重新作图，直至合适。由于加热管的外径为38mm，可取管心距为48mm；以三角形排列计算，b’=(1—1.5)d0=1.5*d0，Di=t(nc-1)+2b’=48*（13-1）+2*38*1.5=690mm，选取加热室壳体内径为800mm鄙厚为12mm；§4·1·4分离室直径与高度的拟定分离室的直径与高度取决于分离室的体积，而分离室的体积又与二次蒸汽的体积流量及蒸发体积强度有关。分离室体积V的计算式为：式中V-----分离室的体积，m3;W-----某效蒸发器的二次蒸汽量，kg/h;P-----某效蒸发器二次蒸汽量，Kg/m3，U-----蒸发体积强度，m3/(m3*s);即每立方米分离室体积每秒产生的二次蒸汽量。一般用允许值为U=1.1~1.5m3/(m3*s)根据由蒸发器工艺计算中得到的各效二次蒸汽量，再从蒸发体积强度U的数值范围内选取一个值，即可由上式算出分离室的体积。一般说来，各效的二次蒸汽量不相同，其密度也不相同，按上式计算得到的分离室体积也不会相同，通常末效体积最大。为方便起见，各效分离室的尺寸可取一致。分离室体积宜取其中较大者。拟定了分离室的体积，其高度与直径符合关系，拟定高度与直径应考虑一下原则：（1）分离室的高度与直径之比H/D=1~2。对于中央循环管式蒸发器，其分离室一般不能小于1.8m，以保证足够的雾沫分离高度。分离室的直径也不能太少，否则二次蒸汽流速过大，导致雾沫夹带现象严重。(2)在条件允许的情况下，分离室的直径尽量与加热室相同，这样可使结构简朴制造方便。（3）高度和直径都适于施工现场的安放。现取分离室中U=1.2m3/（m3*s）；m3。H=1.8m，，D=1.2m§4·2接管尺寸的拟定流体进出口的内径按下式计算式中-----流体的体积流量m3/s；U--------流体的适宜流速m/s，估算出内径后，应从管规格表格中选用相近的标准管。§4·2·1溶液进出口于并流加料的三效蒸发，第一效溶液流量最大，若各效设备尺寸一致的话，根据第一效溶液流量拟定接管。取流体的流速为0.8m/s；所以取ф57X3.5mm规格管。§4·2·2加热蒸气进口与二次蒸汽出口各效结构尺寸一致二次蒸汽体积流量应取各效中较大者。所以取ф76X3.5mm规格管。§4·2·3冷凝水出口冷凝水的排出一般属于液体自然流动，接管直径应由各效加热蒸气消耗量较大者拟定。所以取ф65X3.5mm规格管。

第五章、蒸发装置的辅助设备§5·1气液分离器蒸发操作时，二次蒸汽中夹带大量的液体，虽在分离室得到初步的分离，但是为了防止损失有用的产品或防止污染冷凝液，还需设立气液分离器，以使雾沫中的液体聚集并与二次蒸汽分离，故气液分离器或除沫器。其类型很多，我们选择惯性式除沫器，起工作原理是运用带有液滴的二次蒸汽在忽然改变运动方向时，液滴因惯性作用而与蒸汽分离。在惯性式分离器的重要尺寸可按下列关系拟定：D0=D1；D1：D2：D3=1：1.5：2H=D3h=0.4~0.5D1D0--------二次蒸汽的管径，mD1--------除沫器内管的直径，mD2--------除沫器外管的直径，mD3--------除沫器外壳的直径，mH---------除沫器的总高度，mh---------除沫器内管顶部与器顶的距离，mD1=103mmD2=154.5mmD3=206mmH=206mmh=46.35mm选取二次蒸汽流出管：除雾器内管：除雾器外罩管：§5·2蒸汽冷凝器蒸汽冷凝器的作用是用冷却水将二次蒸汽冷凝。当二次蒸汽为有价值的产品需要回收或会严重地污染冷却水时，应采用间壁式冷却器。当二次蒸汽为水蒸气不需要回收时，可采用直接接触式冷凝器。二次蒸汽与冷凝水直接接触进行热互换，其冷凝效果好，被广乏采用。现采用多孔板式蒸汽冷凝器：

§5·2·1冷却水量根据冷凝器入口蒸汽压强和冷却水进口温度可由图表可查得。-----冷却水量m3/h；-----所需冷凝的蒸汽量，Kg/h实际取§5·2·2由计算可知,进入冷凝器的二次蒸汽的体积流量可计算得到冷凝器的直径D取D=300mm§5·2·3淋水板的设计(1)淋水板数：D<500mm,淋水板取4块(2)淋水板间距：L末0.15m又L0=D+(0.15~0.3)m取L0=0.3+0.2=0.5m。L1=0.6*L0=0.6*0.5=0.3mL2=0.6*L1=0.18mL3=0.6*L2=0.108m<0.15m所以取L0=0.5m。L1=0.3mL2=0.18mL3=0.15m（3）弓型淋水板的宽度：最上面一块=（0.8~0.9）D，其他各淋水板B=D/2+0.05m所以=0.8*0.3~0.9*0.3，取=0.25m，其余取B=0.3/2+0.05+0.2m。(4)淋水板堰高h：D=300mm<500mm,取h=40mm(5)淋水孔径：冷却水循环使用,取d=10mm(6)淋水孔数：淋水孔冷却水流速,,考虑到长期操作易堵,则:N1=n*（1+12%）=134最上一板孔数为N1=n*（1+12%）=134取整为134个其他各板孔数为取整为125个§5·3泵的选型：(1)

真空泵真空度为0.7atm,真空泵排气量---真空系统渗透空气量;---蒸发过程中料液释放的不凝性气体量,常可以忽略;---每小时冷却水能释放的空气量;---蒸发过程中流体的饱和蒸汽压的当量值;---不凝性气体机械夹带的生蒸汽量,常忽略;=0.45kg/h；=0==++)=1.01(0.45+0+0.565)=1.015kg/h=0=++++=.根据排气量和所抽气体基本上是空气的前提条件,选用水喷射泵11/2BA-6系列功率极限真空度材料最大吸气量1.5KW600mmHgA3钢7.4kg/h(2)离心泵q=F/ρ=4557.3/998.2=4.566m3/h查IS型单级单吸离心泵性能表：选型号：IS50-32-125流量：12.5m3/h扬程20m转数：2900r/min轴功率：1.13Kw电机功率：2.2Kw效率60%必须汽蚀余量：2.0m§5·4预热器的选型：将4557.3Kg/h的10%NaOH溶液从20%加热到144.4.4ºC。加热介质采用500Kpa绝压的饱和水蒸气，冷凝液在饱和温度下排出：（1）

拟定流体通入的空间：因加热介质是饱和蒸汽宜通入壳程NaOH溶液通入管内。（2）

拟定流体的定位温度物性数据并选择列管式换热器的型式，饱和蒸汽压强为500Pa,绝压温度T=151.7ºC。因是恒温，饱和蒸汽的定性温度下Tm=151.7ºC，NaOH的定性温度tw=0.5(20+144.4)=82.2ºC.。两流体的温度差Tm-tm=151.7-82.2=69.5ºC两流体的温差大于50ºC，故选用浮式列管换热器。两流体在定性温度下的物性数据如下：物性/流体温度/ºC密度/Kg/m3μ/mPa.sCp/KJ/Kg.ºCλ/W/mºCΓKj/KgNaOH溶液82.21094.70.683.86020.7001

冷凝水120.2

21