多效蒸发设计(以NaOH为例)

1、化 工 原 理 课 程 设 计NaOH 蒸 发 系 统 设 计目 录第一章 前言§1·1 概述第二章 蒸发工艺设计计算§2·1蒸浓液浓度计算§2·2溶液沸点和有效温度差的确定 §2·2·1 各效由于溶液的蒸汽压下降所引起的温度差损失/ §2·2·2 各效由于溶液静压强所因引起的温度差损失§2·2·3 由经验不计流体阻力产生压降所引起的温度差损失§2·3 加热蒸汽消耗量和各效蒸发水量的计算§2·4 蒸发器的传

2、热面积和有效温度差在各效中的分布以及传热系数K的确定§2·5 温差的重新分配与试差计算 §2·5·1重新分配各效的有效温度差，§2·5·2重复上述计算步骤 §2·6计算结果列表第三章 NaOH溶液的多效蒸发优化程序部分§3·1 具体的拉格朗日乘子法求解过程§3·2 程序内部变量说明§3·3 程序内容：§3·4 程序优化计算结果§3·5 优化前后费用比较 第四章 蒸发器工艺尺寸计算§4&#

3、183;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 气液分离器&

4、#167;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校正系数， 阻力系

5、数，F进料量，Kg/h 导热系数，/m2.g重力加速度，9.81m/s2 粘度，Pa.sh高度，m 密度，Kg/m3H高度，mk 杜林线斜率 K总传热系数，W/m2. 加和L液面高度，m 系数L加热管长度，m L淋水板间距，m 下标：n效数1,2,3效数的序号n第n效0进料的p压强，Pai内侧q热通量，W/m2 m平均Q传热速率，W o外侧r汽化潜热，KJ/Kg p压强 R热阻，m2./Ws污垢的S传热面积，m2 w水的t管心距，m w壁面的T蒸汽温度，u流速，m/sU蒸发强度，Kg/m2.h上标：V体积流量，m3/h：二次蒸汽的W蒸发量，Kg/h ：因溶液蒸汽压而引起的W质量流量，Kg/h

6、：因液柱静压强而引起的x溶剂的百分质量,：因流体阻力损失而引起的第一章 前言§1·1概述1蒸发及蒸发流程蒸发是采用加热的方法，使含有不挥发性杂质（如盐类）的溶液沸腾，除去其中被汽化单位部分杂质，使溶液得以浓缩的单元操作过程。蒸发操作广泛用于浓缩各种不挥发性物质的水溶液，是化工、医药、食品等工业中较为常见的单元操作。化工生产中蒸发主要用于以下几种目的：1获得浓缩的溶液产品；2、将溶液蒸发增浓后，冷却结晶，用以获得固体产品，如烧碱、抗生素、糖等产品；3、脱除杂质，获得纯净的溶剂或半成品，如海水淡化。进行蒸发操作的设备叫做蒸发器。蒸发器内要有足够的加热面积，使溶液受热沸腾。溶液在

7、蒸发器内因各处密度的差异而形成某种循环流动，被浓缩到规定浓度后排出蒸发器外。蒸发器内备有足够的分离空间，以除去汽化的蒸汽夹带的雾沫和液滴，或装有适当形式的除沫器以除去液沫，排出的蒸汽如不再利用，应将其在冷凝器中加以冷凝。蒸发过程中经常采用饱和蒸汽间壁加热的方法，通常把作热源用的蒸汽称做一次蒸汽，从溶液蒸发出来的蒸汽叫做而次蒸汽。2蒸发操作的分类按操作的方式可以分为间歇式和连续式，工业上大多数蒸发过程为连续稳定操作的过程。按二次蒸汽的利用情况可以分为单效蒸发和多效蒸发，若产生的二次蒸汽不加利用，直接经冷凝器冷凝后排出，这种操作称为单效蒸发。若把二次蒸汽引至另一操作压力较低的蒸发器作为加热蒸气，并

8、把若干个蒸发器串联组合使用，这种操作称为多效蒸发。多效蒸发中，二次蒸汽的潜热得到了较为充分的利用，提高了加热蒸汽的利用率。按操作压力可以分为常压、加压或减压蒸发。真空蒸发有许多优点：（1）、在低压下操作，溶液沸点较低，有利于提高蒸发的传热温度差，减小蒸发器的传热面积；（2）、可以利用低压蒸气作为加热剂；（3）、有利于对热敏性物料的蒸发；（4）、操作温度低，热损失较小。在加压蒸发中，所得到的二次蒸气温度较高，可作为下一效的加热蒸气加以利用。因此，单效蒸发多为真空蒸发；多效蒸发的前效为加压或常压操作，而后效则在真空下操作。3蒸发操作的特点从上述对蒸发过程的简单介绍可知，常见的蒸发时间壁两侧分别为蒸

9、气冷凝和液体沸腾的传热过程，蒸发器也就是一种换热器。但和一般的传热过程相比，蒸发操作又有如下特点:(1) 沸点升高 蒸发的溶液中含有不挥发性的溶质，在港台压力下溶液的蒸气压较同温度下纯溶剂的蒸气压低，使溶液的沸点高于纯溶液的沸点，这种现象称为溶液沸点的升高。在加热蒸气温度一定的情况下，蒸发溶液时的传热温差必定小于加热唇溶剂的纯热温差，而且溶液的浓度越高，这种影响也越显著。(2) 物料的工艺特性 蒸发的溶液本身具有某些特性，例如有些物料在浓缩时可能析出晶体，或易于结垢；有些则具有较大的黏度或较强的腐蚀性等。如何根据物料的特性和工艺要求，选择适宜的蒸发流程和设备是蒸发操作彼此必须要考虑的问题。(3

10、) 节约能源 蒸发时汽化的溶剂量较大，需要消耗较大的加热蒸气。如何充分利用热量，提高加热蒸气的利用率是蒸发操作要考虑的另一个问题。4蒸发设备   蒸发设备的作用是使进入蒸发器的原料液被加热，部分气化，得到浓缩的完成液，同时需要排出二次蒸气，并使之与所夹带的液滴和雾沫相分离。   蒸发的主体设备是蒸发器，它主要由加热室和蒸发室组成。蒸发的辅助设备包括：使液沫进一步分离的除沫器，和使二次蒸气全部冷凝的冷凝器。减压操作时还需真空装置。兹分述如下：  由于生产要求的不同，蒸发设备有多种不同的结构型式。对常用的间壁传热式蒸发器，按溶液在蒸发器中的运动情况

11、，大致可分为以下两大类：（1）循环型蒸发器   特点：溶液在蒸发器中做循环流动，蒸发器内溶液浓度基本相同，接近于完成液的浓度。操作稳定。此类蒸发器主要有  a.中央循环管式蒸发器，  b.悬筐式蒸发器 c.外热式蒸发器， d.列文式蒸发器 e.强制循环蒸发器。其中，前四种为自然循环蒸发器。（2）单程型蒸发器   特点：溶液以液膜的形式一次通过加热室，不进行循环。  优点：溶液停留时间短，故特别适用于热敏性物料的  蒸发；温度差损失较小，表面传热系数较大。  

12、;缺点：设计或操作不当时不易成膜，热流量将明显下降；不适用于易结晶、结垢物料的蒸发。  此类蒸发器主要有    a.升膜式蒸发器，    b.降膜式蒸发器，    c.刮板式蒸发器本次设计采用的是中央循环管式蒸发器 ：结构和原理：其下部的加热室由垂直管束组成，中间由一根直径较大的中央循环管。当管内液体被加热沸腾时，中央循环管内气液混合物的平均密度较大；而其余加热管内气液混合物的平均密度较小。在密度差的作用下，溶液由中央循环管下降，而由加热管上升，做自然循环

13、流动。溶液的循环流动提高了沸腾表面传热系数，强化了蒸发过程。这种蒸发器结构紧凑，制造方便，传热较好，操作可靠等优点，应用十分广泛，有"标准蒸发器"之称。为使溶液有良好的循环，中央循环管的截面积，一般为其余加热管总截面积的40%100%；加热管的高度一般为12m；加热管径多为2575mm之间。但实际上，由于结构上的限制，其循环速度一般在0.40.5m/s以下；蒸发器内溶液浓度始终接近完成液浓度；清洗和维修也不够方便。第二章 蒸发工艺设计计算§2·1蒸浓液浓度计算多效蒸发的工艺计算的主要依据是物料衡算和、热量衡算及传热速率方程。计算的主要项目有：加热蒸气（生

14、蒸气）的消耗量、各效溶剂蒸发量以及各效的传热面积。计算的已知参数有：料液的流量、温度和浓度，最终完成液的浓度，加热蒸气的压强和冷凝器中的压强等。蒸发器的设计计算步骤多效蒸发的计算一般采用试算法。（1） 根据工艺要求及溶液的性质，确定蒸发的操作条件（如加热蒸气压强及冷凝器的压强），蒸发器的形式、流程和效数。（2） 根据生产经验数据，初步估计各效蒸发量和各效完成液的浓度。（3） 根据经验假设蒸气通过各效的压强降相等，估算个效溶液沸点和有效总温差。（4） 根据蒸发器的焓衡算，求各效的蒸发量和传热量。（5） 根据传热速率方程计算各效的传热面积。若求得的各效传热面积不相等，则应按下面介绍的方法重新分配有

15、效温度差，重复步骤（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

16、.9KPa P=式中 P -各效加热蒸汽压强与二次蒸气压强之差KPa， -第一次加热蒸气的压强KPa-末效冷凝器中的二次蒸气的压强KPa 各效间的压强差可求得各效蒸发室的压强 即P1/=P1- P2/ =P1-2Pi=501.3-2470.9/3=187.4KPa P3/=Pk /=30.4KPa 由各效的二次蒸汽压强从手册中查得相应的二次蒸汽温度和汽化潜热列与下表中: 第一效第二效第三效二次蒸汽压强Pi/(KPa)344.3187.430.4二次蒸汽温度 Ti/()(即下一效加热蒸汽温度)138.2117.966.8二次蒸汽的汽化潜热(即下一效加热蒸汽的ri/)2154.12197.5233

17、2.8 多效蒸发中的有效传热总温度差可用下式计算：有效总温度差 式中 -有效总温度差，为各效有效温度差之和，。-第一效加热蒸气的温度，。-冷凝器操作压强下二次蒸气的饱和温度，。-总的温度差损失，为各效温度差损失之和， =/+/+/式中 /- 由于溶液的蒸汽压下降而引起的温度差损失，/-由于蒸发器红溶液的静压强而引起的温度差损失，/-由于管道流体阻力产生压强降而引起的温度差损失，§2·2·1各效由于溶液的蒸汽压下降所引起的温度差损失/ 杜林规则（dnhringsrule）:某种溶液的沸点和相同压强下标准液体（一般为水）的沸点呈线性关系。在以水的沸点为横坐标，该溶液的

18、沸点为纵坐标并以溶液的浓度为参数的直角坐标图上，可得一组直线，称为杜林直线。附录五位NaOH水溶液的杜林线图。利用 杜林线图，可根据溶液的浓度及实际压强下水的沸点查出相同压强下溶液的沸点，从而可以得出 /值。根据杜林规则也可以计算液体在各种压强下沸点的近似值。此法的义举世：某液体在两种不同压强下两沸点之差（），与水灾同样压强喜爱；两沸点之差（），其比值为一常数，即（）/ （）=，求得k值，其他任意压强下的沸点可一由下式求得，即=-()； 所以，不要杜林线图也可以计算出溶液的 /值。我们的设计是根据Ti/(即相间压强下水的沸点)和各效完成液浓度xi由附录NaOH水溶液的杜林线图查得各效溶液的沸点

19、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各效由于溶液静压强所因引起的温度差损失 由于蒸发器中溶液静压强引起的温度差损失 某些蒸发器在操作时，器内溶液需维持一定的液位，因而蒸发器中溶液内部的压强大于液面的压强，致使溶液内部的沸点较液面处的为高，二者之差即为因溶液静压强引起的温度差损失，为简便起见，日夜内部的沸点可按液面和底层的平均压强来查取，平均压强近似按静力学方程估算：p

20、m=p/+式中 pm蒸发器中 液面和底层的平均压强，pap/二次蒸气的压强，即液面处的压强，pa溶液的平均密度，-液层高度g-重力加速度，根据pm=p/+取液位高度为1米由 NaOH水溶液比重图可得下列数据:NaOH水溶液密度(Kg/m3) Pm1=344.3+=349.6KPa Pm2=187.4+=193.0KPa Pm3 = 30.4+=36.7KPa 根据各效溶液平均压强查得对应的饱和溶液温度为: T=138.8 ; T =118.9; T=72.2根据 = 式中 -根据平均压强求取的水的沸点，-根据二次蒸气压强求得水的沸点 所以1= T- T=138.8-138.2=0.6 2= T

21、 - 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

22、 所以各效溶液沸点：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效的蒸发量.若在焓衡算式计入溶液的能缩热及蒸发器的热损失时，尚需考虑热利用系数一般溶液的蒸发，x（式中x为溶液的浓度变化，以质量分率表示）。 第i效的蒸发量的 计算式为式中 -第i效的加热蒸汽量，当无额外蒸汽抽出时= - 第i效加热蒸气的汽化潜热 -第i效二次蒸气的汽化潜热-原料液的比热 -水的比热，-分别为 第i效及第i-1效溶液的沸点-第i效的热利用系数无因次

23、，对于加热蒸气消耗量，可列出各效焓衡算式并与式（3-2）联解而求得。 第一效的焓衡量式为:W1= 由相关手册查得cp0=3.8644KJ/(Kg.)，cpw =4.2035 KJ/(Kg.)代入 W1= =0.9225D1 (a)同理第二效的热衡算式为: W2=0.94-0.7(0.174-0.125)=0.9057所以W2=第三效的热衡算式为: W3= W3= =0.7456 W2-0.0568W1+237.93 (c) 又W1+ W2+ W3=3038.2 (d) 联立(a),(b),(c),(d)式,解得: W1=1057.9Kg/h W2=1032.6Kg/h W3=947.7Kg/h

24、 D1= 1169.0 Kg/h§2·4蒸发器的传热面积和有效温度差在各效中的分布以及传热系数K的确定 任意一效的传热速率方程为Si=式中 -第i效的传热速率，W。 -第i效的传热系数，W/（m2, ）. -第i效的传热温度差， Si-第i效的传热面积，m2在三效蒸发中，为了便于制造和安装，通常采用各效传热面积相等的蒸发器，即 S1=S3=S3=S若求得的传热面积不相等，应根据各效面积相的原则重新分配各校的有效温度差。方法如下：总的传热系数K； 称为公式（A）式中 ，-管外蒸汽冷凝传热系数与管内液体沸腾传热系数 ，-管外和管内的污垢热阻，-管壁的导热系数，-管壁的厚度， ，

25、-加热管的外径、平均直径与内径，计算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

26、=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 <2000，所以1=10806.6由经验公式: (*) 其中,：- 沸腾液体与其蒸汽间的表面压力,N/m;-液体的导热系数,W/ m.k; P-二次蒸汽压力,Pa; cL-液体的比热 ,J/Kg.K; r -液体的汽化潜热J/Kg; -液体的密度Kg/m3; -蒸汽的密度Kg/m3 - 液体的动力粘度Pa s; -水的动力粘度Pa.s;

27、-管壁与沸腾液体间的传热温差K;C -与表观液面高度百分比h有关的系数;h:-加热管内含纯液柱高度占总管长的百分比; 对一效查得：=0.07943 N/m，P =344.3K Pa，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-3Pa s; =150-144.4=5.6K;带入上式解得： 1=5903.9 带入 公式（A）得K1=1082.3检验 ，带入数据解得tw=149.9，与假设基本相同，故不

28、必重算。所以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-3Pa s，又=2145.1Kj/Kg由公式： 2 =得2=1.13)0.25=9221.7；验证： Re=529.9 <2000，所以2=9221.7由经验公式: 因为t=（135+127.2）/2=131.1对二效查得：=0.08359 N/m，P=187.4 KPa，C=0.33，cL=3738.0J/Kg.K; =0.7090W/

29、 m.k; r =2197.5KJ/Kg; =1142.4Kg/m3; =1.06027Kg/m3; =0.2233×10-3Pa.s;=0.6925×10-3Pa s; =135-127.2=7.8 K;带入上式解得： 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-3Pa s，又=2154.1Kj/Kg根

30、据式 2 =得： 2 =9454.3 ； 验证：验证： Re=493.1<2000，所以2=9454.1因为t=（135.3+127.2）/2=131.25对二效查得：=0.08359 N/m，P=187.4 KPa，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-3Pa s; =135.3-127.2=8.1K;根据公式：得 2=5653.1 由公式（A）计算得到 K2=1057.3检验 带入

31、数据解得tw=135.2，与假设基本相同，故不必重算。所以K2=1057.3第三效根据式 3 = ， Re<2000，验证： Re=假设tw3=110，特性温度t=（117.9+110）/2=113.95，=112.9-110=7.9，由t 值查得冷凝液的导热系数 =0.6846，密度 =947.88Kg/m3 ，粘度=0.2503×10-3Pa s，又=2197.5Kj/Kg由公式： 3 =得3=1.13)0.25=7114.6；验证： Re=817.5 <2000，所以2=7114.6由t=（110+90.2）/2=100.1对三效查得：=0.0855 N/m，P=3

32、0.4 KPa，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-3Pa s; =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值差别不大，优化效果不明显，故用经验值，

33、K1=1860,K2=1280,K3=700分别带入公式： Si=得到 S1=m2 S2= 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计算各

34、效溶液沸点末效溶液沸点和二次蒸汽压强保持不变,各种温度差损失可视为衡值，故末效溶液的沸点=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+

35、0.6=142.22；也可以由 t1 = T1=151.8-9.64=142.16 说明溶液的温度差损失变化不打，不必重新计算，故有效温度差为 46.1 以上计算结果总结如下:效数第一效第二效第三效溶液沸点ti142.22122.7490.2加热蒸汽的温度()151.8134.5107.93有效温度差()9.4612.9323.53§2·5·2·3各效焓衡算: 二次蒸汽的汽化潜热和二次蒸汽温度Ti/() (即下一效加热蒸汽)温度 如下表： 效数第一效第二效第三效二次蒸汽温度Ti/()(即下一效加热蒸汽)温度135.67113.7366.8二次蒸汽的汽化潜

36、热(即下一效加热蒸汽的ri/)2161.32222.32332.8 第一效： =0.940.7×（0.1300.10）=0.919W 1=0.919Kg/h第二效： W2= 第三效： =0.940.7×(0.300.185)=0.8595W3=0.859 5Kg/h W3=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=1032.6Kg/

37、h， W3 =947.7Kg/h比较并计算误差得=0.0103， =0.009， =0.0021相对误差均在误差允许范围之内股计算得各效蒸发面积合理。其各效溶液浓度无明显变化不必再算。§2·5·2·4计算各效传热面积S1= m2S2= = m2 S3= m2 误差计算得： 1-=136.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溶液沸点

38、ti142.22122.7490.2完成液浓度(%)13.018.530.0蒸发水量Wi Kg/h1047.11023.4967.7生蒸汽量D Kg/h1165.4传热面积Si m237.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=40

39、000，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.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） 根据Tpm3，Ta3，T3德具体数值计算相关量b视为不便的量。（4） 在输入笔

40、算的第一效和第二效蒸气的汽化潜热为程序内部循环作准备（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)/100

41、0;求变化后的第一效和第二效蒸气的汽化潜热。当满足条件时退出循环并开始计算面积，同时求出费用，如果不满足条件就自动传递变量第一效和第二效蒸气的汽化潜热，重新进入循环直至满足条件（9）分别计算年生蒸汽费用，蒸发器年折旧维修费用和整孔系统费用。并考虑节省的费用。§3·2 程序内部变量说明程序内部变量说明：a0-计算年蒸汽费用时的常量0.018, a-双精度变量,A3-面积的一维变量, b-计算年蒸汽费用时的常量00.000146,b-双精度变量, B-蒸发器用不锈钢市场钢材价格,40000,c0-水的比热容, c1-电费0.5,c2-动力费附加系数, c-双精度变量,d-双精度

42、变量, D1-生蒸汽-双精度变量dt1,dt2,dt3,-分别为温度差变量 e-循环控制常量0.00001,E-双精度变量, f-溶液流量,fc-蒸发器念折旧及维修率0.15, h-双精度变量,H3-增资系数的一维变量, 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-冷凝水

43、的进口温度20, t1,t2,t3-各效的溶液的沸点,T1-生蒸汽的温度, Tk-冷凝器的蒸汽温度,R1,R2-第一效和第二效蒸气的汽化潜热双精度变量x0,x1,x2,x3-原料液浓度和各效内部的浓度,Tpm3-溶液静平均压强对应的饱和温度的一维变量,Ta3-由杜林曲线查出的各小溶液沸点的一维变量,T3-二次蒸汽温度的一维变量; §3·3程序内容： #include <stdio.h>#include <math.h>main()int i;double a0=0.018,a,A3,b0=0.000146,b,B=40000,c0=4180,c1=0

44、.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,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; fo

45、r(i=0;i<3;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

46、)/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("xiang guan shu ju :n");printf("D1=%

47、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=%fn",D1,w,w1,w2,w3,R1,R2,dt1,dt2,dt3,T1,t1,t2,t3,x0,x1,x2,x3); printf("xia mian ji suan mian ji :n"); A0=m/dt1/3.6; A1=n/dt2/3.6; A2=h/dt3/3.6; printf("mian ji shi :n"

48、); printf("A0=%f,nA1=%f,nA2=%f",A0,A1,A2);printf("nji suan fei yong :n"); j1=o*D1*(a0+b0*T1)/10000; for(i=0;i<3;i+) if(Ai<=100) Hi=1; else if(Ai<=200 && Ai>100) Hi=1.2; else Hi=1.5; for(i=0;i<3;i+) printf("%f,%fn",Ai,Hi);for(i=0;i<3;i+) M+=(440

49、0+B-620)*1.2*(0.667+0.0287*Ai)*Hi; 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=%fn",j,j1,j2,j3); L=j1/j; E=j2/j; N=j3/j; printf("L=%f,E=%f,N=%fn",L,E,N);return;§3·4 程序优化

50、计算结果程序运行结果：D1=1069.216649,w=3038.200000,w1=1068.914989,w2=1000.496284,w3=968.427580,A0=28.070734,A1=33.838044,A2=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 两者比较：总的