多效蒸发工艺设计计算

1、第一章 前言1 - 1 概述、第二章蒸发工艺设计计算2 - 1蒸浓液浓度计算2 - 2溶液沸点和有效温度差的确定2-2-1各效由于溶液的蒸汽压下降所引起的温度差损失/2-2-2 各效由于溶液静压强所因引起的温度差损失2-2-3由经验不计流体阻力产生压降所引起的温度差损失2 - 3加热蒸汽消耗量和各效蒸发水量的计算K的确定2 4 蒸发器的传热面积和有效温度差在各效中的分布以及传热系数2 - 5温差的重新分配与试差计算2 5 1重新分配各效的有效温度差，2-5-2重复上述计算步骤2 - 6计算结果列表第三章NaOH溶液的多效蒸发优化程序部分3 - 1具体的拉格朗日乘子法求解过程3 - 2程序内部变

2、量说明3 - 3程序内容：3 - 4程序优化计算结果3- 5优化前后费用比较蒸发器工艺尺寸计算4444444441 加热管的选择和管数的初步估计1 - 1加热管的选择和管数的初步估计1 - 2循环管的选择1 - 3加热室直径及加热管数目的确定1 - 4分离室直径与高度的确定2接管尺寸的确定2 - 1溶液进出2-2 加热蒸气进口与二次蒸汽出口2 - 3冷凝水出口第五章、蒸发装置的辅助设备5 - 1 气液分离器5 2蒸汽冷凝器5 2 T 冷却水量5-2-2计算冷凝器的直径5-2-3淋水板的设计5 - 3泵选型计算5 . 4预热器的选型第六章主要设备强度计算及校核6-1蒸发分离室厚度设计6-2加热室

3、厚度校核第七章小结与参考文献：符号说明希腊字母：c-比热容，KJ/(Kg.h)0C对流传热半数，W /m2d一管径，mA一温度差损失，cD直径，mTD-一加热蒸汽消耗量，Kg/hT一热损失系数，f一一校正系数，T一阻力系数，F一进料量,Kg/h入导热系数，W /m2. Cg重力加速度，9.81m/s2粘良，Pa.sh高度，mp密度，Kg/m3H局度，mk杜林线斜率K总传热系数，W/m2.CL 液面度，mL加热管长度，mL淋水板间距，mn效数n第n效p 压强，Paq热通量，W/m2Q传热速率，Wr汽化潜热，KJ/KgR热阻，m2.C/WS传热面积，m2t 管心距，mT蒸汽温度，Cu流速，m/sU

4、 烝发强度，Kg/m2.hV 体积流量，m3/hW 蒸发量，Kg/hW 质量流量，Kg/hx 溶剂的百分质量,E加和0 系数下标：1,2,3效数的序号0进料的i 内侧m平均o外侧p压强s污垢的w水的w壁面的上标：:二次蒸汽的:因溶液蒸汽压而引起的：因液柱静压强而引起的:因流体阻力损失而引起的第一章前言1 - 1概述1 蒸发及蒸发流程蒸发是采用加热的方法，使含有不挥发性杂质（如盐类）的溶液沸腾，除去其中被汽化单位部分杂质，使溶液得以浓缩的单元操作过程。蒸发操作广泛用于浓缩各种不挥发性物质的水溶液，是化工、医药、 食品等工业中较为常见的单元操作。化工生产中蒸发主要用于以下几种目的：1 获得浓缩的溶

5、液产品；2、将溶液蒸发增浓后，冷却结晶，用以获得固体产品，如烧碱、抗生素、糖等产品；3、脱除杂质，获得纯净的溶剂或半成品，如海水淡化。进行蒸发操作的设备叫做蒸发器。蒸发器内要有足够的加热面积，使溶液受热沸腾。溶液在蒸发器内因各处密度的差异而形成某种循环流动，被浓缩到规定浓度后排出蒸发器外。蒸发器内备有足够的分离空间，以除去汽化的蒸汽夹带的雾沫和液滴，或装有适当形式的除沫器以除去液沫，排出的蒸汽如不再利用，应将其在冷凝器中加以冷凝。蒸发过程中经常采用饱和蒸汽间壁加热的方法，通常把作热源用的蒸汽称做一次蒸汽，从溶液蒸发出来的蒸汽叫做而次蒸汽。2 蒸发操作的分类按操作的方式可以分为间歇式和连续式，工

6、业上大多数蒸发过程为连续稳定操作的过程。按二次蒸汽的利用情况可以分为单效蒸发和多效蒸发，若产生的二次蒸汽不加利用，直接经冷凝器冷凝后排出，这种操作称为单效蒸发。若把二次蒸汽引至另一操作压力较低的蒸发器作为加热蒸气，并把若干个蒸发器串联组合使用，这种操作称为多效蒸发。多效蒸发中，二次蒸汽的潜热得到了较为充分的利用，提高了加热蒸汽的利用率。按操作压力可以分为常压、加压或减压蒸发。真空蒸发有许多优点：( 1) 、 在低压下操作，溶液沸点较低，有利于提高蒸发的传热温度差，减小蒸发器的传热面积；( 2) 、可以利用低压蒸气作为加热剂；( 3) 、有利于对热敏性物料的蒸发；( 4) 、操作温度低，热损失较

7、小。在加压蒸发中，所得到的二次蒸气温度较高，可作为下一效的加热蒸气加以利用。因此，单效蒸发多为真空蒸发；多效蒸发的前效为加压或常压操作，而后效则在真空下操作。3 蒸发操作的特点从上述对蒸发过程的简单介绍可知，常见的蒸发时间壁两侧分别为蒸气冷凝和液体沸腾的传热过程，蒸发器也就是一种换热器。但和一般的传热过程相比，蒸发操作又有如下特点:(1) (1) 沸点升高蒸发的溶液中含有不挥发性的溶质，在港台压力下溶液的蒸气压较同温度下纯溶剂的蒸气压低，使溶液的沸点高于纯溶液的沸点，这种现象称为溶液沸点的升高。在加热蒸气温度一定的情况下，蒸发溶液时的传热温差必定小于加热唇溶剂的纯热温差，而且溶液的浓度越高，这

8、种影响也越显著。(2) (2) 物料的工艺特性蒸发的溶液本身具有某些特性，例如有些物料在浓缩时可能析出晶体，或易于结垢；有些则具有较大的黏度或较强的腐蚀性等。如何根据物料的特性和工艺要求，选择适宜的蒸发流程和设备是蒸发操作彼此必须要考虑的问题。(3) (3) 节约能源蒸发时汽化的溶剂量较大，需要消耗较大的加热蒸气。如何充分利用热量，提高加热蒸气的利用率是蒸发操作要考虑的另一个问题。(4) 蒸发设备蒸发设备的作用是使进入蒸发器的原料液被加热，部分气化，得到浓缩的完成液，同时需要排出二次蒸气，并使之与所夹带的液滴和雾沫相分离。蒸发的主体设备是蒸发器，它主要由加热室和蒸发室组成。蒸发的辅助设备包括：

9、使液沫进一步分离的除沫器，和使二次蒸气全部冷凝的冷凝器。减压操作时还需真空装置。兹分述如下：由于生产要求的不同，蒸发设备有多种不同的结构型式。对常用的间壁传热式蒸发器，按溶液在蒸发器中的运动情况，大致可分为以下两大类：(5) 1)循环型蒸发器特点： 溶液在蒸发器中做循环流动，蒸发器内溶液浓度基本相同，接近于完成液的浓度。操作稳定。此类蒸发器主要有a.中央循环管式蒸发器，b. 悬筐式蒸发器c.外热式蒸发器，d.列文式蒸发器e.强制循环蒸发器。其中，前四种为自然循环蒸发器。( 2)单程型蒸发器特点：溶液以液膜的形式一次通过加热室，不进行循环。优点：溶液停留时间短，故特别适用于热敏性物料的蒸发； 温

10、度差损失较小，表面传热系数较大。缺点：设计或操作不当时不易成膜，热流量将明显下降；不适用于易结晶、结垢物料的蒸发。此类蒸发器主要有 a.升膜式蒸发器, b.降膜式蒸发器, c.刮板式蒸发器本次设计采用的是中央循环管式蒸发器：结构和原理：其下部的加热室由垂直管束组成，中间由一根直径较大的中央循环管。当管内液体被加热沸腾时， 中央循环管内气液混合物的平均密度较大；而其余加热管内气液混合物的平均密度较小。在密度差的作用下，溶液由中央循环管下降，而由加热管上升，做自然循环流动。溶液的循环流动提高了沸腾表面传热系数，强化了蒸发过程。这种蒸发器结构紧凑，制造方便，传热较好，操作可靠等优点，应用十分广泛，有

11、标准蒸发器”之称。为使溶液有良好的循环，中央循环管的截面积，一般为其余加热管总截面 积的403100%加热管的高度一般为 12m；加热管径多为2575mm之间。但实际上，由 于结构上的限制，其循环速度一般在 0.40.5m/s以下；蒸发器内溶液浓度始终接近完成液 浓度；清洗和维修也不够方便。第二章蒸发工艺设计计算ffi 1并流多效蒸发淡馨浓婚工芭旗程图他 1蒸浓液浓度计算多效蒸发的工艺计算的主要依据是物料衡算和、热量衡算及传热速率方程。 计算的主要项目有：加热蒸气(生蒸气)的消耗量、各效溶剂蒸发量以及各效的传热面积。计算的已知 参数有：料液的流量、温度和浓度，最终完成液的浓度， 加热蒸气的压强

12、和冷凝器中的压强等。蒸发器的设计计算步骤多效蒸发的计算一般采用试算法。(1) (1)根据工艺要求及溶液的性质，确定蒸发的操作条件(如加热蒸气压强及冷凝器的压强)，蒸发器的形式、流程和效数。(2) (2)根据生产经验数据，初步估计各效蒸发量和各效完成液的浓度。(3) (3)根据经验假设蒸气通过各效的压强降相等，估算个效溶液沸点和有效总温差。(4) (4)根据蒸发器的燃衡算，求各效的蒸发量和传热量。(5) (5)根据传热速率方程计算各效的传热面积。若求得的各效传热面积不相等，则应按下面介绍的方法重新分配有效温度差，重复步骤(3)至(5),直到所求得各效传热面积相等(或满足预先给出的精度要求)为止。

13、3.5 104 103F= 320 24=4557.3kg/hXo0.1总蒸发量：W=F (1- X3 )=4557.3 (1- 0.3)=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;F X0FX0X= F -W1 =0.125； X 2= F -W1 -W2 =0.174;X 3=0.3亚2溶液沸点和有效温度差的确定设各效间的压强降相等，则总压强差为：、p /P =Pi-Pk =501.3-30.4=47

14、0.9KPaP -P/K二 P二 n式中 Ap -各效加热蒸汽压强与二次蒸气压强之差KPa,p -第一次加热蒸气的压强KPaP/K-末效冷凝器中的二次蒸气的压强KPa各效间的压强差可求得各效蒸发室的压强即 P1/=P1- ： Pi=501.3-470.9/3=344.3KPaP2/ =P1-2 ： Pi=501.3-2 470.9/3=187.4KPa /P3=Pk =30.4KPa由各效的二次蒸汽压强从手册中查得相应的二次蒸汽温度和汽化潜热列与 下表中：A效第二效第三效二次蒸汽压强Pi/(KPa)344.3187.430.4二次蒸汽温度 Ti/(C)(即下一效加热蒸汽温度)138.2117.

15、966.8二次蒸汽的汽化潜热(即下一效 加热蒸汽的r；)2154.12197.52332.8多效蒸发中的有效传热总温度差可用下式计算：有效总温度差W(T1 一丁K) 一、式中 工At -有效总温度差，为各效有效温度差之和，C。T1-第一效加热蒸气的温度，C。T/K -冷凝器操作压强下二次蒸气的饱和温度，C。工 以总的温度差损失，为各效温度差损失之和，C,.= I/：,/式中 2 A /-由于溶液的蒸汽压下降而引起的温度差损失，c,工一由于蒸发器红溶液的静压强而引起的温度差损失，C,z /-由于管道流体阻力产生压强降而引起的温度差损失，c,8 2 1各效由于溶液的蒸汽压下降所引起的温度差损失/杜

16、林规则（dnhring srue）:某种溶液的沸点和相同压强下标准液体（一般为水） 的沸点呈线性关系。在以水的沸点为横坐标，该溶液的沸点为纵坐标并以溶液的浓度为参数的直角坐标图上，可得一组直线，称为杜林直线。附录五位NaOHK溶液的杜林线图。利用杜林线图，可根据溶液的浓度及实际压强下水的沸点查出相同压强下溶液的沸点，从而可以得出/值。根据杜林规则也可以计算液体在各种压强下沸点的近似值。此法的义举世：某液体在两种不同压强下两沸点之差（tA1 -tA2）,与水灾同样压强喜爱；两沸点之差（ tB1-TB2）, 其比值为一常数，即（tA1 -tA2 ） / （ tB1 -TB2 ） =k,求得k值，其

17、他任意压强下的沸点 tA 可一由下式求得，即t A = tA1 - k （ tB1 TB 2 ）;所以，不要杜林线图也l以计算出溶液的/值。我们的设计是根据 Ti/（即相间压强下水的沸点）和各效完成液浓度 Xi由附录 NaOH水溶液的杜林线图查得各效溶液的沸点tAi分别为：tAi=143.8 C; tA2=126.2 C ; tA3=84.8 C ;贝Ud =143.8-138.2=5.6 C/32 =126.2-117.9=8.3 C ./3 =84.8-66.8=18.0 C/所以 乙 =5.6+8.3+18.0=31.9 C8 2 2各效由于溶液静压强所因引起的温度差损失由于蒸发器中溶液

18、静压强引起的温度差损失“某些蒸发器在操作时，器内溶液需维持一定的液位，因而蒸发器中溶液内部的压强大于液面的压强，致使溶液内部的沸点较液面处的为高，二者之差即为因溶液静压强引起的温度差损失,为简便起见，日夜内部的沸点可按液面和底层的平均压强来查取，平均压强近似按静力学方程估算：DgLpm=p/+2式中Pm一蒸发器中液面和底层的平均压强，pap/一二次蒸气的压强，即液面处的压强，paP一溶液的平均密度，L-液层高度g-重力加速度，；gL根据pm=p/+ 2取液位高度为1米由NaOH水溶液比重图可得下列数据：NaOH 水溶液密度（Kg/m3）E =1089.4,巳=1142.9凡=128371089

19、.4 9.81 13Pm1=344.3+2 10=349.6KPa1142.9 9.81 13Pm2=187.4+2 10=193.0KPa1283.7 9.81 13Pm3 = 30.4+2 10=36.7KPa根据各效溶液平均压强查得对应的饱和溶液温度为：/T pm1=138.8 C ; T pm2=118.9 C; T pm3=72.2 C根据 .-： = tpm -tp式中t pm -根据平均压强求取的水的沸点,tP-根据二次蒸气压强求得水的沸点C/所以产 T pm1- T1 =138.8-138.2=0.6 C/ 2= T pm2- T 2 =118.9-117.9=1.0 C/ 3

20、= T pm3-T 3 =72.2-66.8=5.4 C.-：=0.6+1.0+5.4=7.0 C8 2 3由经验不计流体阻力产生压降所引起的温度差损失由于管道流体阻力产生的压强降所引起的温度差损失在多效蒸发中末效以前各效的 二次蒸汽流到次一效的加热室的过程中由于管道阻力使其压强降低蒸汽的饱和温度也相应 降低由此引起的温度差损失即为 二 根据经验其值可以省略。.：乙 =0C根据以估算的各效二次蒸汽压强t；及温度差损失,即可由下式估算溶液各效溶液的沸点t所以总的温度差损失为J+ 匚 + 乙 =31.9+7.0+0=38.9 C溶液的沸点ti=Ti/+ ：i./. /, /1 = +d=5.6+0

21、.6+0=6.2 C二八；0/ 0/ =8.3+ 1.0+ 0=9.3 c 2222&3=晨0 +A. =18.0+5.4+0 = 23.4 c所以各效溶液沸点：t1=138.2+6.2=144.4 C,t2=117.9+9.3=127.2 C , t3=66.8+23.4=90.2 C8 3加热蒸汽消耗量和各效蒸发水量的计算第i效的始衡算式为：Qi =Dri =(FCp。Wcpw -W2cpw .-WGCpw)-ti)+Wr有上式可求得第i效的蒸发量 四.若在烙衡算式计入溶液的能缩热及蒸发器的热损失时， 尚需考虑热利用系数”一般溶液的蒸发，”可取得0.94-0.7 * (式中 x为溶液的浓度

22、变化，以质量分率表示)。第i效的蒸发量 W的计算式为riti 4 -tiWi = iDi- (FCp0-W1Cpw-W2Cpw -.-WCpw)riri式中Di -第i效的加热蒸汽量，当无额外蒸汽抽出时Di = Wi 1cp0原料液的比热Cpw水的比热t ,分别为 第i效及第i-1效溶液的沸点%第i效的热利用系数无因次，对于加热蒸气消耗量，可列出各效烙衡算式并与式(3-2)联解而求得。iDi = Fcp01第一效的焰衡量式为：Wi=riri:1=0.94-0.7(0.i25-0.i) =0.9225由相关手册查得 Cp0=3.8644KJ/(Kg. C), cpw =4.2035 KJ/(Kg

23、. C )代入R =2ii2.9(a)Wi= ri =0.9225Di 2i54.i 同理第二效的热衡算式为：rt tn2Wi i - (Fcp0 -Wcpw) i /2 W2=r2r22 =0.94-0.7(0.i74-0.i25)=0.9057 所以W2=2i54.ii44.4 -i27.20.9057Wi-(4557.3 3.8644 -Wi 4.2035)- =0.8580叫 i24.852i97.52i97.5第三效的热衡算式为：r3_t2 -t33W2 V (FCp0 训Cpw -W2Cpw) 2T W3=r3r33 =0.94-0.7(0.3-0.i74) = 0.85i8 2i

24、97.5i27.2 -90.20.85i8W2 (4557.3 3.8644 -4.2035Wi -4.2035W2) W3=2332.82332.8=0.7456 W2-0.0568Wi+237.93(c)又 Wi+ W2+ W3=3038.2(d)联立(a),(b),(c),(d)式，解得: Wi=i057.9Kg/h W2=i032.6Kg/h W3=947.7Kg/h Di= ii69.0 Kg/h 8.4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数 K的确定 任意一效的传热速率方程为QiSi= Ki ti式中Qi-第i效的传热速率， WoKi 一第i效的传热系数，W/ (m

25、2, C).At-第i效的传热温度差，CSi-第i效的传热面积，m2在三效蒸发中，为了便于制造和安装，通常采用各效传热面积相等的蒸发器，即 Si=S3=S3=S 若求得的传热面积不相等，应根据各效面积相的原则重新分配各校的有效温度差。方法如下：总的传热系数K;Ko-0, dm di -idi称为公式（A）lP其中，：Ol-沸腾液体与其蒸汽间的表面压力P-二次蒸汽压力,Pa;r液体的汽化潜热J/Kg;4、，一、3r -蒸汽的密度 Kg/m*产-水的动力粘度 Pa.s;,N/m; -L-液体的导热系数，W/ m.k;cl液体的比热，J/Kg.K;3;L-放体的号度 Kg/m匕-液体的动力粘度 Pa

26、 s;*w -管壁与沸腾液体间的传热温差K;式中1a0,必-管外蒸汽冷凝传热系数与管内液体沸腾传热系数RS0, RSi -管外和管内的污垢热阻，九-管壁的导热系数，b-管壁的厚度，do , dm, di- 加热管的外径、平均直径与内径，计算K0值主要在于求取关内溶液的沸腾传热系统。该值受溶液的性质、蒸发器的类型、沸腾传热的形式伊基蒸发操作的条件等许多因素的影响。因此，一般沸腾传热膜系数关联式的准确度较差。6Qi=Di 门=984.3 2113=1.94 10 Kg/h&1 =T1-t1=151.8-145.8=6.0 CdoKo =-;1-LRso 3 4R4：0dmdi查得 Rso=0.86

27、X 10-4m2 .K/W , Rsi=3.44X10-4m2K/W, =17.4W/m K取规格为38父2.5mm的加热管则 d=38mm,d i=38-2.5 2 =33mm,d m= 35.4mm ,b=2.5mm第一效：口。蒸汽在垂直管外的冷凝对流传热系数假设壁温 tW1=150 C,特性温度 t= (151.8+150 ) /2=150.9 C, & =151.8-150=1.8 C,由 t值查得冷凝液的导热系数九=0.6839W/m K ,密度 P=916.14Kg/m 3 ,粘度= =0.1848X10-3,又 r=2112.9Kj/Kg2_3_1_2_3_3g r 4916.1

28、49.81 0.68392112.9 101.13 ()(：0产Jl t =1.132 0.1848 10与 1.8)05=10806.6W/m K ；4： l :t验算：Re= r =398.5 2000 ,所以 4 1=10806.6W / m K 由经验公式：pN*(_)0.5 .( )0.25(*)103 7r口 LC与表观液面高度百分比 h有关的系数； h:-加热管内含纯液柱高度占总管长的百分比对一效查得：0 l =0.07943 N/m , P =344.3K Pa, C=0.33 , CL=3813.1J/Kg.K;l=0.7036W/ m.k; r =2154.1KJ/Kg;P

29、l =1090.5Kg/m 3; Pr =1.87807Kg/m 3;=0.1944 X10-3 Pa.s; l =0.5345 X 10-3 Pa s; At w =150-144.4=5.6K;带入上式解得： :i 1=5903.9W / m K带入公式(A)得 K1=1082.3W /m Kt0 -tw检验twRso二.。tw - t二 iRsi与假设基本相同，故不必重算。所以带入数据解得tw=149.9C,K=1082.3W /m K第二效：假设 tw=135C,特性温度 t= (135+138.2 ) /2=136.6 C, At =138.2-135=3.2 C,由 t 值 查得冷

30、凝液的导热系数 九=0.6853 W/m K ,密度 # =929.06Kg/m 3 ,粘度= =0.2068 X 10-3 Pa s,又 r =2145.1Kj/Kg由公式:-.(- 02=1.13. 2311.13 (3)1%2 =HAt 得929.062 9.81 0.68533 2145.1 103_ 32 0,2068 103.24： l =t)0.25=9221.7W / m K ;验证：Re= r- 由经验公式：=529.9 =1.06027Kg/m3;*=0.2233 X10-3Pa.s; l =0.6925 X 10-3Pa s; tw =135-127.2=7.8 K;带入

31、上式解得： 带入公式（A）-i 2=5443.8W / m K得 K2=1045.7W/m K检验tw重新假设tw=135.3 C,Rsi带入数据解得tw=135.3 ,所以 t=(135.3+138.2)/2=136.25 t=138.2-135.3=2.9由t=136.25 C查得九=0.6853 W / m K ,密度r =2154.1Kj/Kg0 =928.93Kg/m 3 ， 粘度=0.2065 X 10-3Pa s ,又根据式1.13 (-0 2 =23:g r)4T，t得:Ot/02 =9454.3 W /m K ；4.：t验证：验证：Re= 卜=493.12000,所以1 0

32、2=9454.lW/m K因为 t= (135.3+127.2) /2=131.25 C对二效查得：aL =0.08359 N/m , P=187.4 KPa , C=0.33 , CL=3738.0J/Kg.K; =0.7090W/ m.k; r =2197.5KJ/Kg; Pl =1142.4Kg/m 3; 自=1.06027Kg/m 3;=0.2233 X10-3 Pa.s; & =0.6925 X 10-3 Pa s; At w =135.3-127.2=8.1K;二:i；L 小 ：tw cL /L 0.5 0.25二 C(-3)()根据公式：LPr 10 ;?rL得 52=5653.

33、1 W/m K由公式(A计算得到K2=1057.3W /m Ks0si检3tw4%带入数据解得tw=135.2C,与假设基本相同，故不必重算。所以K2=1057.3W/m K第三效根据式1.13 (、23:g r口l，t)44： l ：tRe2000,验证： Re= 出假设 tW3=110C,特性温度 t= (117.9+110 )/2=113.95 C, t =112.9-110=7.9 C ,由 t 值 查得冷凝液的导热系数九=0.6846 W/m K ,密度 * =947.88Kg/m 3 ,粘度= =0.2503 X10-3 Pa s,又 r =2197.5Kj/Kg-231 g r

34、/1.13 (一g)4由公式： a03 =Hit得233947.889.81 0.6846 2197.5 10(；a03=1.132X0.250310 X7.9)0.25=7114.6W/m K ；4： l t验证：Re= r =817.5 2000 ,所以 a0 2=7114.6W / m K由 t= (110+90.2) /2=100.1 C对三效查得：0L =0.0855 N/m , P=30.4 KPa , C=0.33 , CL=3638.30J/Kg.K; l =0.7090W/ m.k; r =2332.8KJ/Kg; Pl =1275Kg/m 3; Pr =0.19326Kg/

35、m 3; =0.2832 X 10-3Pa.s; l =1.555 x 10-3Pa s; & =110-90.2=19.8K;经验公式：LL*.(一)0.25 )i L 0 ：tw CL / ： L 、0.5=C (Q)lPrio3带入上式解得： i 3=4383.8W/m K带入公式(A)得 K3=961.6W/m Kto -tw r Rs0Rsi检3w%带入数据解得tw=iio.oc,与假设基本相同，故不必重算。所以K3=961.6W /m K即 Ki=1082.3 W /m K ,K2=1057.3W/m K ,K3=961.6W /m K 因为参数太多，值查的 不准，计算的 K 值差

36、别不大，优化效果不明显，故用经验值， Ki=1860W/m K ,K2=1280W/ m K ,K3=700W/m K分别带入公式：DrSi= WiK得到D1169.0 2112.9 103/3600 = 49 85S1= t1.K1=1860 (151.8-144.4)- m2D2r21057.9 2154.1 103 /3600244.96mS2=&2 K2 =1280M(138.2-127.2)D3r31032.6 2197.5 103/3600”=32.5S3= I3.K3 =700 (117.9-90.2)m2皿=1_空1 =0.3478误差为：1- Smax49质误差较大，故应调整

37、各效的有效温度差，重复上述步骤。他 5温差的重新分配与试差计算8 5 1重新分配各效的有效温度差，& .大 S2 士 . S3 ：t3平均传热面积：S=2&=38.26m2重新分配有效温度差：S1 .49.85 / - ：t17.4 = 9.64At = S = 38.26CS244.96/ .-：t2-11=12.93S3Ik&2 = S = 38.26C32.51. ：t3 27.7 =23.53=38.26C运 5 2重复上述计算步骤：F X。X= F -W运 5 2 1由所求各效蒸汽量求各效溶液浓度= 0.130；4557.3 0.104557.3 -1057.9X2 =F X0F -

38、W1 -W24557.3 0.104557.3 -1057.9 -1032.6= 0.185X3=0.3 ,85 2 2计算各效溶液沸点末效溶液沸点和二次蒸汽压强保持不变，各种温度差损失可视为衡值，故末效溶液的沸点t3=90.2.而 &3 =23.53 C ,则第三效加热蒸汽温度（即第二效二次蒸汽温 度）为 T 3=T；=90.2+23.53=113.75 C,由于第二效二次蒸汽温度为T2/=113.75 C ,再由X2=0.185查杜林曲线 得:t a2=121.74 C ,/所以 t2= t a2+ &2 =121.74+1.0=122.74 C同理：由 12=122.74C , 2/=1

39、2.93 C /则 T 2= T 1 = t 2+&2 =122.74+12.93=135.67 C再由 T 1/ = 135.67 C, X=0.130 ,查杜林曲线得：t a1=141.62 Ct 1 =141.62+0.6=142.22 C;也可以由t 1 = T 1 -=151.8-9.64=142.16 C说明溶液的温度差损失变化不打，不必重新计算，故有效温度差为46.1 C以上计算结果总结如下效数第一效第二效第三效溶液沸点tc142.22122.7490.2加热蒸汽的温度（C）151.8134.5107.93有效温度差（C）9.4612.9323.538 5 2 3各效始衡算：二次

40、蒸汽的汽化潜热和二次蒸汽温度/（）（即下一效加热蒸汽）温度 如下表:效数第一效第二效第三效二次蒸汽温度Ti/（C）（即下一效加热蒸汽）温度135.67113.7366.8二次蒸汽的汽化潜热（即下一效加热蒸汽的rj）2161.32222.32332.8第一效：“1 =0.94-0.7 X ( 0.130 0.10 ) =0.9192112.9D1W 1=0.919 2161.3 Kg/h第二效：2 =0.94-(0.185-0.130) 0.7 = 0.9015W2 =2161 3142 22 122 740.9015 -W 1+(3864.4 4557.3-4023.5W 1) =139.17

41、+0.8445W 1 Kg/h2222.32222.3-103第三效： 1=0.940.7 X(0.30 -0.185)=0.85952222.3122.74-90.2W1 +(3864.4 4557.3-4023.5W1 -4023.5W2)- W3=0.859 5 2332.82332.8Kg/hW3=211.140.0482 W 1+ 07705W2 因为：W 1+ W2+ W3=3038.2Kg/h得：W 1=1047.1Kg/h ,W2=1023.4Kg/hW3 =967.7 Kg/h ,D=1165.4Kg/h与第一次热量恒算所得结果： W 1=1057.9Kg/h , W2=10

42、32.6Kg/h , W3 =947.7Kg/h 比较并计算误差得 nnn1 =0.0103,2=0.009,3 =0.0021相对误差均在误差允许范围之内股计算得各效蒸发面积合理。其各效溶液浓度无明显变化不必再算。8 5 2 4计算各效传热面积DJiSi= t.K1 =1165.4 2112.9 1 03 / 3600=38.151860 9.64= 37.12 2 mS2=t2,K2W2r2S3=t2.K2 =1280 12.93W1r21023.4 2161.3 103 /3600= 36.27 2967.7 2222.3 103 /3600700 23.53Smin误差计算得： 1-

43、Smax =1 36.27/38.15=0.049 0.05;所以误差允许取平均面积的S=37.18m22 6计算结果列表效数123冷凝器加热蒸汽温度(C)151.8135.67113.7366.8操作压强Pi/(KPa)319.6162.630.430.4溶液沸点t/C142.22122.7490.2完成液浓度()13.018.530.0蒸发水量Wi Kg/h1047.11023.4967.7生蒸汽量D Kg/h1165.4传热面积Si m237.1837.1837.18表中 Pi/ 按 T=T2 =135.67 C查得，P 2 按 T2 =丁3 =113.73 C查得，Pk = Pj =3

44、0.4KPa第三章 NaOH溶液的多效蒸发优化程序部分8. 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,Tpm3=138

45、.8,118.9,72.2, Ta3=143.8,126.2,84.8,T3=138.2,117.9,66.8;(2) 根据公式w=f*(1-x0/x3);计算总的蒸发量(3) (3)根据Tpm3 , Ta3, T3德具体数值计算相关量 b视为不便的量。(4) (4)在输入笔算的第一效和第二效蒸气的汽化潜热为程序内部循环作准备(5) (5)根据公式 w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2;D1=w3*r3/R0;计算分别算出各效蒸发量并求出生蒸汽的量(6) (6)根据公式 x1=f*x0/(f-w1);x2=f*x0/(f-w1

46、-w2); 进一步求出第一效和第二效的浓度(7) (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)分别计算年生蒸汽费用，蒸发器年折旧维修费用和整孔系统费用。并考虑节省的 费用。8 2程序内部变量说明程序内

47、部变量说明：a0-计算年蒸汽费用时的常量0.018, a双精度变量，A3-面积的一维变量，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双精度变量,H3- 增资系数的一维变量,k1,k2,k3 各效总的传热系数L- 双精度变量,M双精度变量,N- 双精度变量,m- 双精度变量,n- 双精度变量,

48、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 原料液浓度和各效内部的浓度,Tpm3- 溶液静平均压强对应的饱和温度的一维变量,Ta3- 由杜林曲线查出的各小溶液沸点的一维变量,T3 二次蒸汽温度的一维变量;3 ,3 程

49、序内容：#include #include main()int i;double a0=0.018,a,A3,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,H3, 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,

50、t1,t2,t3=66.8,T1=151.8,Tk=66.8,Tpm3=138.8,118.9,72.2, Ta3=143.8,126.2,84.8,T3=138.2,117.9,66.8; for(i=0;i3;i+) c+=Tai-Ti;d+=Tpmi-Ti;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-1