《化工系统工程大作业.docx》由会员分享，可在线阅读，更多相关《化工系统工程大作业.docx（13页珍藏版）》请在三一文库上搜索。
1、第一题一管壳式换热器，单壳程、双管程，用渣油预热原油。渣油的质量流量是Wh68250kg / h ，入口温度为Th382 o C，原油的质量流量是Wc175000入口温度为Tc275o ，求原油的出口温度，换热器的总传kg / hC热系数，传热效率。换热器的结构参数如下：传热面积A=130 m2壳体直径D=0.7 m管子规格25 2.5 mm管中心距P=32 mm，正方形排列管长L=6 m管程通道截面积AT=0.0459 m2壳程通道截面积AS=0.0525 m2渣油侧的污垢热阻R 渣油 0.0005 m2K/W原油侧的污垢热阻R 原油 0.0001 m2K/W解：1初设出口温度 Th2 或
2、Tc2 。1）设定热容流率较小一侧流体的出口温度。 Wh = 68250 /hWc=175000 /h显然：Wh WcWhCph WcCpc即：热流体一侧热容流率较小，所以设热流体一侧出口温度。Th2Th1(Th1 Tc1)0.5382 0.5( 382 275) 328.52）进一步估算冷流体出口温度Tc2 =Tc1+Wh/Wc（Th1-Th2）=275+68250/175000（382-328.5）=295.9( )故初设： Th2=328.5( )Tc2=295.9( )2 . 计算物性参数及传热系数；1）定性温度：Tmc=（Tc2+Tc1)/2= （295.9+275)/2=285.5
3、()Tmh= (Th2+Th1)/2= (382+328.5)/2=355.25()2) 计算物性参数 : 密度、比热、导热系数、粘度。相对密度 DD0.9420.248 X0.174( D20 ) 2 0.08410.312 X0.556 exp( X )式中 X=1+t/100 ; tXD 20D20定性温度 ； A密度：1000DD2020时油品相对密度， D20渣油 0.919 ， D20原油 0.850Xh1T 1004.5525Xc1T 1003.8551D h0 . 9420 . 248 X h0 . 174 ( D20h ) 20 . 08410 . 312X h0 . 556
4、exp(X h )X hD 20 hD 20hDh=0.6866同理 Dc=0.6227 。比热 CpCp0.7072 (0.00147 0.000551D20)t 0.318 D20 (0.055 K 0.35) kcal/(kg oC)式中， K 为油品特性因数，可查取或实测对渣油，K 渣油 12.5对原油，K 原油 12.5Cph=0.7857KcalCpc=0.7499 KcalkgCkgC导热系数 0.100810.00054t / D20oh0 . 1008(10 .00054Tmh) Dc 0 .1008 (10 . 00054 Tmc ) D单位Kcal( mCh )20 h2
5、0 c0 .08864 0 . 1003粘度 (单位为 cSt1cSt=10-6 m2/s)注意区别动力粘度和运动粘度，运动粘度单位 :cStexp exp ab ln(t273)1.22 cst式中：aln ln(11 .22 ) bln( t 1273 )bln(11 .22 )t1273ln/ lnt 2z渣油 ;y原油ln(21 . 22 )273t1501 z15001 y90t 21002 z1202 y132hc1 .61710 60 .6013106SI 单位为 :/s3、据换热器的结构尺寸求i ,0, Ka、安排流程：管程： I，h ，渣油热流体，易结垢，粘度大，管内易清洗，
6、热损失少。壳程： O，C 原油 冷流体。b、分别计算管程和壳程的流速。注意单位： 由 （ /h）换算为 u 的单位（m/s）WhWhu0Wcui0 .6021.487A ih A hh 3600cA0N单位是： m/sc、计算Re，Pr 准数。Red i u ihd i i7455 .66ihhRed 0 e u 0c. 90c67285PrC phhC phhhi41 . 1193hhPrC pccC pc u cc11 . 690cc其中： di=0.025-20.0025=0.024 ( p 2d 0 2 )4d eod 04 ( 0 .032 20 .025 2 )40 .027180
7、 .025d、 计算0 ， i ；ai0.023Re0.8 Pr0.33(u / uw)0.14da00.36Re0.55 Pr0.33(u / uw)0.14de30.527Re 0 0 .5cPr0de 0Re 0 0 .550 .36hRe i 0.8Pri0 .023di0 0 .33i 0 . 33Re20001352 .365 .Re2000435 .48e、进行壁温校正：Qc=WcCpc（Tc2-Tc1）Qh=WhCph（Th1-Th2）据：冷侧： TwcTmcTwc=（ Tc1+Tc2） /2+Qc/（0A0 ）热侧： Twh TmhTwh=（Th1+Th2）/2Qh/（i A
8、i ）Qc = WcCpc(Tc2 Tc1)= 3.18 106QH =WhCph(Th1 Th2)= 3.3 106只估算一侧壁温，认为Twc =Twhf、求 K 0：1K0=R0 + 1/ 0+ (d0/ di) Ri + (d0/ di)= 215.41A ond o LA ind i LJ/sJ/s1/ 0g、求= 0.3746NTU 1 (Cmin Cmax)2KANTU( WCp) min2e ()1(C min C max)1(C min C max)2 1(1e ()4、 计算 TH2 ：T H2 = T H1 -（ T H1 T C1）=382-0.3746(382-275)
9、=341.90CT - TH2 0.1H2所以重新计算 TC2、 Tmc、Tmh、物性参数等等。第一次计算数据：4Th1Th2TmhD20hXhDhCph382328.5355.250.9194.55250.6866430.785686Tc1Tc2TmcD20cXcDcCpc275295.9285.450.853.85450.6226930.749927 ha zbzh hAtWh0.08864318.917-2.92961.616731.6167*10-60.045968250 caybyc cAsWc0.10030922.81-3.68710.6013640.6013*10-60.0525
10、175000UidiReiPriiCmax0.6015290.027441.29841.19797435.0817131237.2U0de0Re0Pr 00KCmin1.4869690.027267256.3911.721811352.365215.271653623.08NTUTh20.5218890.5637730.374406341.9386第一次迭代过程及结果：Th1Th2TmhD20hXhDhCph382341.9386361.96930.9194.6196930.6805840.792404Tc1Tc2TmcD20cXcDcCpc275290.624282.8120.853.828
11、120.6256960.747185 hazbzh hAtWh0.08824518.917-2.92961.5274071.5275*10-60.045968250 caybyc cAsWc0.10047822.81-3.68710.6206990.6207*10-60.0525175000UidiReiPriiCmax0.6068840.027946.59739.08365448.6338130757.4U0de0Re0Pr 00KCmin1.4798330.027264848.612.092231367.653220.13854081.57NTUTh20.5291620.5726380.3
12、77652341.59125第二次迭代及结果：Th1Th2TmhD20hXhDhCph382341.5912361.80.9194.6180.6807360.792235Tc1Tc2TmcD20cXcDcCpc275290.756282.8780.853.828780.6256210.747254 hazbzh hAtWh0.08825518.917-2.92961.5295781.5296*10-60.045968250 caybyc cAsWc0.10047322.81-3.68710.6202070.620210-60.0525175000UidiReiPriiCmax0.6067480
13、.027933.53739.13517448.2896130769.4U0de0Re0Pr 00KCmin1.4800110.027264907.812.082821368.242220.040554070.02NTUTh20.5290410.5724810.377602341.5966由第二次可以看出Th2 Th2 = 0.0054 0.1结束计算。故原油的出口温度为 341.597 0 C，换热器的总传热系数为 220.041 W/（ 0C），传热效率是 0.3776.6第二次大作业用黄金分割法求下面问题的最优解，minf （x）=e xx 2，s.t2x3精度要求，将原始区间缩短5 倍。
14、解： a=-2 b=3x1b（ ba）30.618034 （32)0.09017x 2a（ ba）20.618034 （32） 1.09017f （ x1）1.10249f （ x 2） 1.5246f（ x1） 0.9219 f （ x 2） 1.5246a=-2 b=1.09107x1b（ b - a） 1.09017 0.618034 （1.090172） 0.81996x2a（ b - a)20.618034 (1.09017 2)0.09017f （ x1） 1.1124f（ x 2） 1.10249f （ x 20.09017） 1.10249 f（ x10.81966） 1.11
15、24a=-0.81966 b=1.09017x10.09017f（ -0.09017） =1.10249x2 a （ b - a） 0.36067f （ x 2） 0.82729f （ x 2） f （0.36067）0.82729 f （ x 1） f（ 0.09017）1.1.249舍去 a. x1 a=-0.09017 b=1.09017x1 0.36037 f （ x 1） 0.82729 x 2 a（ b- a） 0.63932f （ x 2） 0.93638f （ x1） f（0.36067） 0.82729 f（ x 2） f（0.63932） 0.93638舍去 x 2 .b a=-0.09017b=0.63932若要求将原始区间缩短5 倍。b-a=0.63932+0.9017=0.73102 l n = loEn =50.2=1x*x1 0.36067 f *f （0.36067） 0.8272978