naoh 蒸 發(fā) 系 統(tǒng) 設 計

1、<p>　　化 工 原 理 課 程 設 計</p><p>　　NaOH 蒸 發(fā) 系 統(tǒng) 設 計</p><p><b>　　目 錄</b></p><p><b>　　第一章 前言</b></p><p><b>　　§1·1 概述`&l

2、t;/b></p><p>　　第二章 蒸發(fā)工藝設計計算</p><p>　　§2·1蒸濃液濃度計算</p><p>　　§2·2溶液沸點和有效溫度差的確定 </p><p>　　§2·2·1 各效由于溶液的蒸汽壓下降所引起的溫度差損失/</p&g

3、t;<p>　　§2·2·2 各效由于溶液靜壓強所因引起的溫度差損失</p><p>　　§2·2·3 由經(jīng)驗不計流體阻力產(chǎn)生壓降所引起的溫度差損失</p><p>　　§2·3 加熱蒸汽消耗量和各效蒸發(fā)水量的計算</p><p>　　§2·4

4、蒸發(fā)器的傳熱面積和有效溫度差在各效中的分布以及傳熱系數(shù)K的確定</p><p>　　§2·5 溫差的重新分配與試差計算 </p><p>　　§2·5·1重新分配各效的有效溫度差，</p><p>　　§2·5·2重復上述計算步驟 </p><p>　　&#

5、167;2·6計算結(jié)果列表</p><p>　　第三章 NaOH溶液的多效蒸發(fā)優(yōu)化程序部分</p><p>　　§3·1 具體的拉格朗日乘子法求解過程</p><p>　　§3·2 程序內(nèi)部變量說明</p><p>　　§3·3 程序內(nèi)容：</p>&

6、lt;p>　　§3·4 程序優(yōu)化計算結(jié)果</p><p>　　§3·5 優(yōu)化前后費用比較 </p><p>　　第四章 蒸發(fā)器工藝尺寸計算</p><p>　　§4·1 加熱管的選擇和管數(shù)的初步估計 </p><p>　　§4·1

7、3;1 加熱管的選擇和管數(shù)的初步估計</p><p>　　§4·1·2 循環(huán)管的選擇</p><p>　　§4·1·3 加熱室直徑及加熱管數(shù)目的確定</p><p>　　§4·1·4 分離室直徑與高度的確定</p><p>　　§4&#

8、183;2 接管尺寸的確定</p><p>　　§4·2·1 溶液進出</p><p>　　§4·2·2 加熱蒸氣進口與二次蒸汽出口 </p><p>　　§4·2·3 冷凝水出口 </p><p>　　第五章、蒸發(fā)裝置的輔助設備</p

9、><p>　　§5·1 氣液分離器</p><p>　　§5·2 蒸汽冷凝器 </p><p>　　§5·2·1 冷卻水量 </p><p>　　§5·2·2 計算冷凝器的直徑</p><p>　　§

10、5·2·3 淋水板的設計</p><p><b>　　§5·3泵選型計算</b></p><p>　　§5·4預熱器的選型</p><p>　　第六章 主要設備強度計算及校核</p><p>　　§6·1蒸發(fā)分離室厚度設計</p&

11、gt;<p>　　§6·2加熱室厚度校核</p><p>　　第七章 小結(jié)與參考文獻：</p><p><b>　　符號說明</b></p><p><b>　　希臘字母：</b></p><p>　　c——比熱容，KJ/(Kg.h)　　　　　α――對流傳熱系數(shù)，Ｗ

12、/m2.℃　　　</p><p>　　d——管徑，m　　　　　　　　　 Δ――溫度差損失，　℃　　　　　　　</p><p>　　D——直徑，m　　　　　　　　　η――誤差，</p><p>　　D——加熱蒸汽消耗量，Kg/h　　 η――熱損失系數(shù)，</p><p>　　f——校正系數(shù)，　　　　　　　　 η――阻力系數(shù)，</p>

13、;<p>　　F——進料量，Kg/h　　　　　　　 λ――導熱系數(shù)，Ｗ/m2.℃</p><p>　　g——重力加速度，9.81m/s2　　　 μ――粘度，Pa.s</p><p>　　h——高度，m　　　　　　　　　 ρ――密度，Kg/m3</p><p>　　H——高度，m　　　　　　　　　　　　</p>

14、<p>　　k—— 杜林線斜率　　　　　　　　　　 </p><p>　　K——總傳熱系數(shù)，W/m2.℃　　　 ∑――加和</p><p>　　L——液面高度，m　　　　　　　 φ――系數(shù)</p><p>　　L——加熱管長度，m </p><p>　　L——淋水板間距，m

15、 下標：</p><p>　　n——效數(shù)　　　　　　　　　　　　1,2,3――效數(shù)的序號　　　　　　　　　　　　　　　　　　　</p><p>　　n——第n效　　　　　　　　　　　0――進料的　　</p><p>　　p——壓強，Pa　　　　　　　　　　i――內(nèi)側(cè)</p><p>　　q——熱

16、通量，W/m2　　　　　　　 m――平均</p><p>　　Q——傳熱速率，W　　　　　　　 o――外側(cè)</p><p>　　r——汽化潛熱，KJ/Kg　　　　　 p――壓強 　　　</p><p>　　R——熱阻，m2.℃/W　　　　　　　s――污垢的</p><p>　　S——傳熱面積，m2　　　　　　　 w――水的　

17、</p><p>　　t——管心距，m　　　　　　　　　 w――壁面的　</p><p><b>　　T——蒸汽溫度，℃</b></p><p><b>　　u——流速，m/s</b></p><p>　　U——蒸發(fā)強度，Kg/m2.h　　　　　上標：</p><p>　　V

18、——體積流量，m3/h　　　　　　　′：二次蒸汽的</p><p>　　W——蒸發(fā)量，Kg/h　　　　　　 ′：因溶液蒸汽壓而引起的</p><p>　　W——質(zhì)量流量，Kg/h 〞：因液柱靜壓強而引起的</p><p>　　x——溶劑的百分質(zhì)量,％　　　　?。阂蛄黧w阻力損失而引起的</p><p>

19、<b>　　第一章 前言</b></p><p><b>　　§1·1概述</b></p><p><b>　　1蒸發(fā)及蒸發(fā)流程</b></p><p>　　蒸發(fā)是采用加熱的方法，使含有不揮發(fā)性雜質(zhì)（如鹽類）的溶液沸騰，除去其中被汽化單位部分雜質(zhì)，使溶液得以濃縮的單元操作過程

20、。</p><p>　　蒸發(fā)操作廣泛用于濃縮各種不揮發(fā)性物質(zhì)的水溶液，是化工、醫(yī)藥、食品等工業(yè)中較為常見的單元操作?；どa(chǎn)中蒸發(fā)主要用于以下幾種目的：</p><p>　　1獲得濃縮的溶液產(chǎn)品；</p><p>　　2、將溶液蒸發(fā)增濃后，冷卻結(jié)晶，用以獲得固體產(chǎn)品，如燒堿、抗生素、糖等產(chǎn)品；</p><p>　　3、脫除雜質(zhì)，獲得純凈的溶劑

21、或半成品，如海水淡化。進行蒸發(fā)操作的設備叫做蒸發(fā)器。</p><p>　　蒸發(fā)器內(nèi)要有足夠的加熱面積，使溶液受熱沸騰。溶液在蒸發(fā)器內(nèi)因各處密度的差異而形成某種循環(huán)流動，被濃縮到規(guī)定濃度后排出蒸發(fā)器外。蒸發(fā)器內(nèi)備有足夠的分離空間，以除去汽化的蒸汽夾帶的霧沫和液滴，或裝有適當形式的除沫器以除去液沫，排出的蒸汽如不再利用，應將其在冷凝器中加以冷凝。</p><p>　　蒸發(fā)過程中經(jīng)常采用飽和蒸汽

22、間壁加熱的方法，通常把作熱源用的蒸汽稱做一次蒸汽，從溶液蒸發(fā)出來的蒸汽叫做而次蒸汽。</p><p><b>　　2蒸發(fā)操作的分類</b></p><p>　　按操作的方式可以分為間歇式和連續(xù)式，工業(yè)上大多數(shù)蒸發(fā)過程為連續(xù)穩(wěn)定操作的過程。</p><p>　　按二次蒸汽的利用情況可以分為單效蒸發(fā)和多效蒸發(fā)，若產(chǎn)生的二次蒸汽不加利用，直接經(jīng)冷凝器

23、冷凝后排出，這種操作稱為單效蒸發(fā)。若把二次蒸汽引至另一操作壓力較低的蒸發(fā)器作為加熱蒸氣，并把若干個蒸發(fā)器串聯(lián)組合使用，這種操作稱為多效蒸發(fā)。多效蒸發(fā)中，二次蒸汽的潛熱得到了較為充分的利用，提高了加熱蒸汽的利用率。</p><p>　　按操作壓力可以分為常壓、加壓或減壓蒸發(fā)。真空蒸發(fā)有許多優(yōu)點：</p><p>　?。?）、在低壓下操作，溶液沸點較低，有利于提高蒸發(fā)的傳熱溫度差，減小蒸發(fā)器的

24、傳熱面積；</p><p>　?。?）、可以利用低壓蒸氣作為加熱劑；</p><p>　?。?）、有利于對熱敏性物料的蒸發(fā)；</p><p>　?。?）、操作溫度低，熱損失較小。</p><p>　　在加壓蒸發(fā)中，所得到的二次蒸氣溫度較高，可作為下一效的加熱蒸氣加以利用。因此，單效蒸發(fā)多為真空蒸發(fā)；多效蒸發(fā)的前效為加壓或常壓操作，而后效則在真

25、空下操作。</p><p><b>　　3蒸發(fā)操作的特點</b></p><p>　　從上述對蒸發(fā)過程的簡單介紹可知，常見的蒸發(fā)時間壁兩側(cè)分別為蒸氣冷凝和液體沸騰的傳熱過程，蒸發(fā)器也就是一種換熱器。但和一般的傳熱過程相比，蒸發(fā)操作又有如下特點:</p><p>　　沸點升高 蒸發(fā)的溶液中含有不揮發(fā)性的溶質(zhì)，在港臺壓力下溶液的蒸氣壓較同溫度下純?nèi)?/p>

26、劑的蒸氣壓低，使溶液的沸點高于純?nèi)芤旱姆悬c，這種現(xiàn)象稱為溶液沸點的升高。在加熱蒸氣溫度一定的情況下，蒸發(fā)溶液時的傳熱溫差必定小于加熱唇溶劑的純熱溫差，而且溶液的濃度越高，這種影響也越顯著。</p><p>　　物料的工藝特性 蒸發(fā)的溶液本身具有某些特性，例如有些物料在濃縮時可能析出晶體，或易于結(jié)垢；有些則具有較大的黏度或較強的腐蝕性等。如何根據(jù)物料的特性和工藝要求，選擇適宜的蒸發(fā)流程和設備是蒸發(fā)操作彼此必須要考慮

27、的問題。</p><p>　　節(jié)約能源 蒸發(fā)時汽化的溶劑量較大，需要消耗較大的加熱蒸氣。如何充分利用熱量，提高加熱蒸氣的利用率是蒸發(fā)操作要考慮的另一個問題。</p><p>　　4蒸發(fā)設備   蒸發(fā)設備的作用是使進入蒸發(fā)器的原料液被加熱，部分氣化，得到濃縮的完成液，同時需要排出二次蒸氣，并使之與所夾帶的液滴和霧沫相分離。   蒸發(fā)的主體設備是蒸發(fā)

28、器，它主要由加熱室和蒸發(fā)室組成。蒸發(fā)的輔助設備包括：使液沫進一步分離的除沫器，和使二次蒸氣全部冷凝的冷凝器。減壓操作時還需真空裝置。茲分述如下：  由于生產(chǎn)要求的不同，蒸發(fā)設備有多種不同的結(jié)構(gòu)型式。對常用的間壁傳熱式蒸發(fā)器，按溶液在蒸發(fā)器中的運動情況，大致可分為以下兩大類：（1）循環(huán)型蒸發(fā)器   特點：溶液在蒸發(fā)器中做循環(huán)流動，蒸發(fā)器內(nèi)溶液濃度基本相同，接近于完成液的濃度。操作穩(wěn)定。此類蒸發(fā)器主要有

29、  a.中央循環(huán)管式蒸發(fā)器，  b.懸筐式蒸發(fā)器 c.外熱式蒸發(fā)器， d.列文式蒸發(fā)器 e.強制循環(huán)蒸發(fā)器。其中，前四種為自然循環(huán)蒸發(fā)器。（2）單程型蒸發(fā)器   特點：溶液以液膜的形式一次通過加熱室，不進行循環(huán)。  優(yōu)點：溶液停留時間短，故特別適用于熱敏性物料的  蒸發(fā)；溫度差損失較小，表面?zhèn)鳠嵯禂?shù)較大。 &#

30、160;缺點：設計或操作不當時不易成膜，熱流量將明顯下降；不適用于易結(jié)晶、結(jié)垢物料的蒸發(fā)。  </p><p>　　本次設計采用的是中央循環(huán)管式蒸發(fā)器 ：</p><p>　　結(jié)構(gòu)和原理：其下部的加熱室由垂直管束組成，中間由一根直徑較大的中央循環(huán)管。當管內(nèi)液體被加熱沸騰時，中央循環(huán)管內(nèi)氣液混合物的平均密度較大；而其余加熱管內(nèi)氣液混合物的平均密度較小。在密度差的作用下，

31、溶液由中央循環(huán)管下降，而由加熱管上升，做自然循環(huán)流動。溶液的循環(huán)流動提高了沸騰表面?zhèn)鳠嵯禂?shù)，強化了蒸發(fā)過程。</p><p>　　這種蒸發(fā)器結(jié)構(gòu)緊湊，制造方便，傳熱較好，操作可靠等優(yōu)點，應用十分廣泛，有"標準蒸發(fā)器"之稱。為使溶液有良好的循環(huán)，中央循環(huán)管的截面積，一般為其余加熱管總截面積的40%～100%；加熱管的高度一般為1～2m；加熱管徑多為25～75mm之間。但實際上，由于結(jié)構(gòu)上的限制，

32、其循環(huán)速度一般在0.4～0.5m/s以下；蒸發(fā)器內(nèi)溶液濃度始終接近完成液濃度；清洗和維修也不夠方便。</p><p>　　第二章 蒸發(fā)工藝設計計算</p><p>　　§2·1蒸濃液濃度計算</p><p>　　多效蒸發(fā)的工藝計算的主要依據(jù)是物料衡算和、熱量衡算及傳熱速率方程。計算的主要項目有：加熱蒸氣（生蒸氣）的消耗量、各效溶劑蒸發(fā)量

33、以及各效的傳熱面積。計算的已知參數(shù)有：料液的流量、溫度和濃度，最終完成液的濃度，加熱蒸氣的壓強和冷凝器中的壓強等。</p><p>　　蒸發(fā)器的設計計算步驟多效蒸發(fā)的計算一般采用試算法。</p><p>　　根據(jù)工藝要求及溶液的性質(zhì)，確定蒸發(fā)的操作條件（如加熱蒸氣壓強及冷凝器的壓強），蒸發(fā)器的形式、流程和效數(shù)。</p><p>　　根據(jù)生產(chǎn)經(jīng)驗數(shù)據(jù)，初步估計各效蒸發(fā)

34、量和各效完成液的濃度。</p><p>　　根據(jù)經(jīng)驗假設蒸氣通過各效的壓強降相等，估算個效溶液沸點和有效總溫差。</p><p>　　根據(jù)蒸發(fā)器的焓衡算，求各效的蒸發(fā)量和傳熱量。</p><p>　　根據(jù)傳熱速率方程計算各效的傳熱面積。若求得的各效傳熱面積不相等，則應按下面介紹的方法重新分配有效溫度差，重復步驟（3)至（5），直到所求得各效傳熱面積相等（或滿足預先給

35、出的精度要求）為止。</p><p>　　F==4557.3kg/h</p><p>　　總蒸發(fā)量:W=F(1- )=4557.3(1-)=3038.2kg/h</p><p>　　并流加料蒸發(fā)中無額外蒸汽引出,可設:W2:W3=1:1.1:1.2</p><p>　　而W=W1+W2+W3=3038.2kg/h</p><

36、;p>　　由以上三式可得:W1=920.7kg/h;</p><p>　　W2=1012.7kg/h;</p><p>　　W3=1104.8kg/h;</p><p>　　X1==0.125; X2= =0.174;X3=0.3</p><p>　　§2·2溶液沸點和有效溫度差的確定</p><

37、p>　　設各效間的壓強降相等，則總壓強差為：</p><p>　　=P1-PK/=501.3-30.4=470.9KPa</p><p><b>　　P=</b></p><p>　　式中 P ---各效加熱蒸汽壓強與二次蒸氣壓強之差KPa，</p><p>　　----第一次加熱蒸氣的壓強KPa</p&

38、gt;<p>　　-----末效冷凝器中的二次蒸氣的壓強KPa</p><p>　　各效間的壓強差可求得各效蒸發(fā)室的壓強</p><p>　　即P1/=P1-Pi=501.3-470.9/3=344.3KPa</p><p>　　P2/ =P1-2Pi=501.3-2470.9/3=187.4KPa</p><p>　　P3/=

39、Pk /=30.4KPa</p><p>　　由各效的二次蒸汽壓強從手冊中查得相應的二次蒸汽溫度和汽化潛熱列與下表中:</p><p>　　多效蒸發(fā)中的有效傳熱總溫度差可用下式計算：</p><p><b>　　有效總溫度差</b></p><p>　　式中 -----有效總溫度差，為各效有效溫度差之和，℃。<

40、/p><p>　　-----第一效加熱蒸氣的溫度，℃。</p><p>　　-----冷凝器操作壓強下二次蒸氣的飽和溫度，℃。</p><p>　　-------總的溫度差損失，為各效溫度差損失之和，℃，</p><p><b>　　=/+//+///</b></p><p>　　式中 /---

41、由于溶液的蒸汽壓下降而引起的溫度差損失，℃，</p><p>　　//---由于蒸發(fā)器紅溶液的靜壓強而引起的溫度差損失，℃，</p><p>　　///----由于管道流體阻力產(chǎn)生壓強降而引起的溫度差損失，℃，</p><p>　　§2·2·1各效由于溶液的蒸汽壓下降所引起的溫度差損失/</p><p>　　杜林

42、規(guī)則（dnhring’srule）:某種溶液的沸點和相同壓強下標準液體（一般為水）的沸點呈線性關系。在以水的沸點為橫坐標，該溶液的沸點為縱坐標并以溶液的濃度為參數(shù)的直角坐標圖上，可得一組直線，稱為杜林直線。附錄五位NaOH水溶液的杜林線圖。利用 杜林線圖，可根據(jù)溶液的濃度及實際壓強下水的沸點查出相同壓強下溶液的沸點，從而可以得出 △/值。</p><p>　　根據(jù)杜林規(guī)則也可以計算液體在各種壓強下沸點的近似值

43、。此法的義舉世：某液體在兩種不同壓強下兩沸點之差（），與水災同樣壓強喜愛；兩沸點之差（），其比值為一常數(shù)，即（（）/ （）=，求得k值，其他任意壓強下的沸點可一由下式求得，即</p><p><b>　　=-()；</b></p><p>　　所以，不要杜林線圖也可以計算出溶液的 △/值。</p><p>　　我們的設計是根據(jù)Ti/(即相間壓強

44、下水的沸點)和各效完成液濃度xi由附錄</p><p>　　NaOH水溶液的杜林線圖查得各效溶液的沸點tAi分別為:</p><p>　　tA1=143.8℃; tA2=126.2℃; tA3=84.8℃;</p><p>　　則 =143.8-138.2=5.6℃</p><p>　　=126.2-117.9=8.3℃</p&g

45、t;<p>　　=84.8-66.8=18.0℃</p><p>　　所以=5.6+8.3+18.0=31.9℃</p><p>　　§2·2·2各效由于溶液靜壓強所因引起的溫度差損失</p><p>　　由于蒸發(fā)器中溶液靜壓強引起的溫度差損失 某些蒸發(fā)器在操作時，器內(nèi)溶液需維持一定的液位，因而蒸發(fā)器中溶液內(nèi)部的壓強大于液

46、面的壓強，致使溶液內(nèi)部的沸點較液面處的為高，二者之差即為因溶液靜壓強引起的溫度差損失，為簡便起見，日夜內(nèi)部的沸點可按液面和底層的平均壓強來查取，平均壓強近似按靜力學方程估算：</p><p><b>　　pm=p/+</b></p><p>　　式中 pm—蒸發(fā)器中 液面和底層的平均壓強，pa</p><p>　　p/—二次蒸氣的壓強，即液

47、面處的壓強，，pa</p><p><b>　　—溶液的平均密度，</b></p><p><b>　　-液層高度</b></p><p><b>　　g-重力加速度，</b></p><p>　　根據(jù)pm=p/+取液位高度為1米</p><p>　　由

48、 NaOH水溶液比重圖可得下列數(shù)據(jù):</p><p>　　NaOH水溶液密度(Kg/m3) </p><p>　　Pm1=344.3+=349.6KPa</p><p>　　Pm2=187.4+=193.0KPa</p><p>　　Pm3 = 30.4+=36.7KPa</p><p>　　根據(jù)各效溶液平均

49、壓強查得對應的飽和溶液溫度為:</p><p>　　T=138.8℃ ; T =118.9℃; T=72.2℃</p><p><b>　　根據(jù) = </b></p><p>　　式中 --根據(jù)平均壓強求取的水的沸點℃，--根據(jù)二次蒸氣壓強求得水的沸點℃</p><p>　　所以1= T- T=138.8-138.2

50、=0.6℃</p><p>　　2= T - T=118.9-117.9=1.0℃</p><p>　　3= T-T=72.2-66.8=5.4℃</p><p>　　=0.6+1.0+5.4=7.0℃</p><p>　　§2·2·3由經(jīng)驗不計流體阻力產(chǎn)生壓降所引起的溫度差損失</p><p

51、>　　由于管道流體阻力產(chǎn)生的壓強降所引起的溫度差損失在多效蒸發(fā)中末效以前各效的二次蒸汽流到次一效的加熱室的過程中由于管道阻力使其壓強降低蒸汽的飽和溫度也相應降低由此引起的溫度差損失即為，根據(jù)經(jīng)驗其值可以省略。</p><p><b>　　=0℃</b></p><p>　　根據(jù)以估算的各效二次蒸汽壓強及溫度差損失△，即可由下式估算溶液各效溶液的沸點t</p

52、><p>　　所以總的溫度差損失為</p><p>　　=++ =31.9+7.0+0=38.9 ℃</p><p>　　溶液的沸點ti=Ti/+</p><p>　　5.6+0.6+0=6.2℃</p><p><b>　　℃</b></p><p><b>　　

53、℃</b></p><p><b>　　所以各效溶液沸點：</b></p><p>　　t1=138.2+6.2=144.4℃， t2=117.9+9.3=127.2℃ ， t3=66.8+23.4=90.2℃</p><p>　　§2·3加熱蒸汽消耗量和各效蒸發(fā)水量的計算</p><p&g

54、t;　　第i效的焓衡算式為：</p><p>　　有上式可求得第i效的蒸發(fā)量.若在焓衡算式計入溶液的能縮熱及蒸發(fā)器的熱損失時，尚需考慮熱利用系數(shù)一般溶液的蒸發(fā)，可取得0.94-0.7△x（式中△x為溶液的濃度變化，以質(zhì)量分率表示）。 </p><p>　　第i效的蒸發(fā)量的 計算式為</p><p>　　式中 ------第i效的加熱蒸汽量，當無額外蒸汽抽出時= &

55、lt;/p><p>　　------ 第i效加熱蒸氣的汽化潛熱 ------第i效二次蒸氣的汽化潛熱</p><p>　　-----------原料液的比熱 ---------水的比熱</p><p>　　，--------分別為 第i效及第i-1效溶液的沸點</p><p>　　---------

56、--第i效的熱利用系數(shù)無因次，對于加熱蒸氣消耗量，可列出各效焓衡算式并與式（3-2）聯(lián)解而求得。</p><p>　　第一效的焓衡量式為:W1=</p><p>　　由相關手冊查得cp0=3.8644KJ/(Kg.℃)，cpw =4.2035 KJ/(Kg.℃)</p><p><b>　　代入</b></p><p>

57、　　W1= =0.9225D1 (a)</p><p>　　同理第二效的熱衡算式為:</p><p><b>　　W2=]</b></p><p>　　0.94-0.7(0.174-0.125)=0.9057</p><p><b>　　所以</b></p><p>　

58、　W2=第三效的熱衡算式為:</p><p><b>　　W3= </b></p><p><b>　　W3=</b></p><p>　　=0.7456 W2-0.0568W1+237.93 (c)</p><p>　　又W1+ W2+ W3=

59、3038.2 (d)</p><p>　　聯(lián)立(a),(b),(c),(d)式,解得:</p><p>　　W1=1057.9Kg/h</p><p>　　W2=1032.6Kg/h</p><p>　　W3=947.7Kg/h</p>

60、<p>　　D1= 1169.0 Kg/h</p><p>　　§2·4蒸發(fā)器的傳熱面積和有效溫度差在各效中的分布以及傳熱系數(shù)K的確定</p><p>　　任意一效的傳熱速率方程為</p><p><b>　　Si=</b></p><p>　　式中 ---第i效的傳熱速率，W。

61、 ----第i效的傳熱系數(shù)，W/（m2, ℃）.</p><p>　　---第i效的傳熱溫度差，℃ Si-------第i效的傳熱面積，m2</p><p>　　在三效蒸發(fā)中，為了便于制造和安裝，通常采用各效傳熱面積相等的蒸發(fā)器，即 </p><p>　　S1=S3=S3=S</p><p>　　若求得的傳熱面積不相等，應根據(jù)

62、各效面積相的原則重新分配各校的有效溫度差。方法如下：</p><p><b>　　總的傳熱系數(shù)K；</b></p><p><b>　　稱為公式（A）</b></p><p>　　式中 ，---管外蒸汽冷凝傳熱系數(shù)與管內(nèi)液體沸騰傳熱系數(shù) </p><p>　　，---

63、-管外和管內(nèi)的污垢熱阻，-----管壁的導熱系數(shù)，</p><p>　　----管壁的厚度， ，，-----加熱管的外徑、平均直徑與內(nèi)徑，</p><p>　　計算K0值主要在于求取關內(nèi)溶液的沸騰傳熱系統(tǒng) 。該值受溶液的性質(zhì)、蒸發(fā)器的類型、沸騰傳熱的形式伊基蒸發(fā)操作的條件等許多因素的影響。因此，一般沸騰傳熱膜系數(shù)關聯(lián)式的準確度較差。</p><

64、;p>　　Q1=D1r1=984.32113=1.94Kg/h</p><p>　　=T1-t1=151.8-145.8=6.0℃</p><p>　　查得RS0=0.86×10-4m2 .K/W , RSi=3.44×10-4m2K/W,</p><p><b>　　取規(guī)格為的加熱管 </b></p>

65、<p>　　則d0=38mm,di=38-2.5=33mm,dm=,b=2.5mm</p><p><b>　　第一效：</b></p><p>　　蒸汽在垂直管外的冷凝對流傳熱系數(shù)</p><p>　　假設壁溫tw1=150℃，特性溫度t=（151.8+150）/2=150.9℃，=151.8-150=1.8℃，由t 值查得冷凝液的

66、導熱系數(shù) =0.6839，密度 =916.14Kg/m3 ，粘度=0.1848×10-3，又=2112.9Kj/Kg</p><p>　　1==1.13)0.25=10806.6；</p><p>　　驗算:Re==398.5 <2000，所以1=10806.6</p><p><b>　　由經(jīng)驗公式:</b></p&g

67、t;<p><b>　　(*)</b></p><p>　　其中,：--- 沸騰液體與其蒸汽間的表面壓力,N/m;----液體的導熱系數(shù),W/ m.k;</p><p>　　P-----二次蒸汽壓力,Pa; cL------液體的比熱 ,J/Kg.K;</p><p>　　r -----液體

68、的汽化潛熱J/Kg; ----液體的密度Kg/m3; </p><p>　　-----蒸汽的密度Kg/m3 ----- 液體的動力粘度Pa s;</p><p>　　-----水的動力粘度Pa.s; -----管壁與沸騰液體間的傳熱溫差K;</p><p>　　C -----

69、與表觀液面高度百分比h有關的系數(shù);</p><p>　　h:-----加熱管內(nèi)含純液柱高度占總管長的百分比; </p><p>　　對一效查得：=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

70、; =0.1944×10-3Pa.s;=0.5345×10-3Pa s; =150-144.4=5.6K;</p><p>　　帶入上式解得： 1=5903.9 </p><p><b>　　帶入 公式（A）</b></p><p>　　得K1=1082.3</p><p>　　檢驗 ，

71、帶入數(shù)據(jù)解得tw=149.9℃，</p><p>　　與假設基本相同，故不必重算。所以K1=1082.3</p><p><b>　　第二效：</b></p><p>　　假設tw=135℃，特性溫度t=（135+138.2）/2=136.6℃，=138.2-135=3.2℃，由t 值查得冷凝液的導熱系數(shù) =0.6853，密度 =929.06

72、Kg/m3 ，粘度=0.2068×10-3Pa s，又=2145.1Kj/Kg</p><p>　　由公式： 2 =得</p><p>　　2=1.13)0.25=9221.7；</p><p>　　驗證： Re==529.9 <2000，所以2=9221.7</p><p><b>　　由經(jīng)驗公式:<

73、/b></p><p>　　因為t=（135+127.2）/2=131.1℃</p><p>　　對二效查得：=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.

74、s;=0.6925×10-3Pa s; =135-127.2=7.8 K;</p><p>　　帶入上式解得： 2=5443.8 </p><p><b>　　帶入 公式（A）</b></p><p>　　得K2=1045.7</p><p>　　檢驗 ，帶入數(shù)據(jù)解得tw=135.3℃，<

75、/p><p>　　重新假設tw=135.3℃，</p><p>　　所以 t=(135.3+138.2)/2=136.25℃</p><p>　　t=138.2-135.3=2.9℃</p><p>　　由t=136.25℃查得 </p><p>　　=0.6853，密度 =928.93Kg/m3 ，粘度=0.2065

76、×10-3Pa s，又=2154.1Kj/Kg</p><p>　　根據(jù)式 2 =得： 2 =9454.3 ； </p><p>　　驗證：驗證： Re==493.1<2000，</p><p>　　所以2=9454.1</p><p>　　因為t=（135.3+127.2）/2=131.25℃</p&g

77、t;<p>　　對二效查得：=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;</p><

78、p><b>　　根據(jù)公式：</b></p><p>　　得 2=5653.1 </p><p><b>　　由公式（A）</b></p><p>　　計算得到 K2=1057.3</p><p>　　檢驗 帶入數(shù)據(jù)解得tw=135.2℃，</p><p>　　

79、與假設基本相同，故不必重算。所以K2=1057.3</p><p><b>　　第三效</b></p><p>　　根據(jù)式 3 = ， Re<2000，驗證： Re=</p><p>　　假設tw3=110℃，特性溫度t=（117.9+110）/2=113.95℃，=112.9-110=7.9℃，由t 值查得冷凝液的導熱系

80、數(shù) =0.6846，密度 =947.88Kg/m3 ，粘度=0.2503×10-3Pa s，又=2197.5Kj/Kg</p><p>　　由公式： 3 =得</p><p>　　3=1.13)0.25=7114.6；</p><p>　　驗證： Re==817.5 <2000，所以2=7114.6</p><p>

81、　　由t=（110+90.2）/2=100.1℃</p><p>　　對三效查得：=0.0855 N/m，P=30.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-9

82、0.2=19.8K;</p><p><b>　　經(jīng)驗公式:</b></p><p><b>　　帶入上式解得：</b></p><p>　　3=4383.8 </p><p><b>　　帶入 公式（A）</b></p><p><b>

83、　　得K3=961.6</b></p><p>　　檢驗 帶入數(shù)據(jù)解得tw=110.0℃，</p><p>　　與假設基本相同，故不必重算。所以K3=961.6</p><p>　　即K1=1082.3,K2=1057.3,K3=961.6。因為參數(shù)太多，值查的不準，計算的K值差別不大，優(yōu)化效果不明顯，故用經(jīng)驗值，K1=1860,K2=1280,K3

84、=700</p><p><b>　　分別帶入公式：</b></p><p><b>　　Si=得到</b></p><p>　　S1==m2 </p><p><b>　　S2== </b></p><p><b>　　S3= =

85、m2</b></p><p>　　誤差為： 1-- </p><p>　　誤差較大，故應調(diào)整各效的有效溫度差，重復上述步驟。</p><p>　　§2·5溫差的重新分配與試差計算</p><p>　　§2·5·1重新分配各效的有效溫度差，</p><p&

86、gt;　　平均傳熱面積:S==38.26m2</p><p>　　重新分配有效溫度差:</p><p><b>　　==℃</b></p><p><b>　　==℃</b></p><p><b>　　==℃</b></p><p>　　§2

87、·5·2重復上述計算步驟：</p><p>　　§2·5·2·1由所求各效蒸汽量求各效溶液濃度</p><p><b>　　X1=; </b></p><p>　　X2= ， X3=0.3，</p><p>　　§2·5·2&#

88、183;2計算各效溶液沸點</p><p>　　末效溶液沸點和二次蒸汽壓強保持不變,各種溫度差損失可視為衡值，故末效溶液的沸點=90.2.而 =23.53℃,則第三效加熱蒸汽溫度（即第二效二次蒸汽溫度）為 T3=T2/=90.2+23.53=113.75℃, </p><p>　　由于第二效二次蒸汽溫度為T2/=113.75℃,再由X2=0.185查杜林曲線得:tA2=121.74

89、℃,</p><p>　　所以t2= tA2+=121.74+1.0=122.74℃ </p><p>　　同理：由t2=122.74℃ ，=12.93℃</p><p>　　則 T2= T1/ = t2+=122.74+12.93=135.67℃</p><p>　　再由 T1/ =135.67℃，X1=0.130， 查杜林曲

90、線得: tA1=141.62℃</p><p>　　t1 =141.62+0.6=142.22℃；也可以由 t1 = T1=151.8-9.64=142.16 ℃</p><p>　　說明溶液的溫度差損失變化不打，不必重新計算，故有效溫度差為 46.1 ℃ </p><p>　　以上計算結(jié)果總結(jié)如下:</p><p>　　§2&

91、#183;5·2·3各效焓衡算: </p><p>　　二次蒸汽的汽化潛熱和二次蒸汽溫度Ti/(℃) (即下一效加熱蒸汽)溫度 如下表： </p><p>　　第一效： =0.94－0.7×（0.130－0.10）=0.919</p><p>　　W 1=0.919Kg/h</p><p>　　第二效：

92、 </p><p>　　W2= 第三效： =0.94－0.7×(0.30－0.185)=0.8595</p><p>　　W3=0.859 5[]Kg/h</p><p>　　W3=211.14－0.0482 W 1+ 07705W2</p><p>　　因為： W 1+ W2+ W3=3038.2Kg

93、/h</p><p>　　得： W 1=1047.1Kg/h ， W2=1023.4Kg/h</p><p>　　W3 =967.7 Kg/h ， D=1165.4Kg/h</p><p>　　與第一次熱量恒算所得結(jié)果：</p><p>　　W 1=1057.9Kg/h，W2=1032.6Kg/h， W3 =94

94、7.7Kg/h</p><p><b>　　比較并計算誤差得</b></p><p>　　=0.0103， =0.009， =0.0021</p><p>　　相對誤差均在誤差允許范圍之內(nèi)股計算得各效蒸發(fā)面積合理。其各效溶液濃度無明顯變化不必再算。</p><p>　　§2·5·2

95、·4計算各效傳熱面積</p><p><b>　　S1== m2</b></p><p><b>　　S2= = m2</b></p><p><b>　　S3== m2</b></p><p>　　誤差計算得： 1-=1－36.27/38.15=0.049<

96、;0.05；</p><p>　　所以誤差允許取平均面積的S=37.18m2</p><p>　　§2·6計算結(jié)果列表</p><p>　　表中Pi/ 按T1/ =T2 =135.67℃查得, P2/按T2/ =T3 =113.73 ℃查得, Pk/≈P3/ =30.4KPa</p><p>　　第三章 NaOH溶液

97、的多效蒸發(fā)優(yōu)化程序部分</p><p>　　§3·1 具體的拉格朗日乘子法求解過程</p><p>　　具體的拉格朗日乘子法求解過程如下：</p><p>　?。?）先輸入筆算得相關數(shù)據(jù)（與優(yōu)化計算的過程相比可視為常量的或者不變的量）a0=0.018</p><p>　　b0=0.000146,B=40000，c0=418

98、0,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,</p><p>　　R0=2112.9,r3=2332.8,y=0.6,</p><p>　　x0=0.1, x3=0.3,t0=20, t3=66.8,T1=151.8,Tk=66.8,Tpm[3]={138.8,118.9

99、,72.2},</p><p>　　Ta[3]={143.8,126.2,84.8},T[3]={138.2,117.9,66.8};</p><p>　　根據(jù)公式w=f*(1-x0/x3);計算總的蒸發(fā)量</p><p>　　根據(jù)Tpm[3]，Ta[3]，T[3]德具體數(shù)值計算相關量b視為不便的量。</p><p>　　在輸入筆算的第一效和

100、第二效蒸氣的汽化潛熱為程序內(nèi)部循環(huán)作準備</p><p>　　根據(jù)公式w1=w/(1+r1/r2+r1/r3);w2=w/(1+r2/r1+r2/r3);w3=w-w1-w2; D1=w3*r3/R0;計算分別算出各效蒸發(fā)量并求出生蒸汽的量</p><p>　　根據(jù)公式x1=f*x0/(f-w1);x2=f*x0/(f-w1-w2);進一步求出第一效和第二效的濃度</p>&

101、lt;p>　　根據(jù)根據(jù)拉格朗日乘子法，求出各效溫度差損失dt1、dt2、dt3。</p><p>　　(8)利用經(jīng)驗公式R1=(2466904.9-1584.3*T1-4.9*T1*T1)/1000; R2=(2466904.9-1584.3*t1-4.9*t1*t1)/1000;求變化后的第一效和第二效蒸氣的汽化潛熱。當滿足條件時退出循環(huán)并開始計算面積，同時求出費用，如果不滿足條件就自動傳遞變量第一效和第

102、二效蒸氣的汽化潛熱，重新進入循環(huán)直至滿足條件</p><p>　　（9）分別計算年生蒸汽費用，蒸發(fā)器年折舊維修費用和整孔系統(tǒng)費用。并考慮節(jié)省的費用。</p><p>　　§3·2 程序內(nèi)部變量說明</p><p><b>　　程序內(nèi)部變量說明：</b></p><p>　　a0------計算年蒸汽

103、費用時的常量0.018, a--------雙精度變量,</p><p>　　A[3]----面積的一維變量, b------計算年蒸汽費用時的常量00.000146,</p><p>　　b--------雙精度變量, B--------蒸發(fā)器用不銹鋼市場鋼材價格,40000,</p><p

104、>　　c0--------水的比熱容, c1---------電費0.5,</p><p>　　c2-------動力費附加系數(shù), c----------雙精度變量,,</p><p>　　d---------雙精度變量, D1-------生蒸汽-雙精度變量</p&g

105、t;<p>　　dt1,dt2,dt3,------分別為溫度差變量 e------循環(huán)控制常量0.00001,</p><p>　　E---------雙精度變量, f----------溶液流量,</p><p>　　fc-------蒸發(fā)器念折舊及維修率0.15, h---------雙精度變量,&

106、lt;/p><p>　　H[3]--------增資系數(shù)的一維變量, k1,k2,k3-------各效總的傳熱系數(shù)</p><p>　　L---------雙精度變量, M---------雙精度變量, </p><p>　　N---------雙精度變量, m----

107、-----雙精度變量,</p><p>　　n---------雙精度變量, o---------蒸發(fā)器年工作時間7680,</p><p>　　p---------水的通常密度, R0---------生蒸汽的汽化潛熱,</p><p>　　r1,r2,r3--------各效汽化潛熱

108、的雙精度變量, w-------總的蒸發(fā)量,</p><p>　　w1,w2,w3-------各效蒸發(fā)量, vw --------雙精度變量,</p><p>　　y----------泵的效率, j,j1,j2,j3-------------費用的雙精度變量，</p><p>　

109、　t0---------------冷凝水的進口溫度20, t1,t2,t3-------------各效的溶液的沸點,</p><p>　　T1------------生蒸汽的溫度, Tk---------------冷凝器的蒸汽溫度,</p><p>　　R1,R2-------第一效和第二效蒸氣的汽化潛熱雙精度變量</p>

110、<p>　　x0,x1,x2,x3---------原料液濃度和各效內(nèi)部的濃度,</p><p>　　Tpm[3]---------溶液靜平均壓強對應的飽和溫度的一維變量,</p><p>　　Ta[3]-----------由杜林曲線查出的各小溶液沸點的一維變量,</p><p>　　T[3]-------------二次蒸汽溫度的一維變量;</

111、p><p><b>　　§3·3程序內(nèi)容：</b></p><p>　　#include <stdio.h></p><p>　　#include <math.h></p><p><b>　　main()</b></p><p>&l

112、t;b>　　{int i;</b></p><p>　　double a0=0.018,a,A[3],</p><p>　　b0=0.000146,b,B=40000,</p><p>　　c0=4180,c1=0.5,c2=1,c=0,</p><p>　　d=0,D1,dt1,dt2,dt3,</p>&

113、lt;p>　　e=0.00001,E,f=4557.3,fc=0.15,h,H[3],</p><p>　　k1=1860,k2=1280,k3=700,</p><p>　　L,M=0,N,m,n,o=7680,p=1000,</p><p>　　R0=2112.9,R1,R2,r1,r2,r3=2332.8,</p><p>　　w

114、,w1,w2,w3,vw,y=0.6,</p><p>　　j,j1,j2,j3,</p><p>　　x0=0.1,x1,x2,x3=0.3,</p><p>　　t0=20,t1,t2,t3=66.8,</p><p>　　T1=151.8,Tk=66.8,Tpm[3]={138.8,118.9,72.2},</p><

115、;p>　　Ta[3]={143.8,126.2,84.8},T[3]={138.2,117.9,66.8};</p><p>　　for(i=0;i<3;i++)</p><p><b>　　{</b></p><p>　　c+=Ta[i]-T[i];</p><p>　　d+=Tpm[i]-T[i];<

116、;/p><p><b>　　}</b></p><p>　　b=T1-Tk-c-d;</p><p>　　scanf("%f,%f",&r1,&r2);</p><p>　　w=f*(1-x0/x3);</p><p>　　for( ; ; )</p>

117、<p><b>　　{</b></p><p>　　w1=w/(1+r1/r2+r1/r3);</p><p>　　w2=w/(1+r2/r1+r2/r3);</p><p>　　w3=w-w1-w2;</p><p>　　D1=w3*r3/R0;</p><p>　　x1=f*x0/

118、(f-w1);</p><p>　　x2=f*x0/(f-w1-w2);</p><p>　　m=D1*R0/k1;</p><p>　　n=w1*r1/k2;</p><p>　　h=w2*r2/k3;</p><p>　　a=sqrt(m)+sqrt(n)+sqrt(h);</p><p>

119、　　dt1=b*sqrt(m)/a;</p><p>　　dt2=b*sqrt(n)/a;</p><p>　　dt3=b*sqrt(h)/a;</p><p>　　R1=(2466904.9-1584.3*T1-4.9*T1*T1)/1000;</p><p>　　t1=t3+dt3+dt2;</p><p>　　R

120、2=(2466904.9-1584.3*t1-4.9*t1*t1)/1000;</p><p>　　if(fabs((R1-r1)/r1)<=e&& fabs((R2-r2)/r2)<=e)</p><p><b>　　break;</b></p><p><b>　　else</b></

121、p><p><b>　　{r1=R1;</b></p><p><b>　　r2=R2;}</b></p><p><b>　　}</b></p><p>　　t1=T1-dt1;</p><p>　　t2=t3+dt3;</p><p&

122、gt;　　printf("xiang guan shu ju :\n");</p><p>　　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",</p&g

123、t;<p>　　D1,w,w1,w2,w3,R1,R2,dt1,dt2,dt3,T1,t1,t2,t3,x0,x1,x2,x3);</p><p>　　printf("xia mian ji suan mian ji :\n");</p><p>　　A[0]=m/dt1/3.6;</p><p>　　A[1]=n/dt2/3

124、.6;</p><p>　　A[2]=h/dt3/3.6;</p><p>　　printf("mian ji shi :\n");</p><p>　　printf("A[0]=%f,\nA[1]=%f,\nA[2]=%f",A[0],A[1],A[2]);</p><p>　　printf(&qu

125、ot;\nji suan fei yong :\n");</p><p>　　j1=o*D1*(a0+b0*T1)/10000;</p><p>　　for(i=0;i<3;i++)</p><p>　　{if(A[i]<=100)</p><p><b>　　H[i]=1;</b></p&g

126、t;<p>　　else if(A[i]<=200 && A[i]>100)</p><p><b>　　H[i]=1.2;</b></p><p><b>　　else</b></p><p>　　H[i]=1.5; }</p><p>　　for(i=0

127、;i<3;i++)</p><p>　　printf("%f,%f\n",A[i],H[i]);</p><p>　　for(i=0;i<3;i++)</p><p>　　M+=(4400+B-620)*1.2*(0.667+0.0287*A[i])*H[i];</p><p>　　j2=fc*M/10000;